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Part I 

On the QCD coupling at low 
energy scales 



Chapter 1 

The QCD coupling in 
Perturbation Theory 



The Renormalization Group method is an inherent part of the theoretical 
description of strong interaction processes. As known, the values of the 
QCD coupling extracted from high-energy data are in good agreement with 
the theoretical running coupling a s (Q 2 ) as derived by the RG equations and 
properly normalized at the Z boson mass, thus providing a powerful test of 
asymptotic freedom [J. 

On the other hand, in the low-energy domain, the running coupling a s (Q 2 ) 
develops, at any loop level, unphysical (so-called Landau) singularities that 
contradict the general principles of the local Quantum Field Theory and 
severely complicate the theoretical analysis of hadron dynamics. 
A reasonable prescription to get rid of the Landau singularities in the QCD 
coupling is given by the Analytic Perturbation Theory (APT) approach. In 
this framework, the analycity properties are recovered by imposing a s (Q 2 ) to 
satisfy a dispersion relation with only the unitary cut on the time-like axis, 
while preserving the asymptotic freedom constraint. 

This prescription has been exploited in the framework of a second order 
Bethe-Salpeter-like (BS) formalism adjusted for QCD in recent works [21], for 
the calculation of the meson spectrum in the light and heavy quark sectors. 
It has been there enphasized the relevance of the infrared evolution of the 
APT coupling in the few hundreds MeV region for a reasonable reproduction 

x For a summary of the high energy mesurements of a s (Q 2 ) see [TJ [2]. 



of ground state mesons involving light and strange quarks, at variance with 

a naive truncation prescription. This feature strongly supports the APT 

approach within the BS framework. 

On the other hand, in the present work [3], 0] , the calculations of the meson 

spectrum within the BS formalism have been exploited in order to obtain an 

experimental determination of the strong coupling a° xp (Q 2 ) below 1 GeV, by 

comparison with the data. 

The results turn out to be in satisfying agreeement with the APT coupling 

and its evolution in the 200-500 MeV region. Furthermore, below this scale, 

some hints on the deep IR behavior are further provided, albeit affected by 

large theoretical errors. 

In Sec. 1 the Renormalization Group method in QCD is briefly reviewed, and 
the state-of-the-art on the infrared behavior of the QCD-invariant coupling 
sketched, whereas in Sec. 2 the basis of the Analytic Perturbation Theory 
approach are outlined. Sec. 3 is devoted to the determination of the exper- 
imental values of the QCD coupling from the meson spectrum, the analysis 
method and the discussion of the results. 

Finally, some technical material is exposed in the Appendices. A short 
description of the derivation of the second order BS formalism is given in 
App. A, whereas numerical tables in App. B display all the results in detail. 



1.1 The Renormalization Group method in 
QCD 

In the framework of Quantum Field Theory, an arbitrary change in the renor- 
malization prescription basically results in a reparametrization of the theory 
by a suitable redefinition of couplings and masses, so as to leave physical 
quantities unaffected. The behavior of relevant Green functions under a 
change in the extra mass scale /i introduced by renormalization is specified 
by a set of continuous transformations, properly expressed by the differential 
equations of the renormalization group (RG) q This leads to the notion of 
effective couplings and masses. 

Choosing a renormalization scheme (RS) which performs all subtractions off 
shell, mass parameters can be treated on the same foot as couplings. 



2 For an overview on history and applications see for instance [5] and references therein. 



In the "minimal subtraction scheme" (MS) the renormalization factors ex- 
hibit no mass-dependence and are of the general form 

oo 

Z 3 ^a s ) = l + Y,e- n zf\a s ), (1.1) 

n=l 

where a s = g 2 /47r is the strong coupling and e is the deviation from the 
physical space-time dimension, d = 4 — e . The so-called "modified minimal 
substraction scheme" (MS) amounts to subtracting the whole (ln47r — 7_e)- 
term together with the poles, and it has become a common standard in QCD. 
In MS-like schemes the arbitrary mass scale /i arises to keep the coupling 
dimensionless in 4 — e dimensions, and the renormalized objects turn out to 
be scale dependent. 

In transition from one parametrization (/i 2 ,, a s (/io), ^(/^o)) t° a different one 
(/i 2 , a s (fi 2 ),m(fi 2 )) Green functions are connected by finite renormalization 
transformations, and, in order to have a unique theory under an infinitesimal 
shift of the scale, they must satisfy the (inhomogeneous) RG equation 

RTip^a^rrifj,,/! 2 ) = j T (a»)r(pi] a M , m M , /i 2 ), (1.2) 

where Pi are relevant momenta and the shorthand a^ stands for a s (fi 2 ) (sim- 
ilarly for rri/j). The RG operator is the total derivative over the scale param- 
eter, namely 

R "" 2 ^ + ^ + 7 '" m ^- (L3) 

where the sum over nj quark flavors in the mass-term is understood, and the 
universal functions 



!<9«m ( , _ /i 2 dm^ _ 2 d 



/?((*,,) =fi -p-f, 7m(a M ) = 7TT and Tr(au) = // — ^lnZ r (1.4) 



have been defined. The RG coefficients /3, 7 m and 7r are computed in per- 
turbation theory (PT) with the aid of the relative counterterms. Generally 
speaking, they are gauge-dependent finite functions of couplings and masses, 
once the physical limit e — > is performed; actually, in any mass-independent 
RS (such as minimal subtraction prescriptions or trivially in massless theo- 
ries), the RG coefficients are mass-independent and function of a s alone. 

As far as the QCD coupling is concerned, the beta function in the first of 
Eqs. (11.41) has the usual PT power expansion 

(3(a s ) = -{3 a 2 s - (3 ia 3 s - f3 2 a A s + 0(a 5 s ) , (1.5) 
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where Pi are the 1-loop coefficients. At 1-loop level (i.e. keeping only the first 
term in ( 11.51) ) solution to the first of Eqs. (11.41) reads 



a s (n 2 



a«(/4J) 



1 +/5 a s (/ig)ln(/i 2 /^ 



a 



.(^)E 



ra=0 



-p a s (nl)\n — 
Mo 



(1.6) 
namely, while the theory does not predict the actual size of a s , its scale 
evolution is completely known. Indeed, a s is not an observable by itself, but 
plays the role of the effective expansion parameter for physical quantities, 
and any determination of the strong coupling is then extracted from the 
measurements of observables for which PT predictions exist. Conventionally 
the mass of the Z boson, Mz ~ 91.2 GeV, is used as a reference scale 
to the world average of all (high energy) determinations, which is currently 
a s (Mf) = 0.1176(20) [21 E]. 
Defining the overall scale parameter 



A 2 = fj,Q exp 



A) «*(/■*§) 



Eq. I 11.6P becomes 



a s (n 2 



1 

/5 ln(/i 2 /A 2 ) 



;i-7) 



'1. 



Since 47r/3o = (11 — 2n//3) is a positive number for all presently known 
flavours, a s (/x 2 ) decreases to zero as /i 2 — > (asymptotic freedom), this 
property being not spoiled by higher loop corrections. 

Given the arbitrariness of the renormalization scale \i , the best choice is a 
matter of expediency. The commonly adopted prescription is to fix \x at the 
physical scale of the process at hand. This can be motivated by a glance to 
the homogeneous RG equation satisfied by any (dimensionless) observable 
G, namely |f| 



/'" 



_d_ 
1^ 



G 



l<- 



_d_ 



+ P(a s 



d 

— 

da. 



G = 



;i.9) 



3 Actually quark masses have been neglected, assuming that both the physical and the 
renormalization scales Q 2 and /i 2 are larger than any other relevant mass scale. If non 

(El 



negligible, quark masses give rise to an additional term in Eq. (|1.9|) . similarly to Eqs. 
and (fOj) . 
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where G is typically function of the ratio u Q 2 /^ 2 and of the renormalized 
coupling a s (fi 2 ) , and is given in PT by a power series in a s (fj, 2 ) 

-•(SMS)^(S)(^)Watf> 

Eq. (II. 9p clearly states the scale invariance of any physical observable, that 
is fulfilled only in infinite order. In any truncated order, however, the cancel- 
lation of the scale among the coupling and the coefficients is not complete, 
that is, a residual scale and scheme dependence still persists, and the degree 
of cancellation improves with the inclusion of higher orders. 
The solution of Eq. ( jl.9f) for the PT coefficients gj exhibits a logarithmic 
dependence on Q 2 /fi 2 (starting from the order a 2 ), that can be controlled 
by setting Q 2 = fi 2 . Then Eq. ( 11.101) takes the form 

C(»)=J 0+ 5^W^>Y + O(a;), (1.11) 

7T V 7T I 

and the energy dependence of G translates into the energy dependence of the 
effective parameter a s (Q 2 ) . A change in the value of fi 2 clearly amounts to 
a reorganization of the PT expansion. 

The RG improved perturbative QCD yields a consistent picture of high en- 
ergy strong interaction processes U from a few GeV up to a few hundred GeV. 
On the other hand, the existence of infrared singularities in the RG-invariant 
coupling a s (Q 2 ) spoils the theoretical analysis of low energy hadron dynam- 
ics. Indeed, at 1-loop level a s (Q 2 ), as given by Eq. (11.81) . has a pole-type 
singularity at Q 2 = A 2 , and the effect of including higher order contributions 
does not overcome the problem and only amounts to modifying the singular- 
ity structure. On the contrary, physical observables are known to be analytic 
in the entire q 2 complex plane aside from a cut on the real positive (time-like) 
axis, and their exact expressions should be free of these singularities. 
Many strategies have been suggested to get rid of unphysical singularities 
within PT calculations, so as to handle the strong interaction processes at 



4 Actually, according to the process, more than one scale can be involved, and this can 
be spacelike (q 2 < and Q 2 = —q 2 ) or timclikc (q 2 > and s = q 2 ). 

5 Some care is required for timelike observables since the PT coefficients, in this case, are 
modified by non negligible corrections, proportional to powers of /3qtt , due to the analytic 
continuation fron the space-like to the time-like domain (see Sec. 1.6). 
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low energies. Some of such methods, reviewed e.g. in Ref. [6], originate 
in the general properties of the perturbative power series for the QCD ob- 
servables in the framework of the RG formalism, such as the arbitrariness 
of the renormalization scheme and scale. One can remind, for instance, the 
Brodsky-Lepage-Mackenzie criterion [7J, the "Optimal conformal mapping 
method" [5] (see also Ref. [5]), and the so-called "Fastest Apparent Con- 
vergence" technique [TO]. The latter, e.g. amounts to setting the scale so 
that all the higher (than LO) coefficients vanish. It is related to the effective 
charges, [TTj-[T3], i.e. alternative definitions of the QCD coupling, straight- 
forwardly related to a physical observable, either of spacelike or time-like 
argument. 

In this context it is also worth mentioning the "Optimized Perturbation 
Theory" approach [TU [15] , with the aim of improving the convergence of PT 
expansions by a proper choice of the scale and of RS parameters (such as j3 2 
and fa), on the basis of a minimum sensitivity criterion. 
However, all these definitions can be related to each other by referring back 
to the conventional MS. Most of these techniques yield theoretical predic- 
tions on the low energy behavior of the RG-invariant coupling which are in 
qualitative agreement, suggesting in all cases an IR finite limit. 
Furthermore, recent results of lattice simulations testify to the absence of 
spurious IR singularities in the QCD coupling [T6| [T7j. 

The ghost-pole issue gives rise to severe complications, in particular, as far 
as the bound states problem is concerned, since the characteristic scale Q 
involved (i.e., the momentum transfer in the quark-antiquark interaction for 
mesons) is typically below 1 GeV, according to the state and the mass of 
quarks involved. 

Among the many potential inspired attempts to modify the expression of a s 
in the low energy domain, there exist either trivial tricks as the trucation 
prescription u, as well as more sophisticated models. For instance, one can 
remind the 1-loop definition of the QCD coupling in connection with the 
momentum representation of the Cornell-like potential [18] 



^- am+qvAv *<«=-2?Hr' (L12) 

as Q — ► oo , the usual 1-loop coupling ( 11 .81) is recovered from the first of Eqs. 



6 This simply amounts to assuming a s (Q 2 ) frozen at a maximum finite value in the IR 
region. 
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( 11.121) . whereas as Q — > the latter diverges as A 2 //3oQ 2 , incorporating non- 
perturbative effects Q Estimation of higher order corrections and further 
developments along this line can be found in [19]. On the other hand, a 
number of analyses seems rather consistent with an IR finite coupling and a 
suitable modification of the quark- ant iquark potential to encode confinement 
(see e.g. [211 ED [221 EZj). 

A reliable algorithm to get rid of Landau singularities is provided by the 
Analytic Perturbation Theory (APT) approach [23] , which represents a next 
step, after RG summation, in improving the perturbative results. The APT 
method, as discussed to some extent in Sec. 2, is based on the causality con- 
dition which imposes the analyticity constraint on the RG-invariant coupling 
on the whole cut complex plane C — {q 2 > 0} . PT is then used to evaluate 
the spectral function, i.e. the discontinuity of a s (Q 2 ) across the cut. At 
1-loop, e.g., it reads 



lm[a s (—a)\ = \ , (1-14) 

[1 + (3 a s (/j,l) ln(a///§)] + [vr/5 a s (/ig)] 2 



where Eq. ( 11.61) has been used. Indeed, within the causal approach, Landau 
singularities are suppressed at their very roots, and the IR-safe QCD coupling 
possesses a universal finite limit as Q 2 — > . 

In what follows the basis of the derivation of the QCD coupling from the 
Renormalization Group equations are reviewed, both in the space-like and 
time-like domains. 



7 This can be seen by comparing the second of (|1.12|) with the usual Cornell potential 
in the momentum space representation 



1 4q s (/i 2 ) 1 a 



nQ)=~^-~, d.13) 



and identifying the quantity 2A 2 /3/?o with the string tension a that controls the confining 
part of the potential. 
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1.2 QCD (3 function 



The RG equation for the effective coupling is actually solved in PT with the 

aid of the relation 

lid 
0(a s ) = -a s —\nZ a (1.15) 

where the limit e — > is understood. Thus one needs to compute the renor- 
malization factor for the coupling Z a and this can be accomplished in sev- 
eral ways. One can start from the quark-gluon vertex Zq qg together with 
renormalization factors of quark and gluon propagators Z q and Z g to obtain 
Z a = Z^ qg Z~ 2 Z~ x , but other choices, as ghost-ghost-gluon vertex or trilinear 
and quartic gluon interactions, equally work. In a minimal subtraction pre- 
scription the four-dimensional /3-function is entirely specified by the residuum 
of the simple pole in the e-expansion of Z a { 1 1.1 1) 



PM 



l a 2 A z w ( 

2 "da. a 



oc. 



(1.16) 



Then the MS-coefficients in Eq. (jl.5p reads [29] 



* = i^ 



0i 

02 
03 



(An) 2 

1 
(4tt) 3 

1 
(4tt) 4 



11 - J", 



38 

102 n f 

3 ; 



2857 5033 



2 U 

149753 



-n f 



325 , 

54 f 



(3 



3564C; 



1078361 6508 



162 



27 



Cs) 



"/ 



/ 50065 6472 



V 162 



81 



a) 



n 



./' 



1093 3 
"729"^ 



(1.17) 



where ( u is the Riemann zeta- function, £3 ~ 1.202057. The coefficients 
0j generally depend on the RS, whereas the first two are universal within 
massless schemes. Moreover, in the MS-scheme the /3-function is gauge- 
independent in any order [30], and in an arbitrary mass- dependent scheme 
this feature is preserved only in the first order. 

As already noted, the universal 1-loop coefficient [31] has a positive sign pro- 
vided there is a small enough number of quark fields (rif < 33/2); thus the 
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theory is asymptotically free, that is, the /3-function has a stable UV fixed 
point as its argument decreases to zero. Indeed this coefficient is the sum 
of two contributions, the relevant one with respect to asymptotic freedom 
property being the first, which arises from pure gauge field effects i.e. from 
the nonlinear Yang-Mills interaction terms. The 2-loop coefficient has been 
computed for the first time in [32] and is positive up to rif = 8 . 
Higher order approximations are scheme- sensitive, and it is common practice 
to perform computation within minimal subtraction prescriptions, in which 
the /^-function is unchanged. The first calculation of 3-loop coefficient is 
due to [33], where the ghost-ghost-gluon vertex in the Feynman gauge was 
used, and it turns out to be negative for 6 < rif < 40 . In a more recent 
work [M] the quark-gluon vertex was used instead, providing an indepen- 
dent check in an arbitrary covariant gauge of the previous result and its 
gauge-independence. Finally, the original 4-loop calculation [29] has been 
performed using the ghost-ghost-gluon vertex in a arbitrary covariant gauge, 
and for a generic semi-simple compact Lie symmetry group. The result turns 
out to be gauge- independent, and involves higher order group invariants such 
as quartic Casimir operators; specialized to the standard SU(3) symmetry, 
the 4-loop coefficient is a positive number for every positive rif (see also 
[35]). Finally note that (except for /3 2 ) all four coefficients are positive up to 
n/ = 6. 



1.3 RG-invariant coupling 

The formal solution of the first of Eq. (11.41) is easily worked out 

/i 2 [ a ^ 2) da* 



lnC 2=/ 7^> (1-18 

and yields the evolution of the effective coupling as a function of two dimen- 
sionless arguments t = /i 2 /// 2 , and a = a s ({J%), where /Xq can be viewed as 
a fixed reference scale and fi 2 as a sliding one. It is worth noting that by 
differentiating Eq. f II . 1 8 j) with respect to «o one has 

t— -p(a o ) — )a s (t,a o )=0, (1.19) 

namely the homogeneous RG equation satisfied by the effective coupling (or 
invariant charge). 
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The exact 1-loop solution (11. 6ft or ( 1 1.8ft is obtained by retaining only the first 
term in ( 11.51) . In Eq. ( 11.81) . in particular, the dimensional scale A keeps track 
of the initial parametrization (/i , a s (/io)), and it is scale-invariant; its value 
is not predicted by the theory but must be extracted from a measurement 
of a s at a given reference scale. Emergence of a scale parameter, sometimes 
referred to as dimensional transmutation, breaks naive scale invariance of 
the massless theory, and it is commonly believed to be associated with the 
typical hadron size i.e. to the energy scale where confinement effects set 
in. Roughly speaking, A is the scale at which the (1-loop) coupling diverges 
(Landau pole), and perturbation theory becomes meaningless. Further, A is 
scheme- dependent and receives further corrections at each loop level (but for 
simplicity the same notation is used throughout). 

At the 2-loop level the integration of (jl.ISp leads to a transcendental equa- 
tion, that is, starting from the 2-loop approximation to the /3-function in 
(II. 5ft . one has 

In !L = C + — *— + B 1 In a s - B x In ( 1 + ^-a s ) (1.20) 



Ho Po^s \ A) 

where B\ = (3\/ 0$ and the constant term C from the lower end points can 
be again reabsorbed into the A-parametrization, with the commonly adopted 
prescription (see e.g. [3"B"l 157] ) 

A 2 
ln^^C-Sxln/^o, (1.21) 

K 

which fixes a specific choice for A . Thus one gets the 2-loop implicit solution 

In £ = ^- - Bt In (b 1 + -±-) , (1.22) 

A 2 (3 a s V Po« s / 

from which the 2-loop scaling parameter is immediately read. To achieve an 
explicit expression for the running coupling at this level one should resort to 
the many-valued Lambert function W(() implicitly defined by the equation 

W(()ex V {W(()}=C. (1-23) 

The function W(() has an infinite number of branches Wk(£) k = 0, ±1, ±2 . . . 
such that W*(C) = W_ n ((*) (for more details see [SB])- The exact solution 
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to Eq. (fl~22D for < n f < 8 reads (see e.g. [39] 



l/Bi 



^w-ftkirrbo c -^(iJ (L24) 

where z = fi 2 /A 2 , and W_i(C) is the "physical" branch of the Lambert 
function, since it defines a real values function for ( e (— e _1 , 0) which fulfills 
the asymptotic freedom constraint. As e — > + the asymptotic relations 

W_i(-e)=ln£ + 0(ln|lne|) (1.25) 

W_i (-- + ej =-1- V2ee + 0(e) (1.26) 

hold. Outside this region of the real axis W-i(-e) takes on complex values. 
Though not easy for practical aims, Eq. (11.241) yields the most accurate ex- 
pression for investigating the IR behavior of the running coupling, since it 
has not been derived by means of deep perturbative approximations (aside 
from the truncation of the 2-loop f3 function). 

Actually, a frequently used 2-loop approximate solution, known as the itera- 
tive solution [23], is obtained starting from Eq. (11.201) together with a single 
iteration of the 1-loop formula (II. 8p . that is 



lt V ; ln^ + 5 1 ln(l + 5f l1 - 
where z = /i 2 /A 2 , and A now defined in (11.201) by 



ln^C-^ln^. (1.28) 

This definition is related to the standard one ( 11.211) by the factor 

\n(A/A) = -B 1 \nB 1 . (1.29) 

However, the commonly used 2-loop solution is an asymptotic formula which 
strictly relies on the smallness of a s for fairly large /i 2 , since it amounts to 



8 Note that if 9 < n/ < 16 the principal branch Wq is involved, but here and throughout 
the discussion is focused on n/ < 6 . 
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solving Eq. (11.201) (with the choice (ll.2ip ) where the last term in the r.h.s. 
has been neglected. Again, after one iteration of the 1-loop formula, the 
result is then re-expanded in powers of 1/L, where L = In z and z = /i 2 /A 2 
as before 



a<?Hz) 



Polio, z 



ft In (In z) 

]3l \nz 



1.30) 



Eq. fll.3Up is known as the standard 2-loop running coupling and works only 
in the deep UV regime, i.e. for L ^> 1 . In Fig. 1(a) the fractional differences 
of the two approximate formulas Eqs. fl 1 . 2 T[) and (11.301) with respect to the 
exact 2-loop coupling Eq. fjl.24|) are displayed. 

Under the same assumptions of Eq. (11.301) one can easily derive the 3 and 
4-loop asymptotic formulas. Starting with the approximate implicit solution 
at the 4-loop level, namely 



In— ^ 
Mo 



C+ 



.Pi, 

+ — In a x 



P a s % 



P2P0-PI , /?! 3 - 2/?p/3 1 /? 2 + /?g/3 3 2 ^ n( 3, 
-«s+ -^ttta a a +U{a 8 ), 



0S 



2/5 4 



(1.31) 

the 4-loop coupling in the standard form of an expansion in inverse powers 
of the logarithm L for L ~^> 1 (see e.g. [10] ) reads 



a 



(i) 



(^ 



$L» 



P L 






ffiL 2 



%<-" 



L + -\n 2 
2 



L + 21nL-- 



01 

1 

2 



(In 2 L - In L - l) 



Po 



M, TJ Ps' 



Pi 



1.32) 



The 4-loop solution (ll.32p turns out to be nearly indistinguishable from the 
3-loop curve (see Fig. 1(b)). In Eq. (11.321) the 1-loop solution (jl.8p has been 
emphasized, and being the leading UV term in (11.321) . it defines the asymp- 
totic behavior of the RG-invariant coupling. On the other hand, the 2 and 
3-loop asymptotic expressions are easily read from Eq. (I1.32p by keeping only 
the first two and three terms respectively. 

Exact integration of the truncated 4 or 3-loop /3-function leads to a more 
involved structure than Eq. (I1.3ip . Nevertheless, in [39] a useful solution has 
been still worked out at 3-loop level via the real branch W-i(C) of the Lam- 
bert function together with the Pade' Approximant of the related /3-function. 
Moreover, in [41] a reliable approximation to the higher-loop coupling has 
been suggested, via a power expansion in the 2-loop exact coupling (11.241) . 
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Figure 1.1: (a) Fractional difference between Eqs. (|1.24j) . (|1.30[) and (|1.27|) (solid and 
dashed line respectively), with A = 350 MeV and rif = 4. (b) From pQ: fractional 
difference between the 4-loop and the 1-, 2-, 3-loop approximations to Eq. (|1.32[) . with 
Tif = 5 and normalization conditions at a s (A/§) = 0.119 . 



of the form 



a 



(*) 



(S 



5>i fe) 



[«£V)]' 



1.33) 



n>l 



with k > 3 the loop order, and pn proper functions of the coefficients (3j. 
Comparison [12] of these multi-loop approximants with the relative asymp- 
totic formulas, Eq. (11. 32ft and the 3-loop analogue, reveals the better agree- 
ment of the former ones with the higher-loop exact coupling numerically 
estimated (i.e. starting from the exact implicit solution), even at rather low 
energy scales (see also [12]). 



Finally, as far as the definition of the scaling constant A is concerned, few 
comments are in order. As already pointed out, starting from the 2-loop 
level an arbitrary constant has to be fixed for A to be uniquely defined; be- 
side the commonly accepted convention (11.211) (or (11.28P ). other prescriptions 
have been proposed [H] in order to optimize the 1/L-expansion, while (ll.2ip 
does remain the preferred one as no further terms of order 1/L 2 appear in 
the 2-loop asymptotic formula (11.301) . Thus, at higher-loop levels the scal- 
ing parameter is analogously related to the initial parameterization, and the 
4-loop (and 3-loop) value reproducing the world average a s (M|) = 0.119 is 
A^ =5) = 220 MeV m. 
Furthermore, with regard to the scheme-dependence of the coupling and the 
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scale parameter, restricting the discussion to mass-independent RS (as MS- 
like schemes or trivially any prescription in the massless theory), the renor- 
malized couplings in two such different schemes can be related perturbatively 
at any fixed scale \i 



\JL a \JC a 



1 +V 1 h V 2 

7T \ TV 



1.34) 



Then, it can be verified that the first two coefficients of the relative (3- 
functions do not change as the renormalization prescription is changed, while, 
for instance, the 3-loop ones are related by (3' 2 = /3 2 — f i/3i + v 2 Po — Povf . As a 
result, the runing coupling at each loop-level (e.g. Eq. (11.321) ) depends on the 
RS, through the coefficients /3j with j > 2 and the initial parameterization. 
The latter obviously amounts to suitably adjusting the A , and the relation 
is given by [55] 

In (A/A) = J* (1.35) 

ZirfJQ 

which works through all orders (e.g. Aj^/Ams = exp [(ln47r — 7 e )/2] ~ 2.66). 



1.4 Threshold matching 

Quark mass effects reveal themselves in explicit corrections within higher 
order perturbation theory (see e.g. [IS]), and in the energy dependence of 
the effective (running) quark masses as a result of the RG improvement [57]. 
A direct effect of the quark masses on the evolution of the coupling is through 
the dependence of the /3 coefficients on the number n/ of active quarks, i.e. 
satisfying m/ <S //, where \x is the renormalization scale and m/ the MS 
quark mass (the definition can be also formulated in terms of the pole mass 
Mf that can be related to the former [4*8]). 

Then, in order to have a continuous coupling on the whole domain, the two 
functions a s Uf (yU 2 ) and a s n/ (/i 2 ) must be matched at each heavy quark 
mass scale fx^f^ = 0(m,f) . As a result the scaling parameter A depends also 
on rif (see e.g. [491150]). 

The most straightforward way is to impose continuity of a s by means of the 
matching condition a s nf ( m }) — o; s nf (ml), which works up to next-to- 
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leading order. At 1-loop, it yields 



A (n f ) = A (n / _1) 



A ( n/ -1)-| 2/(33-2^) 

m f 



(1.36) 



that can be extended up to 2-loop, and exhibits explicit dependence on rrif . 
Since trivial matching does not generally hold in higher orders within MS 
scheme, a more accurate formula is required in this case; specifically to obtain 
the global evolution of the 4-loop coupling the proper matching condition 
reads 



(nt-l) (nt) 

a. s = a s 



l + k 2 



a 



(»/)' 



+ h 



a 



(»/)' 



7T 



7T 



1.37) 



with 
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564731 82043 



124416 27648 



Ca 



2633 
31104' 



\n. 



1) 



'1.381 



if /x 



(«/) 



rrif is exactly assumed. With this convention, Eq. ( 11.371) yields the 



relationship between A^ n/ ^ and A*"^ in the MS scheme. One can equally 
fix /i( n/ ) = Mf , that amounts to a proper adjustment of the coefficients in 
ffTBTI) ; for instance [1] starting with (fl~32l and A( n /= 5 ) = 220 MeV, the values 
A („ /=4 ) _ 305 MeV and A (n /= 3) _ 346 Me y are obtained, with thresholds 

fixed at the pole masses Mb = 4.7 GeV and M c = 1.5 GeV. 
Note finally that Eq. (11.371) spoils continuity of a s ; therefore, one can resort 
to a more sophisticated technique [HD] , which relies upon mass- dependent RS 
and yields a smooth transition across thresholds. 



1.5 Landau singularities 



As already noted, the 1-loop coupling (11.81) is affected by a spacelike pole at 
A with residue l//?o • Adding multi-loop corrections does not overcome the 
hurdle; rather, the singularity structure is more involved due to the log-of-log 
dependence, so that a branch cut adds on to the 1-loop single pole in the IR 
domain of the spacelike axis. 
Moreover, at a given loop level, Landau singularities sensibly depend on the 
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approximation used. For instance, the 2-loop iterative formula ( 11.271) has 
a pole at z = 1 (/i = A, see Eqs. (PL281) and (EHD ) with residue l/(2/3 ) , 
and a cut for < z < exp (— -Bi) due to the double logarithm. On the other 
hand, in the 2-loop asymptotic expression Eq. (I1.30p . the singularity in z — 1 
becomes stronger, namely 

af\z)^-^^ z = ^/A^l (1.39) 

and the cut now runs from to 1 . Analogously, the 3- and 4- loop asymptotic 
solutions, as given by Eq. (11.321) . behave as 

, 3) , , B 2 \n 2 (z-1) , 4) , . fl?ln 3 (*-l) . .. , n , 

Po [z - ly po {z - 1) 4 

and are equally affected by an unphysical cut. However, the cumbersome sin- 
gularities structure of the leading Landau singularity, and of the unphysical 
cut as well, are an artifact of the UV approximations introduced. Therefore 
to deal with the IR behavior of the invariant coupling it is necessary to face 
with the exact solution, e.g. Eq. (11.241) at 2-loop; clearly it is singular when 
W_i(C) = -1 that is at z = B~ Bl (^ 2 = B~ Bl A 2 ), corresponding to the 
branch point £ = — 1/e of the Lambert function. In the neighborhood of this 
point, due to the asymptotic behavior (11.261) of W_i(Q , one has 



g(*)*^^ L [z-B;*]- 1 ' 2 1 (1.41) 



i.e., an integrable singularity (the factor B 1 1 in front of A in the singular 
point can be reabsorbed into a proper redefinition of the integration constant, 
through Eq. (0511 ). 

A more detailed investigation about the IR singularity structure of higher- 
order approximations of the QCD coupling has been performed in Ref. [43] . 
on the basis of Eq. (1 1 . 3 3 j) . where, in particular, the location of Landau sin- 
gularities as a function of n/ are studied. 



1.6 Time-like coupling 

The evolution of the QCD coupling has been derived above in the space-like 
region. However, in order to parametrize observables of time-like arguments, 
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one needs an effective coupling defined on the whole real axis. 
While this poses no special problems in high energy processes, at any finite 
energy the issue of which should be the most suitable parameter in the s- 
channel must be carefully considered. The standard practice is to merely 
take over to the time-like domain the space-like coupling at any loop level, 
regardless of the crossing between two disconnected regions, thus exporting 
the IR singular structure in the s-channel. Nevertheless, from many early 
works based upon analysis of e + e~-annihilation into hadrons (see [51] and 
references therein), it is known that this should not be the case except far in 
the asymptotic regime. This is because of the appearance of not negligible 
corrections (7r 2 -terms) in the higher order coefficients of the a s (s)-expansions, 
due to analytic continuation from space-like to time-like axis. The problem 
is related to the IR non analyticity of the invariant coupling, and it receives 
a satisfactory answer in the framework of Analytic Perturbation Theory, in 
which the space-like and time-like channels are treated in a unified description 
(see Sec. 2). 
Referring for definiteness to the e + e~-annihilation ratio into hadrons 

R(s) = " ( f f - ha ^°" S) , (1-42) 

where s = q 2 > , one can start with the hadronic contribution to the photon 
polarization tensor in (space-like) momentum space 

nf (<?) = {g» v q 2 -q»q v )U h {-q 2 ) oc j ' d*xe* x < 0|T(^(x).f (0))|0 > (1.43) 

j M being the quark electromagnetic current operator. As is known, exploiting 
the analytical properties of the two-point correlation function (11.431) in the 
cut complex plane C — {q 2 > 0} , R(s) can be straightforwardly related to 
the discontinuity of Hh(— q 2 ) across the cut 

R(s) = — lim [U h (-s + ie) - U h (-s - ie)\ . (1.44) 

In order to compute II/ l (— q 2 ) on the space-like axis (q 2 < 0) one formally 
works with its first logarithmic derivative (thus avoiding subtraction con- 
stants), the Adler D-function [52] 

D( V) = -^tfl (1 .45) 
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The RG improved PT expansion of the D-function reads 



D PT (Q 2 ) = 3 J2® 2 



f 
f 



^m +d j^m)\ d j<Q^ 



TX \ TX I \ TX 



(1.46) 
where Q 2 = —q 2 , Qf are the quark charges, and the coeffcients in the MS 
scheme read [53] 

d 2 ~ 1.986-0.11571/ 

d 3 ~ 18.244 - 4.216/2/ + 0.08Qn 2 f - 1.24 ( 3 J^ Q) J ( $3 Q/ ) (- 1 - 47 ) 

Then, by integration, Il/ l (Q 2 ) is readily obtained q The result for R(s) is 
then usually recasted as a series in the effective parameter a s (s), naively 
obtained by specular reflection, that is by replacing the space-like argument 
Q 2 = —q 2 > with the time-like one s = q 2 > . Starting from 0(oP s ) the 
series displays corrections proportional to 7r— powers, as a drawback of the 
analytic continuation of the hadronic tensor. Then the ordinary perturbative 
expansion for R(s) (e.g. [M|) reads 



R PT (s) = 3}^Q 2 

f 



IX V TX I V TX 



(1.48) 
r 2 = d 2 ; r 3 = d 3 -S 3 , S 3 = ^ (1.49) 

with d 2 and d 3 the same as in (jl.46p . and bj = (47r) J+1 /5 :( - . The correction 
^3 gives to the 0(a^) coefficient a strongly negative contribution for each n/ 
(roughly S3 ~ 14.3 for nf = 4). Higher order 7r 2 -terms were analyzed in [55] 
and show a remarkable growth; for instance, the fourth order correction is 

6i = di _ r4 = f3. ( r2 + ^L^j (1.50) 

16 V 24 W 

with 5 4 ~ 120 for n/ = 4. A similar treatment also holds for other s-channel 
observables [54], such as the normalized rate for r-decay into hadrons. The 



9 In the massless case the cut spreads over the whole positive axis, while taking into 
account quark masses the cut starts at the two-pion threshold 4m^. . 
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effects of analytic continuation then make the perturbative expansions in the 
time-like region different from Euclidean ones. 

Since the 7r 2 -terms play a dominant role in higher order coefficients, expan- 
sion (J1.48P works only asymptotically at large s, that is when the smallness 
of a s (s) scales down these large contributions. Actually, as already noted in 
pioneer works [S6J EZ], the expression of R(s), resulting from Eq. (jl.44p ex- 
hibits no natural expansion parameter (since both real and imaginary parts 
of a s (Q 2 ) enter into the expression of R), and its choice is essentially a mat- 
ter of expediency. 

Alongside less meaningful attempts, one can remind e.g. the analysis per- 
formed in Ref. [56] of the quantity 



o. 



f 

Jo 



-,1/2 



ln 2 (s/A 5 



7T Z 



1.51: 



as expansion parameter for R(s) . Among the main features, Eq. (11.511) turns 
out to be IR finite, thus avoiding the hurdle of Landau singularities on the 
time-like domain, whereas asymptotically, i.e. for s ^> A 2 e 7r , it reads 



1 



a. 



A, ln(s/A 2 



n 



2 ln 2 (s/A 2 



+ 



1.52) 



resembling the UV behavior of the related space-like coupling (11.81) . However, 
Eq. (11.511) cannot entirely sum up the 7r 2 -terms. 

In order to deal with these corrections, a somewhat different approach, 
known as RKP (Radyushkin-Krasnikov-Pivovarov) procedure, has been sug- 
gested [3H EB] (see also [HH]), which is based upon the analytical properties 



of the polarization tensor Il/ l (— q 2 ) and of the related .D-function (11.451) . that 
can be summarized by the dispersion relations, respectively 

R(s) 



m-r 



D{-q< 



ds- 



-q 



ds 



s — q z 
R(s) 



^2^2 



1.53) 



1.54) 



with R(s) given by (11.441) . and q 2 lying in C — {q 2 = s > 0} . The key point 
here is the inverse of Eq. (11.541) , given by the contour integral 



r ^tJZP^ 



1.55) 
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to be computed along a path in the analyticity region for the D-function. 
Eq. (11.55p can be then generalized to an integral transformation, mapping a 
space-like argument function into a time-like one 



R(s) = $ [D(-q' 



1.56) 



This can be straightforwardly applied to the perturbative expansion (11.461) . 
provided that the integration contour is always kept far enough from IR 
space-like singularities. This yields R(s) as an expansion into the images of 
a s (Q 2 ) and of its powers, through the map $ 



R(s) 




J2 d n® 



n>l 



a s (Q 2 



IX 



1.57) 



where the coefficients d n are the same as in Eq. (11.461) . Eq. (11.571) is to be 
compared with the standard perturbative expansion (J1.48P ; here the 7r 2 -terms 
are entirely summed up, with the drawback that, within this framework, there 
is no uniquely defined expansion parameter. However, it is useful to work 
out its behavior to 0(a s ) for R(s), namely 



a {1) (s) 



<S>[a^(Q 2 )} = — \ ---arctan 



In (s/A 2 



7T 



1.58) 



easily obtained by applying the integral transformation (11.551) to the 1-loop 
space-like coupling (11.81) . The related time-like coupling fl 1 . 5 8 [) is free of 
unphysical singularities at low s (as Eq. (I1.5ip ). and for s ^> A 2 e n can be 
expanded into powers of ir/ ln(s/A 2 ) 



a (1) (s) 



TlpQ 



■Ji 



7T 



ln(s/A 2 ) 31n 3 (s/A 2 ) 



1.59) 



Then Eq. (11.591) can be recasted as a power series in the 1-loop space-like 
coupling 



^\s) 



aU(s) 



nX 




1.60) 



emphasizing that the space-like and time-like couplings differ at 3-loop level. 
The comparison between Eq. ( 11 .60}) and Eqs. ( 11. 48ft and ( 11. 49ft makes it clear 
the RKP resummation of the 7r 2 -terms into the time- like coupling ( 11.581) . 
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The main shortcoming of this recipe is that by applying to Eq. ( 11.581) the 
inverse transformation $ _1 , i.e. relation ( 11.541) . the original input (11. 8ft is 
not recovered. Obviously, this is because integral transformations fl 1 . 5 5 [) and 
its inverse are well-behaved as long as the integrand possesses the correct 
analytical properties in the cut complex plane; actually this is not the case 
for the space-like coupling (11.81) and its higher-loop approximations, and one 
is then forced to compute the integral along a path large enough to avoid the 
IR space-like singularities. 



Chapter 2 



The APT theoretical scheme 



2.1 The "Euclidean" ghost-free coupling 

The analytic approach constitutes the next step, after the RG-summation, in 
improving the perturbative results. Specifically, in addition to the property 
of renormalizability, this method retains a general feature of local QFT, the 
property of causality. The analytic approach has first been devised in the 
context of Quantum Electrodynamics [60], and then extended to the QCD 
case about ten years ago [23]. The basic merits of the analytic approach 
to QCD are the absence of unphysical singularities of the invariant charge 
and the enhanced stability of outcoming results with respect to both higher 
loop corrections and the choice of renormalization scheme [61]. Besides, this 
method enables one to process the spacelike and timelike data in a congruent 
way [62]. A fresh review of the analytic approach to QCD and its applications 
can be found in the paper [63]. 

Usually, in the framework of RG-improved perturbation theory a QCD ob- 
servable D(Q 2 ) of a single argument Q 2 = —q 2 > (spacelike momentum 
transfer squared) can be represented as a power series in the strong cou- 
pling a s (Q 2 ) (cf. Eq. flUB}): 

D PT (Q 2 ) = l + Y J d n [a s (Q 2 )] n , (2.1) 

n>l 

where d n are the relevant perturbative coefficients. However, in the IR do- 
main this expansion becomes inapplicable due to the spurious singularities 

28 
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of the running coupling a s (Q 2 ). For example, the 1-loop expression ( 11. 8ft 
possesses both the physical cut along the negative real semiaxis Q 2 < and 
the unphysical pole at Q 2 = A 2 . 

In the framework of the APT, the power series (12. ip for an "Euclidean" 
observable is replaced [65] by the nonpower expansion 

J D APT (g 2 ) = l + ^,rf n A- t (Q 2 ) (2.2) 



n>l 



over the set of functions 



DO 



Pn(or) 



MQ 2 ) = J ff^ d °- ( 2 - 3 ) 



Here, the spectral function p n (<j) is defined as the discontinuity of the relevant 
power of the perturbative coupling a s (Q 2 ) across the physical cut, namely 

p n (a) = -Im[a s (-a-ie)} n . (2.4) 

7T 

The APT representation ( 12. 2ft for a QCD observable D(Q 2 ) is free of spurious 
singularities. Besides, it displays a better stability, in comparison with the 
perturbative parameterization (12.11) . with respect to both, the higher-loop 
corrections and the choice of the renormalization scheme (see Ref. [63] for 
details). The first-order function Ai(Q 2 ) (12. 3p plays the role of the effective 
Euclidean QCD coupling at a given loop level: 

CO 

a E (Q 2 ) = A 1 (Q 2 ) = /-^jdtr. (2.5) 

o 

The argument Q 2 = —q 2 > now runs over the whole space-like axis, that 
is, aE(Q 2 ) is free of any space-like unphysical singularities by construction; 
moreover, due to the asymptotically free nature of the perturbative input the 
spectral integral ( 12. 5ft needs no subtractions. 

In the 1-loop case this equation can be integrated explicitly [23] 



4 1} (Q 2 )- 



00 



A 2 



ln(Q 2 /A 2 ) A 2 -Q 2 



(2.6) 
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The analytically generated non-perturbative contribution in Eq. (12.61) sub- 
tracts the pole in a minimal way, yielding a ghost-free behavior which avoids 
any adjustable parameter. Furthermore, for Q 2 > A 2 this non-perturbative 
term in Eq. (12.61) can be rewritten as 

namely, as a series of power-suppressed corrections to the 1-loop perturbative 
coupling ( 11.81) . Since the pure perturbative contribution dominates in the 
deep UV region, the asymptotic freedom constraint of the invariant coupling 
is preserved. 

On the other hand, Eq. (12.61) exhibits the infrared finite limit a E (0) = l//?o 
(roughly 1.4, if rif = 3). This value turns out to be independent of the 
QCD scale A, and it is also universal with respect to higher-loop corrections, 
warranting a remarkably stable IR behavior of ( 12.61) . 
Obviously the scaling parameter has now to be re-defined as 

A 2 = H 2 ex V [-<p(p a s (fi 2 ))] , (2.8) 

where the function <fi is the formal inverse of ( 12.61) . that is, it satisfies the 
equation 

^TT + 1 l —^ = x - ( 2 - 9 ) 

(p[x) 1 — exp <p(x) 

with x = Po a E and Q 2 /A 2 = exp (j>(x) . In a similar way, the beta-function 
for the 1-loop coupling ( 12.61) reads [23] 



1 exp <j)(x) 

4> 2 {x) [exp0(x)-lf 



A 1) (x) = -^r+ r ~ V ,T'« (2-10) 



with <fi again given by Eq. (12.91) . Despite the implicit form of (12.101) . the 
symmetry property (3 E (x) = f3 E (1 — x) reveals the existence of a IR fixed 
point at x = 1, corresponding to the IR finite limit (see also [66]). 

2.2 Two- loop and higher orders 

Actually, as discussed in Sec. 1.3, the RG equation for the invariant coupling 
has no simple exact solution beyond the 1-loop level, and asymptotic expres- 
sions are commonly used (Eq. (I1.32p ). with a cumbersome IR nonanalytical 
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structure. As a result, the spectral functions (12.41) become rather involved. 
At the 2-loop level, as already pointed out, to get the most accurate result in 
the IR domain, one needs to start with the exact RG solution as given by Eq. 
(11.241) . In the framework of APT, this amounts to consider the discontinuity 
of (11.241) across the time-like cut [63, [66] 



fj(a) = -lm[a^(-a)] , (2.11) 

7T 



and the exact 2-loop analytic coupling then reads 



e* 



4 2) W= / dt-r—rfSV) (2-12) 



1 V + z 

where t = ln(a/A 2 ) and z = Q 2 /A 2 . In spite of its accuracy, Eq. ( 12.12ft 
cannot be easily handled, and a slightly simpler expression can be obtained 
starting from the 2-loop iterative solution (11.271) which turns out to be a 
better approximation to the exact solution than the usual Eq. (I1.30p . as far 
as the IR domain is concerned (see Sec. 1.3). The related spectral density is 
then given by 



^^-^^-Woimfm f = Ha/A2) (2 - 13) 

I(t) = vr + B x arccos - Bl + I , R(t) = t + B x b V(Bi + t)* + n* 

y/(Bi + t) 2 + 7T 2 Bx 

that can be integrated numerically. Moreover, since the singularity structure 
of Eq. (I1.27P is entirely removed by analytization, the iterative 2-loop analytic 
coupling can be straightforwardly obtained by merely adding to (11.271) two 
compensating terms, that cancel respectively the pole and the cut 



4 2 ^) = 4 2) W + Af+Af (2.14) 

A?V 1 ' 

Af\z 



2 A) 1 - z 
1 pp(- B i) d£ B x 



PoJo t-z [Int + B^n^l-B^lnt)] +tt 2 B 2 

with ^ = a J A 2 . From Eqs. ( 12.141) one can readily verify the IR limit a E it (0) = 
l//?o; for general arguments concerning the universality of the freezing value 
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through all orders see for instance [EU [68] and [69] . 

An useful property of the 2-loop coupling (12. 14ft is the possibility to estimate 

its non-perturbative UV tail by expanding the two compensating terms into 

a power series of 1/z = A 2 /Q 2 for large z [70], in a way similar to Eq. (12.71) . 

namely 

1 oo 

Af W + Af(z) = i-£^ (2.15) 

n=l 

2 7 ^(l+£-ln?) +jr 2 

If riy = 6 is taken, the first coefficient, e.g., in Eq. (I2.15P turns out to 
be c\ ~ —0.54. A numerical comparison between Eqs. (12.141) and (12.121) . 
performed in [66] (see Tab.l), reveals that they differ of about 1.8% in the 
IR region. 

Increasing difficulties arise when dealing with even higher-loop levels. In- 
deed, the exact 3- and 4-loop spectral densities ( 12. 4ft have a rather cumber- 
some structure in terms of the Lambert function. Neverthless, an extensive 
numerical study of the analytic coupling ( I2.5P and its "effective powers" ( 12.31) 
up to 4-loop can be found in Ref . [12] . This analysis reveals that the exact 3- 
loop spectral density can be approximated to a high accuracy (~ 1%) by the 
discontinuity of the 3-loop perturbative coupling (as deduced by Eq. (11.321) ) 
across the physical cut [71] : 

pf{o) = j o (t2 + n2)3 [t(^ 2 - t 2 )J(t) - (3t 2 - 7r 2 )R(t)] . (2.16) 

In this equation t = ln(a/A 2 ), 

J(t) = 2t- Si[iGi(t) + G 2 (t)] + B 2 1 G 1 (t)[2G 2 (t) - 1], (2.17) 
R(t) = t 2 - it 2 - fii[iG 2 (t) - ^G-^t)} + 

B 2 lGi(t)-7r 2 Gl(t)-G 2 (t)-l} + B 2 , 

Gi(t) = arctan(-) , G 2 (t) = - ln(t 2 + vr 2 ) 

2 7T V 71 ^/ 2 

with Bj = f3j/Pl + and (3j the expansion coefficients (11.171) of (3 function . 
The 4-loop analogue of Eqs. (12.161) and (I2.17P turns out to be more involved 
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but of small usefulness, being the difference with respect to the 3-loop ap- 
proximation negligible [63]. In Fig. 12.11 the analytic coupling up to 4- loop is 
compared to the perturbative counterpart. 

One can then argue that the 3-loop analityc coupling obtained by integration 
of the spectral density (12.161) , namely 



„(3) 



«£W 



pfV) 
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Q' 



da , 



(2.18) 



is a satisfactory improvement of the 1-loop result (12.61) for all practical aims. 
By normalizing Eq. (12. 18[) at the Z boson mass to the world average value 
a s {M 2 z ) = 0.1176 ± 0.0020 0, one finds A^ ^ 236 MeV. Evolution at 
the crossing of heavy quark thresholds by continuous matching conditions 
then gives A„ =3 ~ 417 MeV. Eq. (12.181) with the spectral density given 
by Eqs. (12.161) and (I2.17P is actually the expression used in the following 
analysis (see Sec. 3) for the comparison with the experimental values of the 
QCD coupling extracted from the Bethe-Salpeter approach. 

Moreover, it should be noted that an integral representation, similar to Eq. 
(12.141) (and also Eq. ( 12.151) ). of the non-perturbative terms accounting for 
the divergencies of the 3-loop asymptotic solution (11.321) (that is the leading 
singularity in z — 1 ( 11.4011 . and the IR log-of-log generated cut) also exists [70], 
E2j. It turns out [HSj that the asymptotic behavior of a E (Q 2 ) defined by 
Eq. (12.181) . i.e. in the large Q 2 limit, is given by 
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Here 7e is the Euler constant, B 2 = Pofo/P 2 , a s (z) is the perturbative 
counterpart as given by (11.321) . and c\ ~ —0.52 if rif = 6 (see Ref. [69] for 
details). 

high-accuracy 3- and 4-loop analytic cou- 

(2.20) 



Note finally that in Refs. [42 
pling of the form 
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Figure 2.1: Analytic coupling taken from [69j up to 4-loop compared to perturbative 
counterparts (|1.32[) . 



has been also worked out, by exploiting the multi-loop approximation fl 1 . 3 3 [) 
as an input in the dispersion relation (I2.5p . Here a E is given by Eq. (I2.12p 
and the coefficients are the same as in (J1.33J) . Numerical values at low scales 
can be found in [12]. 

For practical applications, one can even resort to simple expressions of the 
form of Eq. (12. 6p . with a special model argument. According to Ref. [72], the 
3-loop Euclidean coupling can be approximated with reasonable accuracy by 
the "1-loop-like" model 
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where 

L2(Q 2 )=^{jP)+B 1 In Jin 2 (g) + 2^, Bl = | . (2.22) 

The relative difference between the exact 3-loop analytic coupling (numeri- 
cally evaluated via the Lambert function) and its approximation ( 12.21D is less 
than 3% for Q > 500 MeV (see Table 1 in Ref. [72] for the details). However, 
this error could reach 5-10% below 500 MeV and it is then unuseful in the 
deep IR domain. In this region the more precise expression (12.181) should 
be used. Moreover, its accuracy in the asymptotic region breaks down when 
taking into account flavor thresholds. Therefore, it has been suggested [73] 
to use Eqs. (I2.2ip and (I2.22p provided that the scaling constant and the co- 
efficient B\ are replaced by adjustable parameters (see tab. Ill in Ref. [73]). 
This yields an accuracy within 1% in the whole space-like region down to 
Q > 1 GeV. 



2.3 The "Minkowskian" domain in APT 

As far as the definition of a reasonable expansion parameter in the time-like 
domain (q 2 = s > 0) is concerned, it should be noted that such a definition 
naturally arises in a self-consistent way within the framework of APT [7^175]. 
Further, it can be regarded as the final step in the RKP-resummation for ir 2 - 
terms outlined in Sec 1.6. 

As already noted, Eq. (I1.54p and its formal inverse (11.551) can be generalized 
to the proper tool for relating s- and t-channel observables. In the APT 
scheme, a standard power expansion for a space-like observables ( 12. ip is re- 
placed by a nonpower one (12.21) . A time-like observables, in turn, within this 
framework, can be also re-expressed as the nonpower expansion 

Rapt{s) = 1 + J2 dn%n{s) (2.23) 

over the set of functions 
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Figure 2.2: (a) Space-like and time-like global analytic couplings in a few GeV domain 
with rif — 3 and A = 350 MeV; (b) "Distorted mirror symmetry" for global expansion 
functions, corresponding to exact 2-loop solutions. 



where the n-th spectral density p n is again given by Eq. (12.41) . This recipe 
is manifestly quite analogous to the RKP non-power expansion (J1.57J) . The 
key point here [76] is that, due to the forced analyticity of the coupling 
and its analytized powers (12.31) . the two sets (12.31) and (I2.24p are put into a 
one-to-one relation by the linear integral transformations (11.54)) and (11.551) . 
namely 

A(Q 2 ) = D[2l n (s)], and 2l n (s) = R[A n (Q 2 )] . (2.25) 

This eventually yields a closed theoretical scheme for representing observables 
of any real argument, both space-like and time-like (for technicalities see |77j). 
The main features of the two functional sets are illustrated in Fig. 12.21 taken 
from [77] . For n > 2 both sets A n and 2l n start with an IR zero and oscillate 
in the IR domain around n — 1 zeros; furthermore, they all display the UV 
asymptotic behavior l/ln 71 ;?, resembling the corresponding powers a™(z) . 
The recursion relations 
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(2.26) 

where I is the loop level, allow one to relate different analytized powers within 
each set, albeit explicit expressions are available in a simple form only at 1- 
loop (see e.g. [72]). Nonethless, the two sets ( 12. 31) and ( 12.241) have been 
numerically evaluated for n = 1,2,3 in [67] up to 3-loop, via the Lambert 
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function (see also [12] and [12]). Furthermore, a generalization of APT to 
fractional powers of a s was also implemented in Ref. [75] . 

As in the space-like domain, the first function of Eq. (j2.24p plays the role of 
the effective coupling of time-like argument, namely 

«a/00= / — px{<r). (2.27) 

Note also that «m(s) can be equally defined by the differential equation [11] 

s-t^m(s) = —pi(s) with <S(oo) = , (2.28) 

as 

as can be verified by differentiating the second of Eq. (I2.25J) with n = 1 , i.e. \j 

i C s+ie dO 2 
^m{s) = —J -^aEiQ 2 ). (2.29) 

Moreover, Eq. (I2.28J) emphasizes the straightforward relation between the 
"time-like /3-function" and the spectral density, thus reviving an old hypoth- 
esis due to Schwinger [79] . 

At 1-loop from Eq. (12.271) one finds again Eq. (11.581) ; however this now leads 
by inversion of Eq. (I2.29|) to the starting space-like coupling (12.61) . being the 
analytic properties preserved within this framework. 

The main feature [71] of Eqs. (12.271) and (12.51) is the common IR freezing 
value <S(0) = a^(0) = l//? , independent of the loop level and of any ad- 
justable parameter (see Fig. 4(a)). Moreover, they exhibit a similar leading 
UV behavior, constrained by asymptotic freedom. Nevertheless this approx- 
imate "mirror symmetry" is broken in the intermediate energy region □, the 
discrepancy being about 9% at 1-loop, and slightly less at a higher loop level 
(see Ref. [H] for numerical comparisons). 

A last remark concerns thresholds effects, that can be included within the 
APT algorithm [75] by modifying the n-th spectral density (12 .4p discontinu- 
ously at the heavy quark thresholds rrif, namely 

p n {a) = p n (a, 3) + ^ °( a ~ m /) \Pn{v, rif) - p n (a, n f - 1)] , (2.30) 

n f >A 



^^Note that this definition for the time- like coupling had been suggested also in [591 [74] , 
2 In Ref. [80] an argument, on the basis of causality principle, against a possible exact 
symmetry ruling the "t-s dual" couplings, has been developed. 
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descending from the trivial matching condition (Sec. 1.4). The global func- 
tions resulting from densities ( 12.301) . A n and 2l n , can be obtained from the 
local ones with rif fixed, by adding specific shift constants c n (rif) , not negli- 
gible in the rif = 3,4 region (technical details are given in [76]; for instance, 
in both the t- and s-channels it turns out Ci(3) ~ 0.02 ). 
The main tests of APT being obviously at low and intermediate scales, a 
number of applications to specific processes, both in the space-like and time- 
like (low and high energy) domains, have been quite recently performed (see 
for instance, Refs. [ED H2J [751 ED (TTJ [73]). As a result, the main advantages 
of the APT approach are better convergence properties of the ghost-free non- 
power expansion with respect to the usual power one, and a reduced scheme 
and loop-dependence. 

Finally, transition from Euclidean (space-like momentum) to the distance 
picture has been also developed in [82], involving a suitable modified sine- 
Fourier transformation, consistently with the APT logic. 



2.4 Modifications of APT and the massive 
case 

Different strategies to incorporate analyticity into the RG formalism, or even 
to implement the above device, exist as well. A number of approaches are 
based on non-perturbative constraints either on the invariant coupling (see 
e.g. Refs. [Ml M E2]) or on the RG (3 function (see, e.g., Refs. [851 ED EZ] 
and also [881 [39]) • 

Among the many attempts to cure the Landau ghost problem it is remark- 
able the existence of models that suggest an IR enhancement of the QCD 
coupling, invoking analycity as well. The most attractive feature is supposed 
to be a straightforward relation with quark confining potential (see Sec. 1.1). 
As an example, it can be reminded the "synthetic coupling" model recently 
developed in [85], which amounts to modify the analytically improved cou- 
pling (12.51) by additional non-perturbative pole-type terms; at 1-loop it reads 



a syn (<2 ) = -5- 



A 2 cA 2 (1 - c)A 2 



+ 



A, Lln(Q 2 /A 2 ) A 2 - Q 2 Q 2 Q 2 + m 2 J 



(2.31) 



where m g = A/yc — 1 . The IR slavery, due to the pole term at Q 2 = , 
is controlled by one dimensionless parameter c G (l,oo), which relates the 
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scaling constant A to the string tension a of potential models (cf. Eq. (11.121) ). 
The pole term at Q 2 = —rn 2 < corresponds to a non vanishing " dynamical" 
gluon mass, while leaving the analytical structure of eq. (12.61) along the space- 
like axis unchanged. The value of the m g parameter has been estimated [85] 
as 400 — 600 MeV. Eq. (12.311) can be derived, analogously to (12.61) . from a 
dispersion relation with a spectral density of the form (12. 4p plus two 5-terms 
properly accounting for the poles. Along with the IR enhancement as 1/Q 2 , 
reproducing the linear confining part of the potential fll. 12|) . construction 
of (12.311) is mainly motivated by the UV asymptotic behavior of its non- 
perturbative contribution, of the form l/(<5 2 ) 3 , decreasing faster than (12. 7p 
as Q 2 — > oo (see also [8"9"]). 

Actually, of some relevance to the following analysis (see Sec. 3), in par- 
ticular, as far as the deep IR behavior of the QCD coupling is concerned 
(i.e. below 200 MeV), is the "massive" modification of the APT approach 
recently developed in [90]. In order to handle the QCD observables which 
do not satisfy the integral representation of the form of Eq. (12. 3p . the APT 
formalism has to be modified appropriately, by taking into account the effect 
of a non- vanishing mass threshold m in the dispersion relations. In this case 
the set of APT expansion functions (I2.3P should be replaced by the set of 
"massive" ones with an adjustable parameter m 

oo 

2 2 \ Q 2 f f ^(y — Am 2 da 



^ Q ' m ^c?f^J^^Tc?V (2 ' 32) 

4m 2 



The APT expansion ( 12.21) should be then replaced by 

D MAPT (Q 2 , m 2 ) = ^' 2 + J2 d n A n(Q 2 > ™ 2 ) • (2-33) 



Q 2 + Am 2 



n>l 



Obviously, in the massless limit Eqs. (I2.33J) and (I2.32p coincide with the 
expressions (12.21) and (12. 3p . respectively. 

Similarly to the case of Eq. (12. 5p in the "massless" APT, the first-order 
function Ai(Q 2 ,m 2 ) (12.321) plays the role of an effective "massive" running 
coupling at the relevant loop level, namely 



-x: 



4m 2 
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It is worthwhile to note that irrespective of the loop level this coupling pos- 
sesses the universal IR limiting value a(Q 2 , m 2 ) — ► at Q 2 = , in qualitative 
agreement with some results from lattice simulations [T71 [16] (see Ref. [90] 
for details). Moreover, as discussed in Sec. 3, the IR behavior of the QCD 
coupling, in particular below 200 MeV, extracted from the meson spectrum 
in BS formalism, can be reasonably described within this "massive" modi- 
fication of APT, that is, by Eq. (j2.34p if a proper value of the adjustable 
parameter m is chosen. 



Chapter 3 

Bound states approach 



3.1 QCD coupling from the meson spectrum 



Below 1 GeV, as discussed in Sec. 1, the running coupling a s is affected at 
any loop level by unphysical singularities that make the RG-improved pQCD 
useless. 

In order to shed some light on the QCD coupling below this scale, one needs 
experimental information to be extracted with the aid of a suitable theoretical 
framework. The latter can be provided by the Bethe-Salpeter (BS) formal- 
ism developed in [23] for the study of the meson spectrum, since the scale 
involved (the momentum transfer in the qq interaction) is typically below 
1 GeV. Many relativistic formalisms for the analysis of the meson properties 
have been developed in the context of QCD, that take confinement into ac- 
count and evaluate the meson (and baryon) spectrum in the light and heavy 
quark sectors. Among the most recent works, one should remind for instance 
Refs. [251 ES E3 M, EU [921 EH] and references therein. 

In the present work the theoretical results on the meson spectrum within the 
framework of a Bethe-Salpeter (BS) formalism adjusted for QCD, have been 
exploited with the aim of extracting "experimental" values of the running 
coupling al xp (Q 2 ) below 1 GeV by comparison with the data. 
The second order BS formalism, developed in [23] and applied to the calcula- 
tion of a rather complete quarkonium spectrum in Refs. [Ml [95] , is essentially 
derived from the QCD Lagrangian taking advantage of a Feynman-Schwinger 
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representation for the solution of the iterated Dirac equation in an external 
field. Confinement is encoded through an ansatz on the Wilson loop corre- 
lator (see Ref. [21] for details); indeed the expression ilnW^ is written as the 
sum of a one-gluon exchange (OGE) and an area terms 

ikiW=(i]nW) OGE + aS. (3.1) 

By means of a three dimensional reduction, the original BS equation takes 
the form of the eigenvalue equation for a squared bound state mass 

M 2 = Ml + Uoge + Uconi , (3.2) 

where M is the kinematic term M = w\ + u>2 = \/ m i + k 2 + y/m 2 , + k 2 , 
k being the cm. momentum of the quark, mi and rri2 the quark and the 
antiquark constituent masses, and U = Uoge + ^Conf the resulting potential. 
As a consequence of ansatz (13. ip . the perturbative part of the potential £/oge 
turns out to be proportional to a s (Q 2 ) , that should be regularized with a 
proper prescription. 

Calculations have been performed in Refs. [9lJ EH] by using both a frozen and 
the 1-loop analytic coupling a E (Q 2 ) as given by Eq. (12. 6ft with an effective 
scaling parameter A^ '!l 3 ~ 200 MeV. 

The results of the two sets of calculations are relatively similar for the heavy- 
heavy quark states. However, for the IS states involving light and strange 
quarks, quite different results have been obtained in the two cases. In the 
case of a frozen coupling the ir and K masses turn out to be too high, ir- 
respective of how small the light quark mass is taken (see Fig. 13.11) ; e.g., if 
the light and the strange quark masses are fitted to the p and the masses, 
one finds m, ~ 500 MeV and m^ ~ 700 MeV, respectively. On the contrary, 
if appropriate values for the quark masses are chosen, the ir, p, K, K*, <fi 
masses can be rather well reproduced when the analytic coupling a E (Q 2 ) 
is used. This occurrence strongly supports the use of the APT prescription, 
outlined in Sec. 2, within BS framework. 

Indeed, the combined BS-APT theoretical scheme provides a substantial 
agreement of the calculated spin averaged c.o.g. masses with the data through- 
out the whole calculable spectrum, and in particular the correct splittings 
l 3 < S'i-l 1 5o is well reproduced at least in the light-light, light-strange and 
heavy-heavy sectors. 
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Figure 3.1: Quarkonium spectrum, three different calculations. Diamonds refer to the 
truncation prescription for the running coupling, squares and circles refer to the calculation 
with the 1-loop APT coupling (|2.6p and two different expressions for running constituent 
masses of light quarks, a solution of the Dyson-Schwingcr equation and a phenomeno- 
logical function of the cm. quark momentum respectively. Horizontal lines represent 
experimental data. 
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In this work a thorough analysis has been undertaken from the reversed point 
of view, that is by comparing theoretical results of the meson spectrum within 
the BS approach, obtained for a definite choice of the parameters, with the 
results of a similar calculation performed by means of a fixed value of a s , 
for every quark- ant iquark state. In what follows af 1 denotes the value that 
reproduces the same theoretical result as by using a E (Q 2 ) , whereas a^ xp 
the value that correctly reproduces the experimental mass. The value a* h is 
clearly an intermediate step object, used to identify an effective Q pertaining 
to each meson state, to be understood as the argument of the related "ex- 
perimental" coupling afv{Q 2 ). 

Since only the leading perturbative contribution in the BS kernel has been 
included, a rough estimate of NLO effects on the extracted af^ p values leads 
to a relative theoretical error which spans from 20% to much less than 4% 
throughout the spectrum, according to the quark masses involved. Further- 
more, since coupling among different quark- ant iquark channels has not been 
taken into account within BS formalism, the theoretical masses are expected 
to reproduce the experimental ones roughly within the half width T/2 of the 
state. If relevant, the experimental error, related to the uncertainty of the 
experimental masses, is added to the theoretical one. 

The results are twofold. On the one hand, the 3-loop APT coupling, nor- 
malized at the world average a s [M\ ) = 0.4476(20) and evolved across heavy 
quark thresholds, reasonably fits a^ xp (<3 2 ) from 1 GeV down to 200 MeV, 
quantitatively confirming the relevance of the APT approach to IR phe- 
nomena down to 200 MeV. On the other hand, below this scale, with the 
limitation due to the large errors, the experimental points give a wee hint 
about the vanishing of a s (Q 2 ) as Q — ► 0. This could correlate with some 
results from lattice simulations [47], and can be theoretically discussed in 
the framework of a recent "massive" modification [90] of APT (see Sec. 2.4) 
which takes into account effects of a finite threshold in the dispersion rela- 
tion 0. 



1 I"t should be noted, however, that the existence of a non vanishing finite limit lower 
than the universal Shirkov-Solovtsov freezing value, as suggested in [27], is still consistent 
with the above results. 
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3.2 BS-model for quarkonium states 



As mentioned, in [94] 195] the meson spectrum is obtained by solving the 
eigenvalue equation for the squared mass operator (13.21) . where the pertur- 
bative and the confinement parts of the potential are respectively 
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Here a^ denote the usual Dirac matrices 7? 7*, erf the 4x4 Pauli matrices 



k+k' 



, Q = k - k' , g j0 = - 



13 ' j 

— — 2"^-, mi and 777.2 are 

constituent masses. 

Eqs. f l3.3H3.5p have been derived from the ansatz ( 13.11) and a 3-dimensional 
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reduction of a Bethe-Salpeter like equation (see Refs. [211 Ell [95] and App. A 
for the details). Actually in the calculation of Refs. [9H EH] only the center 
of gravity (c.o.g.) masses of the fine multiplets were considered as a rule, 
and the spin dependent terms in (13.3H3.5p (spin-orbit and tensorial terms) 
were neglected with the exception of the hyperfine separation term in (13.31) , 
proportional to |o"i-cr2 . Within this limitation, a generally good reproduction 
of the spectrum was obtained for appropriate values of the parameters, as 
shown in Fig. 13.11 Here the results of three sets of calculations are displayed. 
Diamonds refer to the usual perturbative 1-loop coupling (to be replaced 
in Eq. (13.31) ). frozen at a maximum value H , which has been taken as an 
additional adjustable parameter. Squares and circles both refer to the 1- 
loop APT coupling (12.61) with A ~ 200 MeV. For light quarks a running 
constituent mass was used tooo 

It should be stressed that only with the choice (12.61) the V-Sq state has been 
correctly reproduced when light and strange quarks were involved, as in the 
case of 7r and the K mesons. 

In the present work [31 H] a similar calculation with the input (12.61) and a 
slight different choice of the parameters is made (preliminary results were 
given in [96]). First, the string tension has been fixed a priori to the value 
a = 0.18 GeV 2 (consistent with hadron phenomenology and lattice simula- 
tions) and the scale constant to A„ '!_ 3 = 193 MeV. The whole set of remain- 
ing parameters, all the quark masses, are then determined by fitting the it , 
, J/ip and T masses. It turns out m u = m^ = 196 MeV, m s = 352 MeV, 
m c = 1.516 GeV and rrif, = 4.854 GeV. The results for the meson spectrum 
are given in the fourth column of the tables in App. B. 
The effective value for the QCD scale A^ ' e = ^ has been dictated by the com- 
parison with the 3-loop analytic coupling normalized at the Z boson mass 
(see Eqs. ( I2.16|I2.17I) and (I2.18P ) according to the world average. As displayed 
in Fig. 13.21 the relative difference between the two curves is no more than 1% 
in the region of momentum transfer 0.5 < Q < 1.2 GeV, to which the bulk 
of the states used as an input in the calculation belongs. 
Furthermore, as already noted, the set of equations (l3.3H3.5p refer to a single 
definite quark-antiquark channels. So, having correct relativistic kinematics, 



2 Circles refer to a phenomenological running mass as function of the cm. quark momen- 
tum m 2 = m 2 d = 0.17|k| — 0.025|k| 2 + 0.15|k| 4 . Squares refer to a running constituent mass 
resulting from a solution of the Dyson-Schwinger equation with an analytic RG running 
current mass (see App. A for some details). 
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Figure 3.2: Relative difference between the one-loop analytic running coupling a^ (Q 2 ) 
with A^'2 = f93McV and three-loop a§ ] {Q 2 ) with A^ 3) =3 = 4f7McV in the range 
< Q < 1.2 GcV. 



they do not include coupling with other channels like any potential model 
(see App. A). Therefore, one can not expect to have any insight into the 
splitting of over-threshold complicated multiplets which involve mixture of 
different states. Even the position of the c.o.g. mass is expected to be repro- 
duced only within one-half of the width of the state. This has been taken 
into account in the estimate of the theoretical error (see Sec. 4). 

The resolution method of the eigenvalue equation for the operator (13.21 13.31 - 
13.51) used in [9U [95] and in the present work can be summarized as follows. 

a) In the static limit the problem can be reduced to the corresponding one 
for the center of mass Hamiltonian (see App. A) 



Hqm =w 1 +w 2 



4a s 

-— + ar. 

3 r 



(3.6) 



b) The eigenvalue equation for (13.61) is solved for a convenient fixed a s by 
the Rayleigh-Ritz method, using the three dimensional harmonic oscillator 
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basis and diagonalizing a 30 x 30 matrix. 

c) The square of the meson mass is evaluated as (0 a |M 2 |0 a ), a being the 
eigenfunction obtained in step b) (with a the whole set of quantum numbers) 
and the operator M 2 given by Eq. (13.21) . 

d) Prescription c) is equivalent to treat M 2 — Hq M as a first order perturba- 
tion. Consistently the hyperfine separation should be given by 

<*m nl ) 2 - ( l m nl ) 2 = — ^ dkk 2 J dk> k%* nl (k)<p nl (k>) x 




Wi + U>2 w[ + w' 2 f 1 



P£ [ ^« S (Q 2 )^(0, (3-7) 

Wo J-l 



W\W2 V W\W 



where <p n i is the radial part of the complete eigenfunction <p a . 
For the quark masses and string tension a in (I3.6P the same values listed 
above have been used, and as far as a s is concerned, that is supposed to 
be a constant in fl3.6|) . the value a s = 0.35 has been used, which is the one 
typically used in non-relativistic calculations, and is also the freezing value 
adopted in 



3.3 Extracting c^ xp (<2 2 ) from the data 

One focus now on the reversed problem, i.e., the determination of the a e ^ v {Q 2 ) 
values at the characteristic scales of a selected number of ground and excited 
states. 

In order to estimate a^ xp (Q 2 ) a ^ l° w energy one needs first to assign an effec- 
tive Q-value to each state. To this end one first rewrites the squared mass, as 
given by point c) in Sec. 3.2, more explicitly as the sum of the unperturbed 
part, the perturbative and the confinement ones respectively 

m 2 a = (0 a |M o 2 |0a) + (0a|£/oGE|0a> + (<Pa\U C ont\<Pa) • (3.8) 

Here Uoqe is given by the second line of (J3.3P and C/conf by Eq. (13.41) and 
the first two lines of (j3.5p . From the OGE contribution one then extracts 
for each state the fixed coupling value a* h that leads to the same theoretical 
mass as by using q^ (Q 2 ) given by Eq. (12.61) . This can be done by means of 
the relation 

(<Pa\U OG E\<Pa) = (^a^ (Q 2 )0(q; Q)|0 a ) = a* h (0 a |O(q; Q)|0 a ), (3.9) 
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where 0(q; Q) can be drawn again by the second line of Eq. ( 13.31) . The effec- 
tive momentum transfer Q a associated to each bound state is then identified 
by the equation 

4\Ql) = ^. (3.10) 

The next step is to search for the correct (fixed) value of the coupling that 
exactly reproduces the experimental mass of each state. This is defined by 
the relation 

(0 a |M o 2 |0a) + «r P (Q;D(0 a |0( q; Q)|0a) + (<Pa\U C onf\<Pa) = TT&p , (3.11) 

so that, by combining Eqs. ( 13.81) . ( 13.91) and ( 13.111) one finally obtains 

2 2 

777 — 777 

(0a|C(q;Q)|0a) 

This procedure has been applied to a number of light-light, light-heavy and 
heavy-heavy ground as well as excited states. 

Before discussing the uncertainties related to the extracted «° xp (Q 2 ) , few 
comments are in order. 

First note that in the evaluation of Q a in (13.91) the hyperfine splitting has 
been neglected, whereas it has been taken into account in (13.121) . bringing 
possibly to different values of o^ xp for the singlet and the triplet states (when 
there are reliable data for both). 

Furthermore, the sensitivity of the effective Q's derived as above from the 
specific coupling (12.61) . has been checked by analyzing their deviation for 
a 25% shift of A„ '^ 3 around the value 193 MeV, and one finds that the 
average change in the momentum scale amounts to 3% . This makes the 
resulting a| xp (Q 2 ) reliable, at least qualitatively, even in the deep IR region 
(Q < 0.2 GeV), where the discrepancy with respect to massless a E (Q 2 ) is 
sizable. 

Obviously the theoretical meson masses are sensitive to a variation of the 
quark mass parameters (particularly in the case of the tt meson), whereas 
a^ xp and the relative Q a turn out to be much more stable. For instance, in 
the light-light sector, an increase in the light quark mass of 5% amounts to 
a change of about 2% in the value of a^ xp and 0.2% in the relative Q a . 
Finally, a subtle point concerns the choice of the "unperturbed" a s involved 
in the static Hamiltonian (13.61) . Actually, the value adopted is very near to 
the a* h pertaining to the bb(lS) state, but definitively smaller than the typical 
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a* h 's. The point is that the hyperfine splitting is much more sensible than 
the c.o.g. mass to the behaviour of the unperturbed wave function at small 
distance (large momentum), which is specifically controlled by the value of 
the unperturbed a s . As a result, the effective fixed value a s in Eq. (J3.6P that 
reproduces the same splitting as by using the coupling a^ (Q 2 ) turns out to 
be significantly smaller than «* h calculated from the c.o.g. mass. Essentially, 
it was chosen a phenomenological value for the unperturbed a s in order 
to have a good reproduction of the hyperfine splitting so as to reasonably 
reconstruct the c.o.g. of the doublet when one component is missing. It 
was then used the position of the c.o.g. (which is rather stable w.r.t. the 
unperturbed a s ) to extract the final result af cp (Q 2 l ) . 



3.4 Theoretical uncertainties 

First of all, it should be clearly addressed the unavoidable model dependence 
of the above results. Indeed, the formalism exploited in their derivation is 
based on the ansatz (j3.1| . consisting of the sum of two contributions that 
one knows to be asymptotically correct for small and large quark-antiquark 
distances. More sophisticated ansatz also exist (see e.g. [HI], [S2] or [93"]). 
but they turn out to be difficult to implement within BS formalism. 
In the context of the model in hand the sources of error are the istanta- 
neous approximation implied by the three-dimensional reduction of the BS 
equation, the approximations introduced into the resolution of the eigen- 
value equation, the inclusion of only the leading perturbative contribution in 
the BS kernel J, and finally having neglected the coupling between different 
quark antiquark channels. 

Thus, aside from minor effects due to approximations in the resolution of 
the eigenvalue equation and from the three-dimensional reduction, the main 
sources of theoretical error in the whole procedure are expected to arise from 
neglecting the NLO contribution to the BS kernel as well as the coupling 
with other channels, and these have been explicitly estimates. 

The NLO contribution to the perturbative part of the BS-kernel comes from 
four diagrams with two-gluon exchange; two triangular graphs containing a 
four-line vertex of the type g 2 <j)*<j)A tl A tl and two three-line vertices g^d^A^ 
(the spin independent part of the second order BS formalism is quite similar 
to scalar QED), one fish diagram with two four-line vertices, and a crossing 
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(a) 




(b) 




(c) n 




(d) 






Figure 3.3: NLO contributions to the second order BS kernel I. 



box with four three line vertices (see Fig. I3.3J) . If the renormalization scale 
is identified with the momentum transfer Q the fish graphs contribution is 
completely reabsorbed in the renormalization. On the other hand, a some- 
what crude estimate of the contribution of each of the two triangular graphs 
yields 

9m2 (3.13) 



-'triang nu 



4 



4Q 2 + 2m 2 



and for the crossing box graph, similarly 
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These expressions have to be compared with the leading one-gluon term used 
in the present calculation 

/oge ~ I6n-a s — . (3.15) 

Putting all things together, the overall error due to the omission of such NLO 
contributions to the BS kernel is then 



(3.16) 

1 \j \ J OGE / V -<OGE / 

and this produces 

(3.17) 

By using Eqs. (I3.9H3.121) . after some algebra it is easy to recognize that the 
NLO effects on a° xp turn out to be of the same order, that is 

A NLO a,~a* h — , (3.18) 

which is what is assumed in the foregoing. The NLO errors do not exceed 
5% for heavy quark states whereas they are enhanced up to 20% when light 
and strange quarks are involved. 

Finally, since the strength of the neglected coupling with other channels (OC) 
is obviously measured by the width T a of the state, one roughly estimates an 
error of the order of Am a ~ T a /2 in the evaluation of the theoretical meson 
mass m a . On this basis, for each determination of a| xp (Q 2 .) the related 
theoretical error is given by 

(0 o |C>(q;Q)|0 o ) 

Usually the error Am cxp on the experimental mass m exp is much smaller 
than T a /2 . When, however, this is not the case one has to consider also the 
experimental error A exp a s , obtained from (13.191) by replacing m a Y a with 

All other sources of uncertainties aformentioned, that is, the errors implied 
by the three dimensional reduction, and by the approximations introduced 
in the resolution of the eigenvalue equation for M 2 , and the one due to the 
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model dependence, though difficult to be explicitly estimated, can be globally 
taken into account by means of a comparison between the calculated and the 
experimental spectrum. To this aim, one can restrict the considerations to a 
sample of better established data which exclude high orbital excitations (D 
and F states), and introduce an additional conventional error Am on all the 
theoretical masses, defined by 



2 . _L 

1\ S P 



N SP 

^(m a - m oxp )7[(A tot m a ) 2 + (Am) 2 ] ~ 1, (3.20) 

o=l 



where A tot m a is the total error resulting from all the sources explicitly eval- 
uated, i.e., (Atotma) 2 = m aA//7 + (r a /2) 2 + (Am cxp ) 2 , and the sum is re- 
stricted to the selected sample of data (namely, S and P states only). It turns 
out Am ~ 20MeV, and then it can be set Aa s = (2m a Am)/(</> a |0(q; Q)|0 O ) • 
This uncertainty, that turns out to be about 5% on average, has not been 
included in the tables, but is used in the analysis of the following section. 



3.5 BS-model results: concert of low and high 
energy data via APT 

All results are displayed in details in tables I- VII of App. B, and pictorially in 
Fig. 13.41 The first three columns specify the state and its experimental mass 
as given by [2]. The fourth column gives theoretical results for the meson 
masses, and the last three give the effective Q 's, the relative 3-loop APT 
coupling a E (Q 2 ) and the experimental coupling with errors (both theoretical 
and experimental). 

In Fig. 13.41 values of a^ xp at the same Q from triplet and singlet states have 
been combined through a weighted average according to their errors. 

As one may infer from Fig. 13.41 the experimental points exhibit a remarkable 
evolution from 500 MeV down to 200 MeV, where only the safer S and P 
states are involved, in good agreement with the 3-loop analytic coupling 
a E (Q 2 ) properly normalized, i.e., with A;J =3 = 417MeV, and discussed in 
Sec. 2 (cf. Eqs. ( l2~T8l) and (l2~T6l) . 1J2T7J1 ). 
Specifically, defining A tot a s as the total error explicitly evaluated by means 
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of Eqs. (ETBD and (13191) . one finds 

&=N- E« P (^) - 4 3) (QD)V[(A tot « s ) 2 + (Aa s ) 2 ] ~ 0.8 . (3.21) 

a=l 

The agreement quantitatively supports the relevance of the APT approach 
to IR phenomena down to a few hundred MeV. 

On the other hand, at energies below 200 MeV there seems to exist a general 
tendency of af^ p (Q 2 ) to deviate from the APT curve and to approach zero, 
or at least a finite limit, which is less than the universal APT freezing value. 
Note, however, that the analysis of such an extreme IR behaviour is based 
on high orbital excitations (D and F states), that lie well above the strong 
decay thresholds and possess large widths. As a consequence, the theoretical 
reliability of the method is lower at these scales, as apparent from the large 
estimated errors. Moreover, also the discrepancy between a^ (Q 2 ) (used in 
the calculation) and a^ (Q 2 ) (used as a reference term) rises above 10% 
at these scales. In fact, only two states 7r 2 (1670) (interpreted as ss(l 1 D 2 )) 
and / 2 (2150) (ss (1 3 F 2 )), corresponding to Q ~ 120 MeV, generate af p (Q 2 ) 
(marginally) out of the error bands, and the state / 2 (2150), which has been 
observed only once, has never been confirmed. 

Nevertheless, as already noted, such a deep IR behaviour could correlate with 
some lattice results [17], and could be discussed within theoretical models, 
in particular within the recently developed "massive" modification of APT, 
briefly discussed in Sec. 2.4. Specifically, as displayed in Fig. I3.4[ the 1- 
loop coupling a(Q 2 ,m 2 ) (I2.34p with an effective mass m eff ~ (38 ± 10) MeV 
reasonably fits all experimental points down to the deep infrared region. 

Notice that in the selection of states irregular and incomplete multiplets have 
been excluded as a rule. Of this type, e.g., in the light quark sector, are the 3S 
states (m 3 3 Sl —m 3 i So is anomalously large and about twice as m 2 3 5l — m 2 i 5o ), 
1 3 P (m!3p being larger than mxspj, 1 3 D, F, G, H (incomplete). If however 
included in the analysis, all these states would bring the results in agreement 
with the general tendency outlined. 

As far as the renormalization scheme dependence is concerned, note that the 
coupling definition is implicitly embodied in ansatz (13. ip . Specifically, here 
one assumes that after the area term subtraction, i In W is dominated by 
the OGE term, with the fixed coupling a s replaced by the running coupling 
a s (Q 2 ) . This amounts to including all the dressing effects into a s (Q 2 ). It 
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Figure 3.4: Extracted values of af p (Q 2 ) against the 3-loop APT coupling (|Q|) with 



A 



(3) 

I " _ 3 = (417±42)McV (solid), and its perturbative counterpart (dot-dashed). The 



L (l,eff) 



"massive" 1-loop APT coupling (n = 1 in ([232]) ) refers to A^ 'Zj = 204 MeV and m cff = 
(38 ± 10) MeV (dashed). Circles, stars and squares refer respectively to qq , ss and qs 
with q = u,d , diamonds and crosses to cc and 66 ; asterisks stay for qc and qb , whereas 
plus signs for sc and sb . 



is worth noting that the coupling defined in such way is free of unphysical 

singularities by construction. 

At the same time, the 1-loop effective APT coupling a E (Q 2 ), involved in 

the calculation of the meson spectrum, is remarkably stable with respect to 

both the higher loop corrections and the choice of renormalization scheme 

(see Sec. 2. and a detailed discussion of this issue in Ref. [63])- Thus one 

might expect that the same situation should also occur for a cxp (Q 2 ), with the 

possible exception for the deep infrared region, where other nonperturbative 

effects could be relevant. 

Quark self-energy effects have been taken into account by a recursive resolu- 
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tion of the Dyson-Schwinger equation. In the second order BS formalism this 
simply amounts to replacing the current quark masses with the constituent 
masses [95]. 



1.0r 




Q [GeV] 



Figure 3.5: Summary of low (o) and high energy (o) data against the 3-loop analytic 
coupling (|2.5j) (solid curve) and its perturbative counterpart (dot-dashed curve) both nor- 
malized at the Z boson mass. Also shown is the "massive" 1-loop analytic coupling (|2.34l) 
(dashed curve) same as in Fig. 



Finally, Fig. 13.51 displays a synthesis of the results for a s (Q 2 ) extracted from 
the bound states calculation in the BS framework, with a sample of high 
energy data as given by S. Bethke [I], against the 3-loop APT coupling 
a E (Q 2 ) (I2.5P and its massive modification (I2.34p . Also shown in the figure 
is the common perturbative 3-loop coupling with IR singular behaviour that 
is ruled out by the data. As can be seen, the BS-APT theoretical scheme 
allows a rather satisfactory correlated understanding of very high and rather 
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Chapter 4 



Conclusive remarks 



To summarize, as discussed in Sec. 1, straightforward application of the RG 
method to perturbative expansions in QCD eventually gives rise to unphys- 
ical singularities in both the running coupling and the physical observables, 
in the low energy domain. The Landau singularities severely complicate the- 
oretical analyses of hadron dynamics. 

Then, an useful prescription to eliminate unphysical singularities is provided 
by the APT approach, outlined in Sec. 2, which is based upon the analycity 
requirement of the QCD coupling. This prescription has been used in previ- 
ous works in the framework of a Bethe-Salpeter like formalism to compute 
the meson spectrum. 

In this work, the point of view has been then somewhat reversed, namely, 
as discussed in Sec. 3, the comparison between the calculated meson spec- 
trum within the BS framework and the experimental data is used, in order to 
extract information on the infrared behavior of the QCD running coupling. 

The method consists in solving the eigenvalue equation for the squared mass 
operator as given by Eq. fl3.2|) . obtained by a three dimensional reduction 
of the original BS equation. The relativistic potential U then follows from 
a proper ansatz (13.11) on the Wilson loop to encode confinement, and is the 
sum ( 13. 2D of a one-gluon-exchange term Uoge and a confining term Ucorrf • 
The strong coupling occurring in the perturbative part of the potential needs 
to be IR finite since its argument has to be identified with the momentum 
transfer in the qq interaction, and this typically takes values down to few 
hundreds MeV. 
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As a first step, the common perturbative running coupling a s has then been 

replaced by its 1-loop analytized counterpart a E Eq. (12. 6ft as derived by the 

APT algorithm, which avoids the hurdle of the spurious singularities in the 

IR region [23], the whole theoretical scheme being extensively discussed in 

Sec. 2. 

An effective QCD scale Ajjl 3 = 193 MeV has been imposed in order to 

reasonably reproduces the 3-loop APT coupling, normalized at the Z boson 

mass along with world average [2] (i.e., A„ =5 = 236 MeV and A„ =3 = 

417 MeV by continuous threshold matching), specifically in the neighbohood 

of 0.5 GeV, which is the typical scale of most of the imput states. 

Then, it has been taken advantage of the BS results for the meson spectrum, 

both in the light and heavy quark sector, to deduce within this framework 

the fixed coupling value for each state that exactly matches the theoretical 

and experimental mass. 

A key point is finally the comparison of such experimental determinations 

a exp with thg "Euclidean" APT coupling a E (Q 2 ) and its further "massive" 

development. 

The results are twofold. On the one hand, as expected on the basis of the rea- 
sonable agreement between the theoretical and the experimental meson data, 
the 3-loop analytic coupling remarkably fits the extracted points a^ xp (<3 2 ) 
from 1 GeV down to 200 MeV, within the evaluated errors, both theoretical 
and experimental. This confirms and yields a quantitative estimate of the 
relevance of the APT to IR phenomena down to 200 MeV. 
Besides, below this scale, the experimental points exhibit a tendency to devi- 
ate with respect to the APT curve. Despite the sizable errors, this could give 
a hint on the vanishing of a s (Q 2 ) as Q — > in concert with some results from 
lattice simulations [T7] . It should be emphasized, however, that the existence 
of a finite IR limit of a s (Q 2 ) , lower than the universal APT freezing value, 
can not be ruled out by these data. 

Neverthless, a vanishing QCD coupling at zero momentum transfer can be 
theoretically discussed in the framework of a recent "massive" modifica- 
tion [9U] of the APT formalism, which takes into account effects of a finite 
threshold in the dispersion relations. Since in the extremely low Q region 
confinement forces play the dominant role, the overall reasonable agreement 
between the "massive" APT model and the results of the BS formalism would 
suggest a relation between the linear potential, arising from the area term in 
the ansatz (13.11) . and the thresholds effects in the analytic properties of the 



60 



QCD coupling, that needs to be further investigated. 
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Appendix A: Second order Bethe-Salpeter for- 
malism 

In the QCD framework a second order four point quark-antiquark function 
and the full quark propagator can be defined as 

H 4 (x u x 2 ; 2/1, y 2 ) = --Tr co i or (Ai(:n, y t ; A)A 2 (y 2 , x 2 ; A)) (4.1) 

and 

H 2 (x -y) = -J=Tr color (A(x, y; A)) , (4.2) 

where 

(f[A]) = J DAM F [A] J s ^f[A] , (4.3) 

M F [A] = Det n| =1 [l + gYA^iYjdju - mf)' 1 ] and A(x, y; A) is the second 
order quark propagator in an external gauge field. 

The quantity A is defined by the second order differential equation 

(£>„£>" + m 2 CWI - l -g a^F^)A(x, y- A) = -5\x - y) , (4.4) 

(aW = ip^y] an( j £)^ — Q^ _)_ igA^) and it is related to the corresponding 
first order propagator by S(x, y\ A) = [i^ v D v + m curr )A(x, y\ A) , m curr being 
the so-called current mass of the quark. 

The advantage of considering second order quantities is that the spin terms 
are more clearly separated and it is possible to write for A a generalized 
Feynman-Schwinger representation, i.e., to solve Eq. (14.41) in terms of a quark 
path integral [24"l 194"] . Using the latter in (14.11) or (14.21) a similar representa- 
tion can be obtained for H 4 and H 2 . 

The interesting aspect of this final representation is that the gauge field 
appears in it only through a Wilson line correlator W. In the limit x 2 —* x\, 
2/2 - ► Vi or V — > % the Wilson lines close in a single Wilson loop T and if T 
stays on a plane, i In W can be written according to (13.11) as 



i In W = ^a s I dz» I dz v 'D^{z - z') + (4.5) 

a I dz° I dz°'5(z° - z°')\z - z'| f d\{l - [\$j± + (1 - A)^] 2 }" . 
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The area term here is written as the algebraic sum of successive equal time 
strips and dz± = dz — (dz ■ r)r/r 2 denotes the transversal component of dz. 
The basic assumption now is that, in the center of mass frame, (14.51) re- 
mains a good approximation even in the general case, i.e., for non flat curves 
and when x 2 7^ #i, 2/2 7^ Hi or y ^ x. Then, by appropriate manipulations 
on the resulting expressions, an inhomogeneous Bethe-Salpeter equation for 
the 4-point function H±(x\, x 2 ; y±, y 2 ) and a Dyson-Schwinger equation for 
H2{x — y) can be derived in a kind of generalized ladder and rainbow approx- 
imation respectively. This should appear plausible, even from the point of 
view of graph resummation, for the analogy between the perturbative and the 
confinement terms in (14. 5p . In momentum representation, the corresponding 
homogeneous BS-equation becomes 

, /T s . f d 4 u f (, 1„ k + u 1_ k + u 

iFik > = -y (2^m*~" ; 2 p+ ~ 2 p ~~ 

x tff Qp + fc) a" *,(«) a" fff (~P + *j , (4.6) 

where a = 1; a, b = 0, \xv\ the cm. frame has to be understood, P = 
(mB,0); $p(fc) denotes the appropriate second order wave function, that in 
terms of the second order field 4>(x) = (i , ~f fl D fl + m curr )~ 1 ip(x) can be defined 
as the Fourier transform of (O|</>(|) , 0(— |)|P) . 

Similarly, in terms of the irreducible self-energy, defined by H 2 (k) = i(k 2 — 

m curr) _1 + i(k 2 — m^ rr ) _1 i T(k) H 2 {k) , the Dyson-Schwinger equation can be 

written 

dH » /, , k + l k + l\ n - r /7X b 



T(k) = j j^Ia^k-l;—,—^!!^)^. (4.7) 

The kernels are the same in the two Eqs. (14.61) and (14. 7p . consistently with 
the requirement of chiral symmetry limit [22] , being given by 

Io;o(Q;p,p') = 16<K^a s p a p'PD a p(Q) + 



4a jd 3 Ce-^<\C\e(p )e(p' ) J d\{p 2 p' 2 - [\p' Q p T + (1 - X)p p' T } 2 } 



2 



'0 



~f«'~ < **ffiJFfr 
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U,p*(Q;p,p') = -^a sP a (5PQ a - 8^Q p )b aP {Q) + 

/ ,3/- -iQ< >>pPcr ~ SoVp , i x 

+a J d(e pa ww^ eM 

W(0;p.p') = f 50. W» - KQ f W,Qc - kq„)d^(Q) , (4.8) 

where in the second and in the third equation ( = has to be understood. 
Notice that, due to the privileged role given to the cm. frame, the terms 
proportional to a in ( 14. 8ft formally are not covariant. 

In fact, it can be checked that T(k) can be consistently assumed to be spin 
independent and Eq. (14.71) can be rewritten in the simpler form 

with 

R(k,l) = A^a s [(k + l)'"(k + l) v D t , l/ (k-l)+ 

(k - l) v (k - l) v D£(k -l)-(k- lT(k - lYD^k - I)] + 

+a J dhe-^-^r(k + l ) 2 Jl - J^f , (4-10) 

k^ and lj_ denoting as above the transversal part of k and 1. Eq. ( 14. 9 j) 
can be solved by iteration resulting in an expression of the form T(/c 2 ,k 2 ), 
since (14.101) is not formally covariant. Then the constituent (pole) mass m is 
defined by the equation 

m 2 - m 2 CUTT + T(m 2 , k 2 ) = (4.11) 

and the dependence on k 2 , being an artifact of the ansatz (13. ip . is eliminated 
by extremizing m(k 2 ) in k 2 . 

The 3-dimensional reduction of Eq. (14.61) can be obtained by a usual proce- 
dure of replacing H 2 (k) with i(k 2 — m 2 )~ l and I a & with its so-called instan- 
taneous approximation J^ st (k, u) . In this way, one can explicitly integrate 
over uq and arrive to a 3-dimensional equation in the form of the eigenvalue 
equation for a squared mass operator Eq. (13. 2p . with 



(k|[/|k) = (2^V^^7 Iab (k ' k) V^« 12 ' ( } 
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Finally by using Eq. (14.81) one obtains Eqs. (13.3113. 51) . 

Alternatively, in more usual terms, one could look for the eigenvalue of the 
mass operator or center of mass Hamiltonian -£/cm = M — M + V with 
V defined by M V + VM + V 2 = U. Neglecting term V 2 the linear form 
potential V can be obtained from U by the replacement \ W2 >^ w i w 2> _^ 

— -, l The resulting expression is particularly useful for a comparison 

with models based on potential. In particular, in the static limit V reduces 
to the Cornell potential 

4 ry 

V stat = --- + ar. (4.13) 

3 r 

Note that it is necessary to introduce a cut-off B in Eq. (14. 9p . As a con- 
sequence the constituent mass turns out to be a function of the current 
mass and of B, m = m (m CUTT , B). Then if one uses a running current mass 
Tncurr(Q 2 ) , & running constituent mass m{Q 2 ) is obtained, as it has been done 
in [95]. However the singular expression used there 



v 7o/2/3o 

m - m = ™ { Q yA 2 HQ 2 /A 2 ) ) (414) 

is not consistent with Eq. (12. 6ft . and if a more consistent assumption is taken, 
e.g., 

m curr (g 2 )=m(4 1) (g 2 )J , (4.15) 

or the other resulting from the analytization of a similar expression with 
a E (Q 2 ) replaced with the ordinary perturbative c^ (Q 2 ), the dependence 
of m on Q 2 is strongly reduced. For this reason even the light quark mass is 
here treated as a constants to be adjusted with the the dataj. 



^n this way the only role that is left to the DS equation is to justify the difference 
between the constituent and the current masses. 
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Appendix B: Numerical results 



The tables below display the complete set of results as explained in Sec. 3.5. 



The values of all the parameters are a = 0.18 GeV 2 , 
= 196 MeV (q = u,d), m s = 352 MeV, m c 



A 



(l,eff) 



193 MeV 



m 



1.516 GeV and 

m& = 4.854 GeV . Meson masses are given in MeV. The last column dis- 
plays the experimental coupling a| xp (<5a) with the theoretical error Anlo 
due to the next-to-leading order terms neglected, the theoretical error Ar 
from the half width T/2 and the experimental error A exp respectively. 



f Center of gravity masses of the incomplete multiplets estimated in analogy with other 
multiplcts. 

Table 4.1: qq (q = u, d) 



States (MeV) 

1 3 S! 

1ASS 
2% 
2% 
2ASS 



in 



cxp 



™ t h 



Q 



,.( :! > 



^ x p ± A 



NLO 



± A r ± A, 



cxp 



7T 

P (770) 



tt(1300) 
p(1450) 



134.9766 ±0.0006 
139.57018 ±0.00035 
775.5 ± 0.4 

638 

1300 ± 100 

1459 ± 11 

159 



138 136 401 0.522 



749 
613 

1223 448 
1363 
139 



0.502 



0.534±0.122±-±- 
0.517 ±0.122 ±0.048 ± - 

0.451 ±0.114 ±0.152 ±0.081 
0.427 ± 0.114 ± 0.062 ± 0.010 



1'P1 


bi(1235) 


1229.5 ± 3.2 


1234 209 0.637 0.688 ± 0.155 ± 0.124 ± 0.006 


FDs 


tt 2 (1670) 


1672.4 ±3.2 


1595 144 0.701 0.544 ± 0.151 ± 0.364 ± 0.009 


Table 4.2: ss 


States 


(MeV) 


^exp 


n ( 3 ) 


af p ± A NLO ± A r ± A cxp 


l a Si 
2 3 Si 


(/.(1020) 
(/.(1680) 


1019.460 ± 0.019 
1680 ± 20 


1019 418 0.514 
1602 454 0.500 


0.525 ±0.098 ±0.002 ±- 
0.435 ± 0.096 ± 0.068 ± 0.019 



^i 


/ii(1380) 


1386 ± 19 


3 P 2 


/2(1525) 


1525 ±5 ] 


3 Pi 


/i(1510) 


1518 ±5 I 1521 


3 Po 


/o(1500) 


1507 ±5 J 



1472 216 0.631 0.824 ± 0.098 ± 0.083 ± 0.032 



1484 



0.603 ± 0.098 ± 0.070 ± 0.009 



T^o 



1 3 F 4 
1 3 F 3 
1 3 F 2 



^(1870) 



1842 ± 8 



1807 149 0.695 0.658 ± 0.079 ± 0.318 ± 0.023 



2165 f 2070 118 0.733 0.452 ± 0.064 ± 0.137 ± 0.024 



/ 2 (2150) 



2156 ±11 
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Table 4.3: qs (q = u, d) 



States 



(McV) 



m 



cxp 



TOth 



.,(3) 



af xp ± Anlo ± A r ± A, 



cxp 



1% 




497.648 ±0.022 1 ... 
> 495 
493.677 ±0.016 J 


491 


409 


0.518 


0.529±0.122±-±- 


l 3 Si 


( K* (892)° 
\ AT*(892) ± 


896.00 ±0.25 1 
891.66 ±0.26 j ™- 11 


887 






0.526 ±0.122 ±0.017 ±- 


1ASS 




398 


396 








2 3 Si 


A'*(1410) 


1414 ± 15 


1485 


451 


0.501 


0.571±0.117±0.102±0.013 



1^P1 



A"i(1270) 
K*(U30)° 
J Ff|(1430) ± 
A"i(1400) 



1272 ± 7 
1432.4 ±1.3 1 
1425.6 ±1.5 
1402 ± 7 



1355 213 0.634 0.820 ± 0.129 ± 0.081 ± 0.012 



V 1417.7 1367 



0.583 ±0.129 ±0.133 ±0.007 



1 3 P 


A o *(1430) 


1414 ±6 












1 3 D 3 
1 3 D 2 
l 3 Di 


iff (1780) 
A 2 (1770) 
if*(1680) 


1776 ±7 
1773 ±8 
1717 ±27 I 


>1763 


1712 


150 


0.694 


0.617 ± 0.113 ± 0.273 ± 0.031 


1 3 F 4 
1 3 F 3 
1 3 F 2 


-fQ(2045) 
A" 3 (2320) 
^(1980) 


2045 ±9 1 
2324 ±24 ' 
1973 ±25 J 


> 2121 


1973 


116 


0.736 


0.248 ± 0.095 ± 0.413 ± 0.071 
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Table 4.4 


: cc 






States 


(MeV) 


^-exp 


wth 


Q 


(3) 


a«P ± A NLO ± A r ± A exp 


1% 


Vc(lS) 


2980.4 ± 1.2 


2980 


561 


0.464 


0.467 ±0.025 ±0.008 ±0.001 


l 3 Si 


JMis) 


3096.916 ± 0.011 


3097 






0.467 ±0.025±-±- 


1ASS 




117 


118 








2% 


Tjc(25) 


3638 ± 4 


3595 


500 


0.483 


0.446 ± 0.023 ± 0.007 ± 0.004 


2 3 Si 


i>(2S) 


3686.093 ± 0.034 


3653 






0.455 ±0.023±-±- 


2ASS 




48 


58 








3 3 S! 


^(4040) 


4039 ± 1 


4030 


483 


0.489 


0.485 ± 0.022 ± 0.049 ± 0.001 


4 3 Si 


^(4415) 


4421 ± 4 


4337 


474 


0.492 


0.384 ± 0.022 ± 0.042 ± 0.006 


l^i 


h c (lP) 


3525.93 ± 0.27 


3532 


269 


0.592 


0.631 ± 0.012 ±- ±- 


1 3 P 2 


Xc2(lP) 


3556.20 ±0.09 ] 










l 3 Pi 


Xcl(lP) 


3510.66 ±0.07 \ 3525.3 


3537 






0.640 ±0.012 ±0.002 ±- 


l 3 Po 


Xco(lP) 


3414.76 ±0.35 J 










2 3 P 2 


Xc 2 (2P) 


3929 ± 5 ] 










2 3 Pi 




\ 3915+ 


3929 


274 


0.589 


0.644 ± 0.013 ± 0.027 ± 0.009 


2 3 P 


X(3872) 


3871.2 ±0.5 1 










1 3 D 3 
1 3 D 2 




I 3820+ 


3822 


190 


0.654 


0.707 ± 0.008 ± 0.030 ± 0.006 


l 3 Di 
2 3 D 3 
2 3 D 2 


^(3770) 


3771.1 ±2.4 J 












I 4183+ 


4150 


198 


0.646 


0.606 ±0.009 ±0.132 ±0.008 


2 3 Di 


^(4160) 


4153 ±3 J 











Quantum numbers of h c {lP) and X(3872) mesons are not well established. 
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Table 4.5: bb 



States 



2 3 Si 
3 3 S! 
4 3 Si 
5 3 Si 
6% 



(MeV) 



mth 



.X 3 ) 



^ xp ± A 



NLO 



± Ar ± A c 



i 3 Pi 

l 3 Po 
2 3 P 2 
2 3 Pi 
2 3 P 



T(15) 

T(2S) 

T(35) 

T(45) 

T(10860) 

T(11020) 



9460.30 ± 0.26 

10023.26 ± 0.31 

10355.2 ± 0.5 

10579.4 ± 1.2 

10865 ± 8 

11019 ± 8 



9461 


951 


0.381 


9987 


630 


0.445 


L0321 


552 


0.466 


L0588 


517 


0.478 


L0820 


497 


0.484 


L1034 


506 


0.481 



0.378 ± 0.006 ±- ±- 

0.416±0.004±-±- 

0.433±0.003±-±- 

0.493 ± 0.003 ± 0.013 ± 0.002 

0.424 ± 0.003 ± 0.078 ± 0.011 

0.508 ±0.003 ±0.057 ±0.012 



X62(1P) 9912.21 ±0.26 ±0.31 

Xfei(lP) 9892.78 ±0.26 ±0.31 

Xbo(lP) 9859.44 ±0.42 ±0.31 

X b2 (2P) 10268.65 ±0.22 ±0.50 

Xfei(2P) 10255.46 ±0.22 ±0.50 

Xb (2P) 10232.5 ±0.4 ±0.50 



9899.87 



9880 387 0.528 



0.519±0.002±-±- 



10260.24 10231 343 0.549 0.524 ± 0.002 ± - ± 0.001 



Table 4.6: Light-heavy quarkonium systems 



States 



(MeV) 



''CXp 



m t h 



.,(3) 



af p ± A NLO ± A r ± A 



cxp 



qc 

l 3 Si 
1ASS 



D ± 

D a 

J D*(2010) ± 

D*(2007)° 



1869.3 ±0.4 
1864.5 ±0.4 
2010.0 ±0.4 
2006.7 ±0.4 
141 ± 1 



1867.7 



2008.9 



1843 


459 


0.498 


0.488 ±0.082 ± - 


-± 


2000 






0.499 ±0.082 ±- 


-± 


157 











1P2 
1P1 



qb 



l 3 Si 
1ASS 



£>^(2460) ± 
D^(2460) 
Di(2420) ± 
£>i(2420)° 



2459 ± 4 
2461.1 ±1.6 
2423.4 ±3.1 
2422.3 ±1.3 



2443 232 0.619 0.651 ± 0.051 ± 0.028 ± 0.005 



B* 



5279.0 ±0.5 
5279.4 ±0.5 
5325.0 ± 0.6 
46 ± 3 



5279.1 5246 516 0.478 

5311 
64 



0.456±0.036±-± 
0.471 ±0.036±-± 
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Table 4.7: Light-heavy quarkonium systems 



States 


(MeV) 


77^cxp 


mth 


Q 


(3) 


< p ± Anlo ± Ar ± A cxp 


sc 
1% 

l 3 Si 
lASS 




1968.2 ± 0.5 

2112.0 ± 0.6 

144 ± 1 


1959 
2109 
149 


472 


0.493 


0.494 ±0.055±-±- 
0.497 ±0.055±-±- 


1P 2 
lPi 


£» a2 (2573) ± 

D s i(2536) ± 


2573.5 ±1.7 1 
2535.35 ±0.34 ±0.5 J 


2545 


236 


0.616 


0.626 ± 0.036 ± 0.009 ± 0.002 


sb 

1% 

lASS 


s 8 * 


5367.5 ± 1.8 

5412.8 ± 1.7 
47 ± 4 


5343 

5408 
65 


535 


0.472 


0.457±0.024±-± 0.001 
0.473 ± 0.024 ±- ±0.001 


IP 


5^(5850) 


5853 ± 15 


5830 


255 


0.602 


0.592 ± 0.013 ± 0.042 ± 0.027 
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Part II 

Progress in the determination 
of the gluon polarization 



Chapter 5 

Polarized parton distributions 
from DIS 



One of the fundamental properties of the nucleon structure is the spin dis- 
tribution among its quark and gluon constituents. The spin structure of the 
nucleon has been extensively investigated by inclusive polarized lepton scat- 
tering off polarized protons and neutrons, since the discovery of the EMC 
spin effect [TJ. Despite the naive expectation that the nucleon spin is car- 
ried by quarks, the experimental results revealed that only a small fraction 
is actually due to quarks (the so-called "proton spin crisis"). Considerable 
efforts, both experimentally and theoretically, have subsequentely gone into 
understanding where the remaining fraction of the nucleon spin resides (see 
e.g. Ref. [2] for a recent review). 

To this end, the determination of the first moments of the polarized parton 
distribution functions (pdfs) is necessary. Attention is currently devoted to 
the reconstruction of pdfs at all values of the Bjorken x and the momentum 
transfer Q 2 . Of special interest is the determination of the gluon polariza- 
tion, in order to evaluate the total contribution carried by the gluons to the 
nucleon spin. 

The possibility of inferring the polarized gluon density from scaling violations 
has been extensively studied. However, phenomenological analyses of the 
current inclusive DIS data only yield weak constraints on the gluon first 
moment, and leave the x-shape completely undetermined. This issue has 
been thoroughly investigated in the present work, by performing a NLO 
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analysis in perturbative QCD of world data from inclusive DIS experiments, 
as discussed in detail in Sec. 6. 

Actually, some light may be shed on the gluon content of the nucleon by 
observing specific events that receive leading contributions from gluon initi- 
ated subprocesses, that is, by directly measuring the polarized gluon density. 
In particular, as discussed in Sec. 7, in polarized fixed target experiments, 
the gluon polarization can be probed through the photon-gluon fusion pro- 
cess, by selecting open-charm events or hadron pairs with high transverse 
momentum, and both strategies are currently used at the COMPASS exper- 
iment (CERN). Indeed charmed mesons in the final states provide a clean 
tag of photon-gluon fusion, since there are essentially no competing processes 
and a clean perturbative analysis is possible. Therefore, a phenomenological 
study of open-charm photoproduction from deep-inelastic muon scattering 
off nucleons, with longitudinally polarized beam and target, has been also 
performed, with reference to upcoming results on experimental asymmetries 
from COMPASS. The results are presented in Sec. 7 and, specifically, the 
constraints that can be obtained on the polarized gluon density at the COM- 
PASS kinematics are discussed. 

In what follows the main theoretical tools are briefly reviewed. Specifically 
Sec. 5.1 recalls the basic ingredients for the study of polarized lepton-nucleon 
scattering, i.e. the polarized structure functions g\ and g^. In Sec. 5.2 the ex- 
pectations of the naive quark-parton model are summarized, whereas Sec. 5.3 
is devoted to the main aspects of the perturbative QCD analysis, that is, the 
relations between structure functions and polarized parton distributions. In 
Sec. 5.4 the way experimental information on the relevant structure function 
gx are extracted from the experimentally measured quantities is described, 
and, finally Sec 5.5 sketches how target mass corrections are taken into ac- 
count in the present analysis. 



5.1 Structure functions and their moments 

Polarized structure functions are the form factors which parametrize the 
cross section for deep-inelastic scattering of polarized leptons off a polarized 
hadronic target. If the momentum transfer Q 2 involved is much smaller 
than the Z boson mass, the antisymmetric part of the hadronic tensor can 
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be expressed in terms of two spin-dependent structure functions, g\ and g 2 , 
i.e. |31I5 

iWF = -^ J d A xe iqx <p,s\J^(x)r ] (0)\p,s > 



ime^qp 



-9i(x,Q 2 ) + 
pq 



s a pq - p a qs 
(pq) 2 



92(x,Q 2 



(5-1) 



where p M , m and s M are, respectively, the target four-momentum, mass, and 
spin (normalized as s^s^ = — 1). Neglecting weak interaction effects, J M in 
Eq. (15.11) is the electric current. The deeply inelastic region is identified by 
the Bjorken limit 



Q 2 = ~<? 



oo , with x 



2pq 



fixed 



(5.2) 



in which the invariant mass of the hadronic system in the final state is much 
larger than the nucleon mass, namely 



W 2 = m 2 + Q 2 



x 



x 



> VTi 



In the case when the nucleon spin s is purely longitudinal, that is, 



m 



Sl = A p q 

pq 



A = ±- 



1 



w-a/I + 4m 2 x 2 /Q 2 
the antysimmetric tensor Eq. (15.11) can be written as 

4m 2 x 2 



W% = Xme^ pa 



QpPo 
pq 



9i(x,Q 



Q 2 



92(x,Q' 



(5.3) 



(5.4) 



(5.5) 



In the Bjorken limit the product Am — > 1, thus leaving a finite residue as 

m — > 0. Furthermore, the factor 4m 2 x 2 /Q 2 vanishes, and the structure 

function g 2 decouples. As a result only g± is asimptotically relevant. 

In case of transverse polarization both structure functions equally contribute, 

but the whole cross section is strongly suppressed by an overall target mass 

factor. 



The light-cone expansion of the current product in Eq. (15. ip implies that the 
first moment of the relevant structure function g\ for a nucleon target, at 
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leading twist, is given by [5] 



Ti(Q 2 



dx gi(x, Q z 



^GiCf)* 



i=l 



(5.6) 



where Oj are given by matrix elements of the axial current for the i-th flavor, 
namely 

mdiSfj, =<p, s|^i7 M 75^»b,s >, (5.7) 

e, is the electric charge, and Ci(Q 2 ) are perturbatively calculable coefficient 
functions. Then, assuming that only the three lightest flavors are activated, 
Eq. (J5.6P can be recast as 



ITW) = y 2 



C NS (Q 2 ) ±o 



3 as 



^c s (g 2 )a 



(5.i 



where the plus (minus) sign refers to a proton (neutron) target, and the 
singlet and nonsinglet matrix elements are defined by 



a 
a 3 



a d - 2a s 



(5.9) 



a u + a d + a s 

respectively Cs(Q 2 ) and Cns(Q 2 ) are Q 2 -dependent singlet and nonsinglet 
QCD coefficient functions; to order a s , in the MS scheme e.g., one has 



O.s 



C S (Q 2 ) = C m (Q 2 ) = 1 - - + 0(a 



7T 



(5.10) 



whereas Cg 7^ C*ns at higher orders [6]. QCD corrections turn out to be 
sizable in the Q 2 -Tsmge of the present experiments and are thus relevant for 
the comparison with the data [7]. 

The singlet and nonsinglet components of r\ can be then extracted using 
Eq. (15.81) . Indeed, the triplet and octet currents are conserved and therefore 
scale invariant, and the related matrix elements can be derived from any other 
process. The triplet axial charge is equal to the axial coupling measured in 
nucleon /3-decay 

9A 

9v 



a 3 



F + D = 1.2601 ±0.0025 



(5.11) 



where F and D are weak hyperon decay constants in flavor SU(3) symme- 
try [HI [9] . The octet matrix element is then given by 



a s = 3F- D = 0.588 ± 0.033 



(5.12) 



S3 



On the other hand, the singlet axial charge a® is scale dependent due to 
the anomalous nonconservation of the singlet axial current [TO] , and can be 
determined from a measurement of the proton first moment, by rewriting 
Eq. ([EH]) as 

Cs(Q 2 )ao(Q 2 ) = 9r?(Q 2 ) - Ic NS (Q 2 )(3F + D) . (5.13) 

Since the two terms on the right hand side are roughly of the same order 
(typically Ti ~ 0.1), the value of a^ arises from a large cancellation between 
them (see e.g. [1]). Present data indeed indicate that a$ is compatible with 
zero, but values as large as ao(10 GeV 2 ) ~ 0.3 are not excluded. 

Finally, a fundamental prediction of the theory is the Bjorken sum rule 

r?(Q 2 )-r?(Q 2 ) = ic NS (Q 2 )a 3 , (5.14) 

o 

that was originally derived from current algebra and isospin symmetry [TT] . 
Eq. ( I5.14p relates the integral over all x of the nonsinglet polarized structure 
function g\ — g™, at fixed Q 2 , to the well- measured /3-decay coupling a 3 . A 
comparison with the data, using the computed Cns(Q 2 ) thus allows a direct 
test of isospin in this channel, as well as the predicted scale dependence. 
In addition, it is possible to use the measurements to extract a relatively 
accurate determination of the strong coupling a s (Q 2 ) 



5.2 The parton model predictions 

In the naive quark-parton model the spin-dependent structure function g\ can 
be simply expressed as the charge-weighted difference between momentum 
distributions for quark helicities aligned parallel (q 1 ') and antiparallel (g^) to 
the longitudinally polarized parent nucleon 



gi(x) = -J2e 2 Aq + (x), Aq + (x) = Aq l (x) + Aq l (x). (5.15) 



i=l 



The polarized quark distributions are defined as 

Aft(x) = qj(x) - q\(x) , (5.16) 



84 



and similarly for antiquarks. Assuming n/ = 3, one can define the flavor 
singlet and nonsinglet combinations of polarized quark densities, namely 



AS = Au+ + Ad + + As+ 

Ag 3 = Au + - Ad + 

Aq 8 = Au + + Ad + - 2As + . 

Then the structure function g^ can be cast into the form 



p-^i \ 



< e z > 



AE(x) ± -Aq 3 (x) + -Aq 8 (x) 



(5.17) 



(5.1* 



where the plus (minus) sign stays for a proton (neutron) target, and the av- 
erage charge is < e 2 >= 2/9 for rif = 3. 

In the naive parton model gluons do not contribute to the nucleon spin, and 
the first moments of singlet and nonsinglet quark densities are straightfor- 
wardly related to the matrix elements of axial currents, such that 



a 



/ dxAT,(x) , a 3 = / dxAqz(x) and a 8 = dxAq, 
Jo Jo Jo 



•s(x). (5.19) 



In this context, the first moments of Aq^~(x) can be interpreted as the con- 
tribution of the quark flavor % to the nucleon spin, and thus the singlet axial 
charge a^ as the net quark helicity. 

Then, under the assumption that the strange sea in the nucleon is unpolar- 
ized, namely 



dxAs + (x) 



< 



dxAq + (x) 



(5.20) 



q=u,d 



the quark helicity coincides with the nonsinglet (octet) axial charge, a^ ~ a§. 
Thus, from Eq. (15.121) one roughly has 



a ^ 0.6 



(5.21) 



whereas the value extracted from the data (see Eq. (15.131) ) is compatible with 
zero. The Ellis- Jaffe sum rules [13] then give the parton model predictions 
for the proton and neutron first moments 



rr = ±t^3 + ^ 



12 



36 



(5.22) 



N5 



that are not consistent with the measured values, at the level of several stan- 
dard deviations. Thus, against the Ellis- Jaffe prediction and intuition, quarks 
seem to carry a very small fraction of the total nucleon helicity. Actually, it 
is only at the naive parton level that the first moment of the singlet part of 
the structure function g\ corresponds to the total helicity fraction carried by 
quarks. 



5.3 Polarized structure functions and pdfs in 
QCD 

When including QCD corrections, also the gluon density contributes and 
the structure function g 1 is related to the scale-dependent polarized quark 
and gluon distributions by the convolution with appropriate coefficient func- 
tions pn dg 

9i(x, Q 2 ) = — |— [C NS g> Ag NS + C s ® AS + 2n f C g <g> Ag] , (5.23) 

where < e 2 >= nj 1 Y^a=\ e h n f being the number of active flavors with 
electric charge e$, and ® denotes the convolution product with respect to x 

f®g= f-f(-}g{y)- (5.24) 

Jx y \yJ 

In Eq. (15.231) Ag(x, Q 2 ) is the polarized gluon distribution, whereas the sin- 
glet and nonsinglet quark distributions are defined as 

AS(x,g 2 ) = ^(Ag, + Ag i ) (5.25) 

"Z, / e 2 \ 

Ag NS (x, Q 2 ) = J2 (7^7 " l ) ( A( H + A ?i) . 



«=i 



where Aq t and Agj are the scale-dependent quark and antiquark polarized 
densities of flavor i, defined as in Eq. (j5.16p . 

The perturbative part of the Q 2 -dependence of the polarized parton densities 
is given by the Altarelli-Parisi evolution equations [T9] . The singlet quark 



8(3 



and gluon distributions mix according to 



8 ( AS \ ««(*) ( ^ p l 2n/AP a 5 ^ x> ( AS 



flln. V A» ; - 2. V ^! ^//' J ® U 9) ' <5 ' 26) 

while the nonsinglet quark evolves independently as 

) -Ag NS = ^APf ® Ag NS . (5.27) 



dint ™° 2tt 

In Eqs. ( I5.26P and (J5.27P AP(x, Q 2 ) denote the spin-dependent splitting func- 
tions, and t = ln(Q 2 /A 2 ) with A the QCD scale parameter. The factorization 
and renormalization scales are both taken equal to Q 2 , so that all the scale 
dependence appears through t. 

Eqs. (15.261) and (I5.27P are valid for all orders of perturbative QCD; the coef- 
ficient functions C and the evolution kernels AP may each be expanded in 
powers of a s . At NLO 

C(x, a s ) = C {0 \x) + ^-C {1 \x) + 0{a 2 s ) (5.28) 

2ir 

AP(x, a.) = AP(°'(x) + ^IapW(x) + 0(a 2 s ) . (5.29) 

Z7T 

It is also convenient to introduce anomalous dimensions 

7(71, a s )= / dxx n ~ 1 AP(x,a s ) (5.30) 

Jo 

i.e. the Mellin transforms of the splitting functions. One can analogously 
define moment-space coefficient functions C(n, a s ) and parton distributions 
Ag NS (n,Q 2 ), AS(n,Q 2 ) and Ag(n,Q 2 ). 

In accordance with the partonic picture, in Eq. ( 15. 28ft C*ns( x ) = C$ (x) = 
5(1 — x), while C g (x) = 0, so that at order a° g\ decouples from Ag, 
and is just a linear combination of polarized quark distributions, whose Q 2 
dependence is entirely specified by Eqs. ( 15.261) and ( I5.27P (the leading order 
splitting functions are calculated in Ref. [19]). 

Beyond leading order, splitting functions and coefficient functions are no 
longer universal, hence even though the scale dependence of the observable 
structure function gx is determined uniquely, at least up to higher order cor- 
rections, its separation into contributions due to quarks and gluons is scheme 
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dependent and thus essentially arbitrary. The NLO coefficient functions 
may be modified by a change of the factorization scheme which is partially 
compensated by a corresponding change in the NLO anomalous dimensions, 
hence both are required in a full NLO order computation. The complete set 
of NLO coefficient functions and splitting functions can be found in Ref . [6] . 
The modified minimal subtraction (MS) factorization scheme is commonly 
used in the analysis of polarized parton distributions. In the MS scheme 
the first moment of the gluon coefficient function vanishes, the gluon density 
does not contribute to the first moment of g\ and the scale dependent singlet 
axial charge is thus equal to the singlet quark first moment 

a (g 2 )=AS M s(l,g 2 ). (5.31) 

Alternatively, one can define a different scheme, the Adler-Bardeen (AB) 
scheme [TOj, such that the singlet quark first moment is conserved at all orders 
and can be identified with the total quark helicity, whereas the polarized 
gluon density is defined as in the MS scheme. In the AB scheme the gluon 
polarization directly contributes to the singlet axial charge (and thus to the 
first moment of g\ ) , which is now written as 

a (Q 2 ) = AE A b(1) - n f ^^-Ag(l, Q 2 ) . (5.32) 

Eqs. (15.311) and fj5.32j) then yield the relation between the first moments of 
the singlet quark distributions in the two schemes. They differ by the prod- 
uct a s (Q 2 )Ag(l, Q 2 ), which is due to the anomalous nonconservation of the 
singlet axial current [TO], and is scale invariant at leading order; this implies 
that Ag(l,Q 2 ) increases as l/a s (Q 2 ) with Q 2 , and the gluon contribution 
in Eq. (15.321) is not asymptotically suppressed by powers of a s . As a result, 
this scheme dependence does not vanish at large Q 2 , and the definition of 
the singlet quark first moment is therefore maximally ambiguous. 

In the AB scheme, which will be adopted in the following phenomeno logical 
analysis (see Sec. 6), the first moment of the gluon coefficient function at 
NLO reads \m\ 

C g (l,a s ) = -^ + 0(a 2 ), (5.33) 

and the first moment of the structure function g\ is given at NLO by 

Y X {Q 2 )= f dx 9l (x,Q 2 )= (5.34) 

Jo 
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< e 2 > 



1 - 2i) (Ag NS (l, Q 2 ) + AS(1, Q 2 )) - n,£.Ag{l, Q 2 ) 

IT / Z7T 



2 

Different factorization schemes were also discussed in Refs. [20], and the 
dependence of the results of phenomenological analyses of spin-dependent 
pdfs on the choice of the scheme was studied. 

It should be finally noted that a number of theoretical models attempt to ex- 
plain the quark helicity distribution within the nucleon and the unexpected 
smallness of the singlet axial charge. A possible interpretation is to assume a 
large and negative contribution from the strange sea polarization As + . This 
suppression of the axial charge might be explained by invoking nonperturba- 
tive mechanisms based on instanton-like vacuum configuration |14J . In this 
case As = As. Another scenario is possible, where the smallness of the sin- 
glet axial charge is due to the intrinsic strangeness, i.e. the C-even strange 
combination is large, but As differs significantly from As. Specifically, in 
the Skyrme models of the nucleon [IB] the strange distribution (in particular 
its first moment) is large, while the antistrange distribution is much smaller 
and does not sizably contribute to the axial charge [T5j . 
Finally, another possible mechanism has been proposed to understand the 
nucleon spin structure (see [12] and references therein), which includes the 
gluon contribution. Indeed, in the Adler-Bardeen scheme, one can assume, on 
the basis of Eq. f)5.32p . that a cancellation between a large (scale-indepedent) 
singlet quark distribution and a large gluon polarization takes place. In this 
case As + <C Au + , Ad + , as expected in the parton model. 



5.4 Phenomenology of g\ 

Experimental information on the structure functions g±(x, Q 2 ) and g2{x, Q 2 ) 
are extracted from measurements of spin-asymmetries. Longitudinally po- 
larized leptons are scattered off a hadronic target that is polarized either 
longitudinally or transversely. The longitudinal (A\\) and tranverse (A±) 
asymmetries are formed by combining data taken with opposite beam helic- 
ity 

a U _ g-TT pi-* _ p.?-* 

Au = — m TT > A ± = —> r~ ■ ( 5 - 35 ) 

The symbols cr^ and a'^ denote the cross sections for the lepton-nucleon 
scattering with their parallel and antiparallel helicity states, respectively, 
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whereas cr - * and a^ represent the scattering cross sections for transversely 
polarized nucleon targets. 

The asymmetries A\\ and A±_ are related to the photon absorption cross 
section asymmetries Ai and A 2 by the relations 

A\\ =D(A 1 + V A 2 ) , A ± = d (A 2 - CAi) , (5.36) 

where 

A(x,g 2 ) = ^ P, A 2 (x,Q 2 ) = T ° T . (5.37) 

°l/2 + °3/2 (T l/2 + ^3/2 

Here, crT, 2 and <rj, 2 represent the absorption cross sections of virtual trans- 
verse photons for the total helicity of the photon-nucleon system of 1/2 and 
3/2 respectively; a LT denotes the interference term between the transverse 
and longitudinal photon-nucleon amplitudes. 

The factor D in the first of Eqs. (15.361) is interpreted as the depolarization 
of the photon with respect to the primary lepton beam 

D - y(i + i 2 y/2)(2-y) 

"y 2 (l+7 2 ) + 2(l + J R)(l- Z /- 7 V/4)' l ' } 

and depends upon the kinematic factors 

E' 2mx 

y = 1 — — and 7 = — -== (5.39) 

where y is the fractional energy lost by the lepton. The depolarization factor 
also depends on the ratio R(x, Q 2 ) of the cross sections for the longitudi- 
nally polarized photon to the transverse one, and is experimentally measured. 
Moreover, 

7(1 - y - -j 2 y 2 /A) 

" = (l+TV2)(l-»/2) - (5 ' 40) 

The kinematic factors in the second of Eqs. (I5.36P are given respectively by 



zyi-,,- 7 y/4 c = 7(i-y/2) (541) 

Then, the asymmetries A\ and A 2 can be related to the polarized structure 
functions g\ and g 2 by 

a f n2 N 9i{x,Q 2 )-^ 2 g 2 {x,Q 2 ) , 2 , l(gi(x,Q 2 )+g 2 (x,Q 2 )) 

(5.42) 
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where Fi(x, Q 2 ) is the unpolarized structure function of the nucleon. 
If both the longitudinal and transverse asymmetries, Au and A±, are mea- 
sured, one can extract both structure functions, gi(x, Q 2 ) and g 2 {x, Q 2 ) from 
experimental data with minimal assumptions, using the relations 



9i(x,Q 2 ) 
92(x,Q 2 ) 



F^x.Q 2 



A 



A, 



(l + 7 2 )(l + r/C) 

Fi(x,q 2 ) r/c 

(l + 7 2)(l + r / C) LU 



:i+70^-fo-7)^ 



i 



4l 



'/ 



7 






(5.43) 
(5.44) 



Actually, the transverse asymmetry is generally poorly determined, and in 
most cases only the longitudinal asymmetry A\\ is measured. In this case, 
the relevant structure function g^ can be expressed via the observable An and 
the unknown structure function g 2 , namely 



9i(x,Q* 



Fi(x,Q 2 ) A|, | 7 ( 7 -rj) 
(1 + 777) D 777 + 1 



g2(x,Q 2 



(5.45) 



An experimental determination of g\ thus relies on a theoretical assumption 
on g 2 . In Eq. (15.451) F 1 is expressed in terms of the unpolarized structure 
function F 2 , usually extracted from unpolarized DIS experiments, i.e. 



F^Q 2 



F 2 (x,Q 2 



2x(l + R(x,Q 2 )) 



7 



(5.46) 



Then, as discussed in Sec. 6, parametrizations for the unpolarized struc- 
ture function ^(x, Q 2 ) and the ratio R(x, Q 2 ) of longitudinal and transverse 
virtual photo-absorption cross sections must be used in order to extract ex- 
perimental information on the structure function g±(x, Q 2 ). 
Finally, the experimentally measured counting rate asymmetry A exp is related 
to the cross section asymmetry Ay by 



.4. 



cxp 



ftPtP b A\ 



(5.47) 



where P& is the beam polarization, P t the polarization of the target nucleon, 
and ft the target dilution factor, i.e. the fraction of polarized nucleons in the 
target material. 

Relations (15.451) and (15.461) simplify if the nucleon mass m is neglected. In- 
deed, in this case the kinematic factor 7^0, and the structure function g 2 
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decouples. This is generally a reasonable approximation in the small-x or 
large-Q 2 region, where the factor 7 2 is of the order of 10~ 2 — 10~ 3 . Then, one 
simply has 

ft^ ^iTflg?)) *- with Ai ^- (5 - 48) 

It should be finally noted that some of the experimental measurements con- 
cern the asymmetry A\ (or Am), while others concern values of the combina- 
tion of An and Aj_ which corresponds to g\jF\. The two quantities coincide 
for m = (aside from the depolarization factor in the case of Am), but they 
do not when mass corrections are included. Then, in the case A\ is measured, 
the related theoretical asymmetry is given by \g\ — (4m 2 x 2 /Q 2 )go\l ' F\ (or by 
the similar relation obtained by Eq. (I5.45P if A\\ is measured), in the second 
case simply by g\jF\. 



5.5 Target mass corrections 

A large part of experimental data in polarized deep-inelastic scattering are 
taken at relatively low values of Q 2 . Specifically, Q 2 is usually around few 
GeV 2 for data points in the small-x region (data at Q 2 < lGeV 2 are also 
available, but they are usually not included in perturbative analyses). In this 
kinematical region, contributions suppressed by inverse powers of Q 2 arise 
when taking into account the finite value of the nucleon mass m, and could 
play a relevant role. 

Another source of power-suppressed terms originate from the operator prod- 
uct expansion of the hadronic tensor W^ u in Eq. (15.11) . i.e. from matrix 
elements of operators of non-leading twist (see e.g. [5]). These contributions 
are usually referred to as dynamical higher twist, and their effect can not 
be calculated in perturbation theory, so that it is difficult to assess their 
relevance in any phenomenological analysis of polarized pdfs. 

When considering target mass corrections, in order to obtain experimental 
information on g\ from the measured asymmetry A\\, the complete relation 
Eq. (15.45J) must be used (or the first of Eq. (15.421) if the data concern Ai), 
that is, all the kinematical factors must be retained. Moreover, also the struc- 
ture function g 2 now contributes when comparing the theoretically predicted 
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observable g\ to the data. However, experimental data on the structure func- 
tion g 2 are restricted to a limited range in the (a;, Q 2 )-plane and are affected 
by large uncertainties (see e.g. [21]). Therefore, a model for g 2 is usually 
invoked in phenomenological analyses. 

To take target mass corrections into account, following Ref. [22], the mo- 
ments of the polarized structure functions g\ and g 2 can be expressed in 
terms of the matrix elements of the twist-2 and 3 operators appearing in the 
light-cone expansion of the forward scattering amplitude. Specifically, at the 
first order in m 2 /Q 2 , the n-th moments of g\ and g 2 read [22] 



m 2 n(n + 1) 



gi(Q) = an + 7^ / , ' {na n+2 + Ad n+2 ) + O [ — ) (5.49) 



m 



Q 2 (n + 2) 2V "^ "^ \Q 4 

„,^,, n — 1,, . m 2 n(n — 1) r , . . , _ f m A 

Here, a n and d n are given by matrix elements of the twist-2 and 3 polarized 
operators 

< p,s|z n_1 [i)m a D^ . . .D fl "- l X i tp] s \p,s >= -manM^ 1 -^ 1 (5.50) 

< p, s\i n - 1 [i) l5l x D a D^ . . . D^- 2 \i%l)] s , \p, s >= mdnM^ 1 -^- 2 

where m, p and s are mass, four-momentum and spin of the target nucleon, 
the symbol [■ ■ -]s means complete symmetrization in the indices a, jjl\ . . . /i n _i, 
whereas [. . .]s' denotes antisymmetrization with respect to A and a and sym- 
metrization over the other indices. The symbol M^ 1 -^- 1 denotes the most 
general rank-n tensor which can be formed with one spin four-vector s and 
n — 1 momentum four-vectors p; similarly ji^^ 1 — ** n - 2 j s antisymmetric in A 
and a and symmetric in all other indices. 

If m = is taken in Eqs. (I5.49J) the n-th moments of the structure functions 
(aside from QCD perturbative corrections) reduce to the simple form [23] 



m 



9io(Q 2 ) = a n + ( — I (5.51) 

One can then use Eq.s (15.511) to eliminate the matrix elements a n and d n 
from Eqs. (15.491) in favour of the n-th moments of the structure functions at 
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zero nucleoli mass, i.e. g™ and g^o- Then, to order 0(m 2 /Q 2 ) one has [22] 



9i(Q 2 ) = 9io(Q 2 ) + 

92(Q 2 ) = qUQ 



m 2 n(n + 1) 
Q 2 {n + 2) 2 
m 2 n(n — 1) 



(n + 4K +2 (Q 



4(n + 2) 
(n+1) 



^^^^(Q 2 )-^^ 



^20 +2 (Q 2 

(5.52) 



If 4 = 

lation 



Q 2 (n + 2) 2 

0, the structure functions obey the so-called Wandzura-Wilczek re- 

1 



gUQ z 



n 



n 



gUQ 2 )- 



(5.53) 



Actually, twist-3 contributions are not power-suppressed with respect to the 
leading twist; however, if d n -C a n is assumed, one can then use Eq. ( 15.531) to 
express g™ and g% in Eqs. f 1 5 . 5 2 j) in terms of the zero-mass structure function 
g™ alone, namely 



g n 2t\Q 2 ) 



n+l 



n 



girw 



(5.54) 



with g 1Q given by Eq. (15.231) . Thus, when including mass corrections in the 
analysis, the theoretically predicted structure function g\ (or its n-th mo- 
ment) as given by Eq. (15.231) is corrected by the first of Eqs. (I5.52p . together 
with the hypothesis (I5.54[) on the structure function #2 at zero mass. 



Chapter 6 

Global analysis of inclusive DIS 
data 



The problem of extracting relevant physical quantities from the existing in- 
clusive DIS data is considered here in detail. The analysis is mainly aimed 
at assessing the relevance of inclusive measurements to the determination of 
the spin-dependent gluon density, its first moment and the x-shape. 

To begin with, some details on the analysis method are given in Sec. 6.1, 
and the boundary conditions for parton distribution functions are discussed 
in Sec. 6.2. Sec. 6.3 is devoted to the study of the impact of the most 
recent inclusive DIS data on the extraction of polarized pdfs. In particular, 
the effect of very low-energy data is studied in Sec. 6.4, in connection with 
higher-twist corrections. 

The results on the gluon distribution function are given in Sec. 6.5., and their 
dependence upon input densities extensively investigated. Finally, the main 
consequences of the above analysis on the structure function gi and moments 
of quark and gluon polarizations are summarized in Sec. 6.6. 



6.1 Analysis method 

The nucleon structure function g\ is constructed at the measured points in 
the (x, Q 2 ) plane by convoluting the hard coefficient functions with the so- 
lutions of the NLO evolution equations. A parametrization for the polarized 
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pdfs in the singlet, nonsinglet and gluon sectors defined in Sec. 5, namely 
AS, A^ns an d Ag, as a function of a number of parameters is thus assigned 
at a given initial scale Qq. The whole set of parameters is then optimized by 
minimizing Q| the total x 2 5 namely 

i V £<?i data 0, g 2 ) ; ' l ' 

where the sum runs over the available data points. Experimental data for g\ 
are derived from the measured longitudinal spin asymmetry as discussed in 
Sec. 5.4, and are usually given along curves Q 2 {x) in a restricted region of the 
final-state invariant mass for each experimental set of data. The uncertainty 
5gf ata includes both statistical and systematic errors added in quadrature, 
and gi(x, Q 2 ) denotes the theoretically predicted observable. 
A full next-to-leading order analysis of the available experimental data on 
the polarized structure function g± of the nucleon has been performed in the 
Adler-Bardeen factorization scheme, as discussed in Sec. 5.3. 
The AP evolution is performed in the Mellin space of the n-th moments, 
and, by the inverse transformation, the structure function g\ in terms of 
the Bjorken-x is finally recovered for the fitting procedure, i.e. the x 2 - 
minimization. 

In order to extract g\ data from the experimental asymmetries, as discussed 
in Sec. 5.4, one needs a parametrization for both the ratio R = o^ja^ of the 
photoabsorption cross sections, and the unpolarized structure function F 2 . 
Because of the slightly better coverage of the small- x region, due to the 
COMPASS data (see Sec. 6.3 below), the parametrization R1990 [47] used 
in previous analyses [12], [39] for the ratio R = ul/ctt has been replaced by 
the most recent one R1998 [48]. The latter parametrization indeed extends 
the kinematic range to lower and higher values of x (0.005 < x < 0.86), by 
adding to the fit a considerable number of data [49] on R(x, Q 2 ). 
Furthermore, the commonly used parametrization of the unpolarized struc- 
ture function F 2 , which also enters into the experimental determination of 
<7i, given by NMC [50] has been replaced by the unbiased one [51], obtained 
by training sets of neural networks on the experimental data (see [5*2"]). 
However, no sizable effects are observed as a result of both these implemen- 



1 The x 2 minimization procedure is performed by the CERN program library MI- 
NUIT [25]. 
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tations, that is, all the parameters of pdfs turned out to be always largely 
consistent within the errors, and the quality of the fit essentially unchanged. 

Finally, it should be noted that target mass corrections have been included 
as a rule up to the first order in m 2 /Q 2 , as explained in Sec. 5.5, throughout 
this analysis. Indeed, when considering the finite size of the nucleon mass 
m, also the structure function g 2 enters into the theoretical determination of 
g 1 . Then, due to the lack of precise experimental data on g 2 , the Wandzura- 
Wilczek hypotesis has been invoked in the calculations (see Sec. 5.5). 



6.2 Parametrization of pdfs 

The functional form of spin- dependent parton distribution functions at the 
initial scale Ql is fixed according to the conventional parametrization 

Af(x,Ql) = N fVf x a f(l-xff(l+ V x s f) (type-A) (6.2) 

where Af denotes AS, Aq^s, and Ag. The overall factor Nf is defined by 
the normalization condition 

/ dxAf(x,Q 2 ) =r)f, (6.3) 

Jo 

i.e. such that the parameters r]f, fitted to the data, directly reflect the value 
of the first moment of pdfs at the initial scale, which is taken as a rule at 
Ql = 1 GeV 2 . 

The nonsinglet quark density at Ql = 1 GeV 2 is defined by the linear com- 
bination of the triplet and octect distributions 

Aq NS (x, Ql) = ±-Aq 3 (x, Ql) + -Aq 8 (x, Ql) , (6.4) 

where the plus (minus) sign stays for a proton (neutron) target. The distri- 
butions Ag 3 and Ag 8 are assumed to have the same x-dependence, while the 
parameter 77g, corresponding to the first moment of Agg, is fixed to the value 
?7s = 0.588 ± 0.025 by using SU(3) symmetry and octet baryon decay con- 
stants j8]. The parameter 773 is instead fitted to the data. At higher scales, 
new contributions are generated dynamically as heavy quark thresholds are 
crossed. 
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Since the large-x behavior of the polarized gluon distribution can not be de- 
termined by inclusive DIS data, i.e. by scaling violations, the value of /3 g 
needs to be fixed a priori. Typical values of j3 g examined in the following 
analysis are restricted to the range 4 < j3 g < 10, and a correlation with the 
gluon first moment at the initial scale r\ g is observed, in that increasing the 
(fixed) value of (3 g leads to a slight decrease of r\ g . Then, both parameters 
have been followed out as new data are included and the whole fit improved. 
Moreover, all the exponents 5f in the last factor of Eq. (16. 2p have been fixed, 
as in Ref . [12] , respectively to the values g 

<5 S = S g = 1 and 5 NS = 0.75 . (6.5) 

A different functional form has been also adopted for the initial parametriza- 
tions, to evaluate the impact of the input densities on the final results. Specif- 
ically, while keeping the initial scale at Q% — 1 GeV 2 , the less singular input 

AS = Nxr]xx as (In 1/xf* (type-B) (6.6) 

A/ = NfTjf [(In l/x) a f + 7/ x s * (In 1/xf'] A/ = Ag NS , A^ 

is used for an alternative fit, with the same normalization condition (16.31) . 
Similarly to the standard parametrization, the exponents 5f are fixed to the 
values S g — 1 and <5ns = 0.75. At variance with the initial densities (I6.2p . in 
Eq. (16.61) the rise at small x is at most logarithmic, and then softer than any 
power. Since In 1/x ~ (1 — x) as x — > 1, the large- x behavior in Eq. (16.61) is 
similar to Eq. (16.21) for the singlet distribution, whereas both terms contribute 
in the gluon and nonsinglet cases. 

The total number of parameters determined from the fit amounts, in both 
cases, to at most eleven, and no limits are imposed on their variation. Also 
the sign of all parameters are left free, including the first moment at the 
initial scale, although the data always choose rjf to be positive. 
Finally note that, in order to avoid any bias produced by further constraints, 
the input densities Eqs. (16.21) and (16. 6p are not explicitly forced by the positiv- 
ity condition \Af(x, Q 2 )\ < f(x,Q 2 ), at variance with many recent analyses 



2 The choice for Sns in particular is discussed in [26] , and is made in order to avoid 
asymptotic behaviors in the nonsinglet channel more singular than x~ 05 as x — ► 0, 
corresponding to the maximum saturation of the partonic constraint \Aq-^s(x,Q 2 )\ < 
1ns( x iQ 2 )- Indeed, a correlation between the parameters (5ns and cxns has been pointed 
out in [26j . such that cuns approaches zero as <5ns decreases from unity. 
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(see e.g. [271 [28]) . where, on the other hand, the initial parametrizations 
for the polarized gluon, valence and sea quark distributions are taken pro- 
portional to the related unpolarized quantities, with proportionality factors 
given by powers of x. However, the positivity bound has been checked a 
posteriori, and, in the case of fits of type-A, turns out to be always largely 
fulfilled (see also [221). 



6.3 Update 

The main purpose of the following analysis is to investigate the impact of 
the recent experimental data on the polarized parton distributions, with a 
special focus on the gluon polarization. The whole analysis is performed at 
NLO in the Adler-Bardeen factorization scheme, and the normalization of 
the strong coupling a s (M§) = 0.118 is assumed throughout. 
The starting point is a global fit of type-A (Eq. (16.21) ). of inclusive DIS data, 
on the basis of 176 experimental points used in the previous analyses [121 ES]) 
as summarized in Tab. 1. 



Data set (target) 



x-range 



Q -range (GeV 



No. data 



Ref. 



SMC (p) 
E143 (p,d) 
E155 (p) 
HERMES (p) 
SMC (d) 
E155 (d) 
E142 (n) 
HERMES (n) 
E154 (n) 



0.005 
0.031 
0.015 
0.023 
0.005 
0.015 
0.035 
0.033 
0.014 



- 0.479 

- 0.749 
-0.75 
-0.66 

- 0.480 
-0.75 

- 0.466 

- 0.464 

- 0.564 



58 
9.52 
34.72 
- 7.36 
54.8 
34.79 
1.1-5.5 
1.22 - 5.25 
1.2 - 15 



1.3 

1.27 

1.22 

1.01 

1.3 

1.22 



12 

56 

21 

20 

12 

21 

8 

9 

11 



m 

m 

[32] 
[33] 

EDI 



[35J 
[36] 



Tab. 1 Initial set of inclusive DIS data (176 experimental points). 



Note that the effects of the E155 proton data [32] alone have been analyzed 
in recent works [38], and the conclusion was reached that their accuracy 
significantly contributes to reducing the statistical errors of pdfs. 
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As can be seen by the scatter plot in Fig. 16.11 the covered region of the 
(x, Q 2 )-plane is roughly restricted to 0.01 < x < 0.8 for the Bjorken x (with 
the exception of very few SMC points at lower x), and to 1 < Q 2 < 58 GeV 2 



for the momentum transfer, 
rule to all the data. 



Indeed the cut Q 2 > 1 GeV 2 is imposed as a 
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Figure 6.1: Distribution in the (x, Q 2 )-plane of the set of 176 DIS data up to 2000, as 
summarized in Tab. 1. 



The results of a fit of type-A on this set of data are shown in the second 

column of Tab. 3. By comparing such results with those of the previous 

analyses [121 EH], globally the same picture arises, and, in particular, the 

x-shapes of pdfs are similar to those of Refs. [121 EH], and likewise not very 

precisely determined. 

A sizable first moment of the polarized gluon density at the initial scale is 

found, and, both the gluon and the nonsinglet quark distributions rise as 

x —* 0, whereas a flat behavior is observed in the singlet sector. 

Note that the parameters 7 9 and 7s, which control the shape of the gluon 

and singlet distributions at intermediate x, can not be disentangled in this 
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analysis. As observed in [20], this is due to the mixing of Ag and AE in 
the AP evolution. Therefore 7 5 = 7s is taken, thus leaving ten parameters 
to be determined by the x 2 minimization. Finally, the errors given here are 
statistical errors from the fit. 

A bulk of more precise data has become recently available. These are in 
particular 15 points at rather high energy from COMPASS U, and a number 
of lower energy points from HERMES and JLAB as summarized in Tab. 2. 
Furthermore, a set of complementary data from SMC has been also included, 
that was not considered in the previous analyses [121 ESj • A large amount 
of very low-energy data from CLAS [IH H6] will be considered separately, in 
connection with higher-twist effects (see Sec. 6.3). 



Data set (target) 


x-range 


Q 2 -range (GeV 2 ) 


No. data 


Ref. 


COMPASS (d) 
HERMES (p,d) 
JLAB 

SMC (p,d) 


0.0046 - 0.566 
0.0264- 0.7311 
0.33 - 0.60 
0.0043 - 0.121 


1.10- 55.3 
1.12- 14.29 
2.71 - 4.83 
1.09- 23.1 


15 
74 
3 
16 


[44j 



Tab. 2. New set of inclusive DIS data up to 2006. 



As shown in Figs. 16.21 and 16.31 such data slightly improve the coverage of 
the (x, Q 2 )-plane, and in particular the low-rc region down to x > 0.0046 
(COMPASS). 

The results of a fit of type-A, including all these most recent data is displayed 
in the third column of Tab. 3. Due to the increased number of experimental 
points and their accuracy, an overall reduction of the statistical errors from 
the fit is observed. Aside from a somewhat flatter behavior at small-a; of the 
polarized singlet distribution, essentially the same shape of pdfs turns out 
in both the singlet and nonsinglet sectors, which are displayed at the initial 
scale in Fig. 16.41 In particular, the new small-x data points confirm the rise 
of the nonsinglet distribution Ag^s, already determined in the previous anal- 
yses [T2"| |2"6] . The singlet and nonsinglet first moments are fairly unaffected 
by the new experimental data. 



3 Note that the latest COMPASS data [40] at very low x and Q 2 have been excluded 
altogether in this analysis, all lying below the cut Q 2 > 1 GeV 2 imposed to the data set. 
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Figure 6.2: Distribution in the (x, Q 2 )-plane of the new set of DIS data up to 2006, as 
summarized in Tab. 2. 
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Figure 6.3: Distribution in the (x, Q 2 )-planc of the complete set of data. 
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Parameters (Q 2 , = 1 GeV 2 ) 


176 data 


284 data 


a s 

Pz 

7s 

Vg 
a g 

Pg 

lg 

V3 

&NS 

Pns 
Ins 


0.347 ± 0.025 
0.732 ± 0.504 
4.307±0.626 
10.044 ±16.356 
0.856±0.666 
-0.772 ± 0.443 
4. (fixed) 
10.044 ±16.356 
1.119 ±0.035 
-0.364 ± 0.307 
2.982± 0.291 
9.749± 15.426 


0.344 ± 0.011 
1.421 ± 0.167 
3.212 ± 0.918 
-0.673 ± 0.805 
0.402 ± 0.062 
-0.300 ± 0.194 
10. (fixed) 
-0.673 ± 0.805 
1.095 ± 0.024 
-0.276 ± 0.165 
3.186 ± 0.157 
7.356 ± 6.271 


rVd.o.f 


0.936 


0.879 


a (10GeV 2 ) 
r?(10GeV 2 ) 


0.142±0.023 
0.113±0.002 


0.234 ± 0.010 
0.121 ± 0.001 



Tab. 3. Fits of typc-A, respectively to the set of 176 data points of Tab. 1 (second column) 
and to the full set including data of Tab. 2 (third column). 



The most relevant feature here is the sizable reduction of the gluon first 
moment and of its statistical error. Actually, such a low value of r\ g is also 
correlated to the large fixed value of j3 g , that seems, however, the value 
preferred by the data. Indeed, the same fit has been repeated by exploring 
lower (fixed) values of j3 g , and the minimum \ 2 is found, for this set of data, 
with (3 g = 10. The correlation between the parameter that controls the large- 
x shape of Ag and the first moment turns out to be a distinctive feature of 
the type-A parametrization, and is not observed in the case of the type-B fit 
(see Sec. 6.5.2). A deeper discussion of the impact of each set of experimental 
points on r\ g is given in Sec. 6.5.1. 

The resulting polarized gluon density (type-A) is displayed in Fig. 16.51 with 
increasing scale Q 2 (red lines), and compared with the analogous result of 
the reference fit to 176 data points (black lines). 

Finally, the value of the singlet axial charge ao(Q 2 ) is considerably increased, 
and is given in Tab. 3 at Q 2 = 10 GeV 2 for both fits. Also shown in Tab. 3, 
as an example, is the first moment of the polarized structure function gx 
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for a proton target at the same scale and computed over the full arrange. 
Complete results for gi are summarized in Sec. 6.6.1. 



6.4 Low-energy data 

A large amount of low-energy and very precise data points from the CLAS col- 
laboration [351 ES], both on proton and deuteron targets, has been left aside in 
the previous analysis in order to disentangle the effects of high and low energy 
data on polarized pdfs. Indeed, CLAS data are restricted to a limited re- 
gion of the final-state hadronic invariant mass, roughly 1 GeV< W < 3 GeV, 
where non-perturbative effects could be relevant. This kinematical region in- 
volves large Bjorken-x at moderate values of the squared momentum transfer 
Q 2 , and is characterized by the presence of nucleon resonances which con- 
tribute to higher-twist effects in the structure functions. 
In order to perform a consistent perturbative analysis of world data, a quan- 
titative criterion for selecting experimental points from CLAS has been stud- 
ied. To this end, a lower bound on the invariant mass W has been prelimi- 
narly estimated on the basis of the relevance of higher-twist contributions to 
the moments of the polarized structure function g-±. 

As a next step, the precise location of the kinematic cut to be imposed on 
the whole set of data has been determined by explicitly evaluating residual 
higher-twist effects. 



6.4.1 Selection of CLAS data 

In order to fix a lower cut on the (x, Q 2 ) plane to the full set of CLAS data, 
a recent analysis [53] of the higher moments of the proton structure function 
g\ has been exploited. All the available data, both in the DIS and resonance 
region [35], down to very low Q 2 , are used in [53] to estimate the higher-twist 
contributions (twist-4 and 6) to the n-th moments of g\ 

M n (Q 2 ) = 5 Vn (Q 2 ) + HT n (Q 2 ) . (6.7) 

with n = 3,5,7. Here the leading and higher-twist terms, Sr] n (Q 2 ) and 
HT n {Q 2 ) respectively (see Eqs. (32) and (45) in Ref. [53]), have been deter- 
mined by fitting a number of parameters to the data. As expected, the total 
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Figure 6.4: Singlet and nonsinglet distributions (proton target) at Q 2 
fit of type-A on the whole set of available data. 
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Figure 6.5: Best-fit forms (type-A) of xk.g(x) in both fits, respectively to the set of 176 
data points (black curves) and to the full set including data of Tab. 2 (red curves), with 
incrcsing Q 2 . 
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higher-twist term turned out [53] to be sizable mainly for Q 2 ~ few GeV 2 and 
is still non-negligible even at C} 2 — 10 GeV 2 for the higher moments. More- 
over, it has been pointed out in [53] that the total higher-twist contribution 
is significantly larger in the polarized case than in the unpolarized one (see 
also [53]). 

On the basis of the above analysis, the ratio HT n (Q 2 ) / 5r] n (Q 2 ) can be easily 
evaluated for each nona wide range of Q 2 . This fact has been used here to 
select a lower bound on the momentum transfer, Qq, for each n, by requiring 
HT n (Q 2 )/5 Vn (Q 2 ) < 10% for Q 2 >Q 2 . 

Then, in order to define a proper cut to CLAS data in the (x, C} 2 )-plane, to 
each n a corresponding value of x must be also assigned. To this end, the 
integrand x n ~ l gi(x, Qq) of the n-th moment of the proton structure function 
g\ has been evaluated at the quoted Q\, and the position of the peak, x , 
has been taken as a characteristic value of x. 



C\J 
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CLAS data 



W=2.5 GeV 
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Figure 6.6: Set of CLAS data cut at W > 2.5 GeV used in the present analysis. 



By combining these values of xo and Qq for each n — 3, 5, 7, the corresponding 
invariant mass Wq has been computed according to Eq. (I5.3I) . It has been 
found W = 1.93, 2.48, 2.96 GeV for the moments n — 3,5, 7, respectively. 
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Then, taking an average over the three values of Wq, one roughly finds Wo ~ 
2.5 GeV as a lower bound on the final-state invariant mass. One may thus 
conclude that data points below this threshold can not be safely included 
in the analysis without taking into account higher-twist corrections, whose 
contribution amounts to at least 10% of the leading twist term in this region 
of the (x, Q 2 )-plane. 

As a first result, this leads to a preliminary selection on the full set of CLAS 
data. Indeed, only 148 experimental points lying above the cut Wo = 2.5 GeV 
are retained in what follows, and their distribution in the (x, Q 2 )-plane is 
displayed in Fig. 



6.4.2 Higher-twist effects 

As a result of the above analysis, a lower bound should be consistently 
imposed on the whole set of available data. Indeed, all the data below 
Wo = 2.5 GeV must be certainly discarded, and the precise position of the 
global cut has been investigated by studying the related effects of residual 
higher twist corrections on the results. 

Three possible cuts are compared in Fig. 16.71 with the whole set of data, 
including the selected CLAS data (cut at W > 2.5 GeV). As can be seen, a 
cut at W > 3 GeV would entail the exclusion of all the experimental points 
from CLAS (with the exception of only one point on proton target), and of 
a number of new HERMES data as well. Higher cuts would mean a sizeable 
loss of information, whereas imposing the lower bound at W — 2.5 GeV does 
not significantly reduce the set of data but higher twist effects could still play 
a relevant role in this region. 

This issue has been quantitatively evaluated by introducing a phenomenolog- 
ical term for the structure function gi, so as to take into account higher-twist 
contributions in the Bjorken-x space, namely 

9l (x,Q 2 ) - gi ( x ,Q 2 )[l + ^\ . (6.8) 

The multiplicative factor in Eq. (16.81) depends on one dimensional parameter 
c to be fitted to the data, in order to study the role of lower energy points in 
the determination of higher-twist corrections. A broad range of possible cuts 
Wo has been investigated; specifically, the same fit of type-A, with one more 
fitting parameter c, has been performed on the whole set of data (including 
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Figure 6.7: Distribution in the (x, Q 2 )-plane of the whole set of data including selected 
CLAS data, compared with lines of fixed final-state invariant mass W = 2.5, 3 and 3.5 GcV. 



the selected subset of CLAS data), progressively cut at higher values W of 
the invariant mass. 

The values of the HT parameter c with their errors, as determined by these 
fits, are displayed in Fig. I6.8I as a function of the cut W imposed. As a 
result, the higher- twist contributions turn out to be gradually less relevant 
as the full set of data is restricted to higher values of W. Indeed, as can be 
seen by Fig. I6.81 the largeness of the error bar increases with the value of the 
cut Wo, making the higher- twist parameter more consistent with zero. 
A glance to Figs. I6.7I and I6.8I reveals that the choice W = 3 GeV, where 
the higher- twist parameter c differs from zero by no more than 1.5 standard 
deviations, seems a reasonable one, in that it justifies a purely perturbative 
analysis while reducing the loss of experimental information. 

In Tab. 4 are shown, as an example, the results of two fits of type-A with 
11 fitted parameters and different cuts on W, i.e., no cut at all (except for 
CLAS data cut at W > 2.5 GeV) and W > 3.5 GeV over all the available 
points (which excludes CLAS data altogether). 
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w 

Figure 6.8: Values with errors of the parameter c fitted to the set of data gradually cut 
at higher W. 



Parameters (Qq = 1 GeV 2 ) 


no cut 


W > 3.5 GeV 


Wz 


0.401 ± 0.014 


0.397 ± 0.021 


Q.T. 


1.357 ± 0.273 


1.347 ± 0.304 


& 


2.655 ± 1.233 


2.435 ± 1.198 


7s 


-1.015 ± 0.448 


-1.191 ± 0.320 


Vg 


0.444 ± 0.066 


0.494 ± 0.075 


la 


-0.388 ± 0.184 
8. (fixed) 
-1.015 ± 0.448 


-0.451 ± 0.179 
6. (fixed) 
-1.191 ± 0.320 


V3 


1.129 ± 0.027 


1.137 ± 0.041 


&NS 


-0.236 ± 0.158 


-0.327 ± 0.218 


(3ns 


3.119 ± 0.163 


2.903 ± 0.300 


Ins 


7.114 ± 5.630 


8.979 ± 10.051 


c 


-0.561 ± 0.210 


-0.315 ± 0.650 


rVd.o.f 


0.937 


1.056 



Tab. 4. Fits of type-A, respectively to the full set of data with no cut except for CLAS 
(W > 2.5 GeV) (second column), and with a global cut at W > 3.5 GeV (third column). 



109 



Comparing the second column of Tab. 4 with the third one of Tab. 3 it can 
be seen that CLAS data do not sizably affect the shape of quark and gluon 
densities, whereas the main effects of the lower energy points, all cut in the 
third column of Tab. 4, is essentially the improvement in the determination 
of higher-twist contributions to g\. Similar conclusions were also reached in 
Ref. [55] where the impact of a much extended subset of CLAS data [IB] (i.e. 
W > 2 GeV) on polarized pdfs has been studied. 



6.5 Gluon polarization from g\ data 

As discussed in the previous section, in order to perform a consistent per- 
turbative treatment of world data, avoiding systematic errors induced by 
higher-twist effects, the final set of experimental points retained is the one 
selected by the cut on the invariant mass at W > 3 GeV, which amounts to 
238 data points and covers the region above the solid curve in Fig. 16.71 In- 
deed, in this kinematical region the higher-twist contributions can be safely 
neglected. 

The corresponding experimental asymmetries on proton, neutron and deute- 
ron targets included in the analysis are given as a function of x in Figs. 16.91 
16.101 and 16.111 respectively. 

As already noted, all the CLAS data (except for one point) have been ex- 
cluded, together with the three points from JLAB and also a number of lower 
energy HERMES data. Such a cut has obviously no effect on the higher en- 
ergy data from COMPASS and SMC. 

The information that can be extracted specifically on the polarized gluon 
density by a global fit of the present inclusive DIS data are discussed in the 
following sections. 



6.5.1 Best-fit results and impact of the data on the 
gluon first moment 

The results of a type-A fit to this set of experimental points is shown in 
Tab. 5. Here, higher-twist corrections are switched off (c = fixed), whereas 
target mass corrections are included as a rule. As a cross check, it has been 
verified that no sizable differences can be traced with respect to the results 
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Figure 6.9: Asymmetry data on proton target summarized in Tabs. 1 and 2, cut at 
W>3GeV. 
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Figure 6.10: Asymmetry data on neutron target summarized in Tabs. 1 and 2, cut at 
W > 3 GeV. 
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Figure 6.11: Asymmetry data on dcuteron target summarized in Tabs. 1 and 2, cut at 
W > 3 GeV. 



of a similar fit performed on the same set of data where the higher-twist 
parameter c has been fitted (c = —0.76 ± 0.48). 

Comparing the results of Tab. 5 with those of the third column of Tab. 3 
(including the whole set of data except for CLAS) makes clear that the main 
effect of the cut on the quark and gluon polarizations amounts to an overall 
increase of the statistical errors due to the loss of a number of experimental 
points, mainly from HERMES. On the other hand, the systematic error due 
to higher-twist effects is reduced by excluding low-energy data. 
The x-shapes of the quark singlet and nonsinglet distributions are very similar 
to those displayed in Figs. 16.41 and slightly differ only for the normalization. 
Specifically, the nonsinglet quark distribution rises at small-x, whereas the 
Adler-Bardeen singlet quark turns out to be flat at the initial scale. The 
best-fit form of quark distributions is displayed in Fig. 16.121 (black lines). 
The first moment of the polarized gluon density (and its statistical error) at 
the initial scale is slightly increased, but is neverthless much smaller than 
the previous determinations (see Tab. 3, second column), in agreement with 
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many current analyses (see Sec. 6.5). 

The parameter 7 9 in Tab. 5 has been disentangled from 7s, but clearly the 
present data can not give precise constrains on its value. Again, the behavior 
of Ag as x — > 1 can not be determined by inclusive DIS measurements and 
the parameter (3 g has to be fixed. The range 4 < (3 g < 10 has been explored 
and the minimum x 2 is found for f3„ = 6. 



Parameters (Ql = 1 GeV 2 ) 


type-A (W > 3 GeV, 


238 data) 


m 


0.391 ± 0.021 




as 


1.154 ± 0.334 




fc 


1.918±0.972 




7s 


-1.232 ± 0.192 




Vg 


0.555 ±0.146 




a g 


-0.484 ±0.293 




Pg 


6. (fixed) 




7 5 


3.829 ±10.443 




V3 


1.127 ± 0.034 




®NS 


-0.321 ± 0.271 




Pns 


2.970 ± 0.336 




Ins 


9.475 ± 13.376 




X 2 /d.o.f 


0.932 


a (10GeV 2 ) 


0.246 ±0.025 




r?(10GeV 2 ) 


0.125 ± 0.003 





Tab. 5. Best-fit result of type-A on the selected data set, and c = fixed. 



data set 



Vg (Q 



I = 1 GeV 2 



f3 g (fixed) 



rVd-o.f. 



176 set 

+ SMC2 

± COMPASS 

± HERMES-06 

+ JLAB 

238 set (W > 3 GeV) 



0.86 ± 0.67 
0.69 ± 0.52 
0.53 ± 0.14 
0.43 ± 0.07 
0.40 ± 0.06 
0.56 ± 0.15 



4. 
6. 



10. 
6. 



0.94 
0.96 
0.93 

0.88 
0.88 
0.93 



Tab. 6. Impact of each set of experimental data on the value of the first moment rj g of the 
gluon distribution. Also shown are the corresponding value of j3 g and the related % 2 . 
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The x-shape of the best-fit result of type-A (black curves) for the polarized 
gluon distribution is displayed in Fig. 16. 131 with increasing momentum trans- 
fer. 

Finally, the impact of each set of data on the first moment of the gluon polar- 
ization is summarized in Tab. 6, together with the corresponding best value 
of j3 g and the related x 2 '• As can be seen, all the measurements contribute to 
decreasing r) g , the much relevant effect being that from the high-energy data 
by COMPASS, lying well above the cut imposed on the (x, Q 2 )-plane. 

6.5.2 Dependence of Ag on input densities 

The assumed functional form of input densities sizably affects the determi- 
nation of relevant quantities, such as the x-shape of the gluon polarization. 
A way of evaluating this effect is to perform the whole analysis by chang- 
ing the initial parametrization of pdfs. This has been done here by using 
the functional form given by Eq. (16. 6p (type-B) for an alternative fit on the 
selected set of data cut at W > 3 GeV used in the type-A fit (Tab. 5). 



Parameters (Ql = 1 GeV 2 ) 


type-B (W > GeV, 


238 data) 


rm 


0.417 ± 0.023 




as 


2.496 ± 0.425 




Pa 


3.512±0.467 




Vg 


0.801±0.190 




a g 


1.050 ± 0.486 




Pg 


3.113±0.868 




lg 


-1.447±0.647 




V3 


1.219 ±0.038 




&NS 


1.733 ± 0.121 




Pns 


4.655± 0.138 




Ins 


-0.269± 0.060 




rVd.cf 


0.926 


a (10GeV 2 ) 


0.218±0.027 




r?(10 GeV 2 ) 


0.129±0.004 





Tab. 7. Best-fit result of type-B on the selected data set, and c = fixed. 



114 



The analysis is performed in the Adler-Bardeen scheme, and the results 
for the pdf parameters, the singlet axial charge ao(Q 2 ) and the first mo- 
ment Ti(Q 2 ) of the proton structure function at Q 2 = 10 GeV 2 are given 
in Tab. 7. These results may be compared with those of the related fit of 
type-A (Tab. 5). 

In the quark sector first moments turns out to be compatible within the er- 
rors with those from the type-A fit. The resulting x-shapes of singlet and 
nonsinglet distributions are also very similar to the best-fit forms of type-A, 
and are displayed in Fig. 16.121 at the initial scale (red lines) . 
A stronger sensitivity to the initial functional form is observed in the gluon 
sector. The central value of the gluon first moment turns out to be larger 
than that of type-A, even though still compatible within the errors. 
On the other hand, the x-shape of Ag is considerably changed with respect 
to type-A, and both are displayed in Fig. 16.131 with increasing momentum 
transfer Q 2 . Note that, at variance with the fits of type-A, the behavior of 
Ag as x — > 1 has not been fixed here a priori. 

A quantitative estimate of the systematic uncertainty produced by the choice 
of input densities may be performed by analyzing distinctive features of the 
x-shape in both type-A and B gluon distributions, such as the position of 
the peak and the corresponding value of the function. This study has been 
performed in Sec 7, in connection with pseudo-data analyses for charm photo- 
production at the COMPASS experiment that aim at measuring Ag directly 
(see Sec. 7). 

Finally, being the quality of the fit fairly the same in both cases, no one of 
the assumed functional forms is neatly preferred by the data, that is both 
choices of input densities give equally good fits of inclusive world data, al- 
though the gluon x-shapes are significantly different. This result makes clear 
the large bias produced by the initial parametrization on the spin-dependent 
gluon distribution. This fact, sometimes hidden within most of the avail- 
able phenomenological studies, mainly performed by imposing standard in- 
put densities (and constrained by the positivity condition), has been here 
clearly stated and taken into account in the estimation of the theoretical 
errors (see Sec. 6.6.2 below, and also Sec. 7.3). One may thus conclude 
that, even though the better coverage of the (x, Q 2 )-plane, mainly due to 
the precise COMPASS data, allows a more reliable determination of the first 
moment of the gluon polarization by scaling violations, its x-shape is still 
largely unconstrained by the present inclusive DIS data. 
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Figure 6.12: Best- fit results of type- A and B (black and red lines respectively) for the 
polarized quark singlet and nonsinglet distributions (proton target) at Q 2 = 1 GeV 2 , over 
the selected set of available data (W > 3 GeV). 
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Figure 6.13: Best- fit results for the polarized gluon distribution with increasing Q 2 , over 
the selected set of available data (W > 3 GeV). The black curves denote the fit of type- A, 
and the red ones the fit of type-B. 
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6.6 Phenomenological implications 

The main results of the above analysis as far as phisically relevant quantities 
are concerned are summarized in the following sections. Specifically, the 
best-fit forms of structure functions are presented, and the final estimate for 
the singlet axial charge and the polarization of the gluon and the quark flavor 
singlet combination in the nucleon are discussed. 



6.6.1 Best-fit results for g\ 

The best-fit structure functions g\ (type-A) for proton, neutron and deuteron 
targets, are displayed in Figs. 16.141 16.151 and 16.161 respectively, in the range 
of momentum transfer covered by the data (i.e. Q 2 = 1, 10 GeV 2 ), and com- 
pared with the experimental points with their errors used in the analysis. 
Note that the sharp rise at large-x observed in all the three curves at the 
initial scale (Q 2 = 1 GeV 2 ) is due to the inclusion of target mass correc- 
tions, and it disappears if m 2 = is taken exactly for the nucleon mass. 
Such an effect, however, is not physical since the approximation adopted for 
the calculation of target mass corrections becomes unreliable at very large 
x (see [22]). Indeed, because of the m 2 /Q 2 suppression factor, this unphys- 
ical peak is washed out at higher Q 2 , which is actually the experimentally 
interesting region for large-x. 

Figs. I6.17[ 16.181 and 16.191 compare the best-fit forms of g\ for each target, 
corresponding to both fits, type-A and B respectively, at Q 2 = 10 GeV 2 . 
Sizable differences can be traced in the small-x region in the two fits. This 
fact is due to the way the functional forms of pdfs extrapolate the behavior 
observed in the last measured points, emphasizing the impact of the initial 
parametrization on the small- x behavior of the structure functions. In all 
cases g\ is driven negative at small- x by the rise of the polarized gluon den- 
sity in this region. 

The values of the first moment of g± obtained by integration over both the 
full rr-range and only over the measured range, are given in Tab. 8, at typ- 
ical values of Q 2 for each target. Both fits, type-A and B respectively, are 
considered. The truncated moments turn out to be generally quite close to 
each other. 



117 



X 

Q- 



0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




-0.01 



1 1 1 1 — i — i — n - 1 — 

Q 2 =10.GeV 2 

Q 2 =1.GeV 2 

CLAS :-+--- 
SMC i x 

E143 t— X-- 

E155 ^-&- 

HERMES ------- 

HERMES06 :- --•© - 



- 1 1 — I — I I I I 



- 1 1 1 — i — i — r 




0.001 



0.01 



0.1 



Figure 6.14: Best-fit result (type- A) for the proton structure function compared with 
the data points used in the analyis (W > 3 GeV) . 
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Figure 6.15: Best-fit result (type- A) for the neutron structure function compared with 
the data points used in the analyis (W > 3 GeV) . 
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Figure 6.16: Best- fit result (type- A) for the deuteron structure function compared with 
the data points used in the analyis (W > 3GeV). 
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Figure 6.17: Best-fit results of type-A (solid) and type-B (dot-dashed) for the proton 
structure function. 
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Figure 6.18: Best-fit results of type-A (solid) and type-B (dot-dashed) for the neutron 
structure function. 
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Figure 6.19: Best-fit results of type-A (solid) and type-B (dot-dashed) for the deuteron 
structure function. 
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rf(10GeV") 


1^(10 GeV") 


r«(5GeV") 


type-A (full range) 
type-A (meas. range) 


0.125 ± 0.003 
0.129 


0.038 ± 0.003 
0.047 


-0.047 ± 0.004 
-0.029 


type-B (full range) 
type-B (meas. range) 


0.129 ± 0.004 
0.128 


0.035 ± 0.003 
0.044 


-0.058 ± 0.005 
-0.032 



Tab. 8. First moment of the structure function g\ for proton, deuteron and neutron 
targets. 



6.6.2 Polarized quark and gluon densities 

Moments of quark singlet and gluon polarizations as well as the value of the 
singlet axial charge may be evaluated by taking the average over the central 
values in the two fits, type-A and B respectively, and the final estimates thus 
read 



AE(1) = 0.40 ± 0.02 (exp) ± 0.06 (th) 
A#(l, 1 GeV 2 ) = 0.68 ± 0.12 (exp) ± 0.29 (th) 
a (10GeV 2 ) =0.23 ± 0.02 (exp) 



(6.9) 



: °:^(th). 



The first error quoted in Eqs. (16. 9p is the one related to the statistical errors 
from the fit, the second error is the theoretical one. 

Indeed, as discussed in [El [26], the main sources of theoretical uncertainty in 
the determination of quark and gluon polarizations, and of the related phys- 
ical quantities, stem from the truncation of the perturbative expansions, as 
well as from the choice of the assumed initial parametrization of pdfs as em- 
phasized in Sec. 6.5.2. 

The error on the gluon first moment Ag(l, 1 GeV 2 ) due to the choice of input 
densities can be evaluated by comparing the results of the type-A and B fits, 
i.e. by taking the half- difference of the two central values as an estimate of 
the related uncertainty. It turns out roughly ~ ±0.12. Similarly for the sin- 
glet quark and the axial charge one has ~ ±0.05 and ~ ±0.01 respectively. 
The uncertainty induced by the truncation of the perturbative series for the 
coefficient functions and the evolution kernels is reflected by the dependence 
of the results on the renormalization scale fiR and on the factorization scale 
fip respectively, that are usually both fixed to the value of Q 2 of each data 
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point. A way of studying the effects of higher orders is by changing the value 
of hr and jIf around the chosen values. On the basis of the previous anal- 
yses [561 [T2l [26] . the theoretical error on Ag(l, 1 GeV 2 ) due to higher order 
corrections has been estimated to be roughly ~ ±0.26. Similarly for AE(1) 
and a (10GeV 2 ) one has ~ ±0.03 and ~ 1°;^ respectively. 
Another source of theoretical uncertainty is due to higher-twist corrections, 
and, as extensively discussed in Sec. 2.4, the related systematic error has 
been strongly reduced in the present analysis by discarding very low energy 
data, i.e. below the cut at the final-state invariant mass W = 3 GeV. How- 
ever, by exploting results of Sec. 2.4.2, an estimate of higher-twist effects may 
be attained by comparing the results of the same fit of type-A to the whole 
set of available data (Tab. 4, second column) with the one performed on the 
selected subset cut at W > 3 GeV (Tab. 5). The related error on the gluon 
first moment Ag(l,lGeV 2 ) amounts to ~ ±0.06, and it is thus negligibly 
small if compared to the theoretical errors from the higher orders and the 
initial parametrization. For AE(1) and ao(10GeV ) one finds ~ ±0.01 and 
~ ±0.02 respectively. 

The overall theoretical uncertainties in Eqs. (16.91) is then obtained by com- 
bining the errors quoted above. Other sources of theoretical error, such as 
the uncertainty on the value of a s (M§) and 77s, and the one related to the 
positions of heavy quark thresholds, all entail very small effects [56] and have 
been neglected in the final estimates. Also the impact of target mass correc- 
tions is found to be small (see [55J 1ST]). 

In order to compare these results with those of current analyses, it should 
be recalled that in the Adler-Bardeen factorization scheme the polarized sin- 
glet quark density AE(1) is scale independent, whereas in the MS scheme, 
adopted in most of the recent works [271 EH EH EH] , it coincides with the (non 
conserved) singlet axial charge a>o(Q 2 ) (see Eq. ( 15.321) ). The polarized gluon 
density Ag should be the same in the two schemes at NLO accuracy. A rea- 
sonable agreement between the above results and the most recent estimates 
is found. As an example, in Ref. [27] the quoted values read 



AE(1, 1 GeV 2 ) = 0.25 ± 0.10, Ag(l, 1 GeV 2 ) = 0.47 ± 1.08 , (6.10) 

that are very similar to their previous estimates |38j . 

Furthermore, in Refs. [HJ [55] two solutions with either Ag > or Ag < 

have been found; the first moment of Ag at the initial scale turns out to be 
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small in absolute value, i.e. \rj g \ < 0.3 [HI [55] for both solutions, but the 
shapes of the distributions are very different. It is finally worth noting that 
in the analysis of Ref. [H] the positive solution for Ag seems to be favored 

by the data. 
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Figure 6.20: Best- fit result (type- A) for xAg(x) compared at Q 2 
other paramctrizations. i.e. AAC [38], BB [58] and GRSV [59] , 



1 GeV 2 with the 



Then, even though the most recent global fits to world data fairly agree within 
the errors as far as the first moment of the polarized gluon density is con- 
cerned, they produce very different shapes of Ag, as shown in Fig. 16.201 Here, 
the best-fit result of type-A (Tab. 5) is compared at the scale Q 2 = 1 GeV 2 
with other available parametrizations [38 ] 158], 159] of the spin-dependent gluon 
distribution. One may thus concludes that to the present accuracy inclusive 
DIS data are not capable to discriminate among the existing different scenar- 
ios, and thus to pin down the x-shape of the gluon polarization. As discussed 
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extensively in the next section, some light may be shed on this issue by direct 
measurements of Ag, rather than by scaling violations, that is, by including 
in the analysis experimental data from exclusive reactions in which gluons 
contribute at leading order. 



Chapter 7 

Open-charm photoproduction 
at COMPASS 



The spin-dependent gluon distribution Ag in the nucleon seems to be hardly 
constrained by the inclusive DIS measurements. Actually, as discussed in 
Sec. 6, a global fit to world data including the most recent and precise ex- 
perimental results on the polarized structure function g±, allows a better 
determination of its first moment, that turns out to be much smaller than 
previous estimates (see e.g. [121 [60]). However, the polarized gluon density 
exhibits a sizable dependence on the assumed functional form (see Sec. 6.4.2), 
and, as a result, the x-shape is still largely unknown. 

At present, experimental constraints on the gluon polarization are expected 
to come from exclusive reactions, aimed at measuring Ag directly, rather 
than through scaling violations. Such measurements entail the observation 
of specific events that receive leading contributions from gluon initiated sub- 
processes, such as, for example, production of heavy-flavored hadrons and 
jets with large transverse momentum. In particular, as discussed below, the 
open-charm photoproduction represents a promising approach for extract- 
ing Ag(x). This channel is currently experimentally studied at COMPASS, 
and first results on the spin asymmetry are upcoming. Therefore, a com- 
plete analysis of polarized photoproduction of open-charm at the COMPASS 
kinematics has been performed in the present work, in order to assess the 
relevance of such measurements in constraining the x-shape of Ag(x). 

A brief overview of the current scenario on direct measurements of the gluon 
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polarization is sketched in Sec. 7.1. Sec. 7.2 is devoted to the theoretical 
study of the COMPASS experiment, and finally phenomenological results 
are discussed in Sec. 7.3. 



7.1 Direct determination of Ag 

Reactions dominated by gluon initiated processes at the partonic level can be 
investigated with both electromagnetic and strong probes, that is, in lepto- 
and photo-production reactions and hadro-production respectively. 



7.1.1 Experiments at polarized proton-proton colliders 

Inelastic proton-proton collisions with polarized beams at high energies pro- 
vide the possibility to investigate the gluon spin distribution by using strong 
interaction among partons from the colliding hadrons. 

New perspectives have been opened up by the RHIC-Spin (Relativistic Heavy 
Ion Collider) machine at Brookhaven National Laboratories, where a num- 
ber of high-pr processes can be generated by proton beam collisions at c.m.s. 
energies of v5 = 200 — 500 GeV, which are sensitive to Ag, such as prompt 
photon and heavy flavor production, jet, single-inclusive and di-hadron pro- 
duction (see e.g. [61]). Indeed, in all cases, the polarized gluon density plays 
a prominent role already at LO of QCD, through the underlying hard pro- 
cesses of gluon-gluon fusion and gluon-quark scattering. 
Complete NLO corrections to polarized QCD hard-scattering processes, such 
as high-pr pion and jet production by polarized proton collisions, are also 
available [6§1 M, EI] • 

In particular, in the NLO analysis [62j of high-p^ pion hadroproduction, the 
predicted asymmetry A^ ] L , i.e. the ratio of polarized and unpolarized cross 
sections for the specific reaction, exhibits sizable differences for two distinct 
sets of spin-dependent pdfs [59] , which both provide good fit to inclusive DIS 
data, but differ significantly in the polarized gluon density (see also [6"3]). 
The gluon polarizations considered in Ref. [62] are the maximal gluon po- 
larization "GRSV-max" - based on the maximal saturation of the positivity 
contraint, i.e. the assumption Ag — g at the input scale - and the "standard" 
distribution of the GRS V analysis [SH] ■ The latter in particular is compared 
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in Fig. l6.20l with the best fit result (type-A) of Sec. 6 and with other available 
parametrizations. As a result, observed spin asymmetries should be capable 
of discriminating among different scenarios. 

First experimental results from PHENIX and STAR detectors at RHIC on 
double helicity asymmetry at c.m.s. energy V5 = 200 GeV, in pion [65] and 
jet [66] production, probing the kinematic region in the parton momentum 
fraction 0.03 < x < 0.3, turn out to be consistent with a moderate gluon 
polarization, such as the "standard" GRSV distribution ^U\. Specifically, 
the observed asymmetries are small, and a large positive gluon contribution 
to the nucleon spin, as given by the maximally polarized gluon "GRSV- 
max" [59], is clearly disfavored by the data [651 [66] . 

Besides, the PHENIX 7r data points have been also included in a recent 
analysis of inclusive DIS world data by the AAC Collaboration [27], and the 
conclusion was reached that the ensuing uncertainty on Ag is significantly 
reduced (by about 60%). Nonethless, within this analysis, along with a pos- 
itive gluon polarization, a negative gluon distribution in the small-x region - 
type-1 and -3 respectively [27] - is still consistent with the data, and resides 
outside the estimated error band of type-1 Ag. Further insights into the 
low-a; shape of the polarized gluon density may be supplied by higher energy 
runs (VS = 500 GeV) at RHIC. 



7.1.2 Lepton- nucleon scattering in polarized fixed tar- 
get experiments 

Polarized lepton-nucleon interactions still play a crucial role for investigating 
the spin structure of the nucleon. In order to obtain observables which are 
sensitive to the polarized gluon density, one has to consider reactions that 
are less inclusive than DIS, in which one measures one or more outgoing 
final-state particles resulting from spin-dependent gluon interactions. The 
largest cross sections typically stem from the photoproduction regime, which 
is characterized by scattered leptons at very small angles, and thus by ex- 
changed photons almost on-shell. 

A crucial point here is the clear-cut sensitivity of cross sections and spin 
asymmetries in photoproduction reactions, such as for example heavy fla- 
vored hadron production, at both collider and fixed target energies, to the 
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shape and size of the spin-dependent gluon distribution, as was shown in 
earlier studies (see e.g. [671 168] ). 

At present, experimental results on lepton-nucleon interactions are supplied 
by low-energy fixed target experiments, where a beam of longitudinally polar- 
ized leptons is scattered off longitudinally polarized nucleon targets at c.m.s 
energies \/S of at most few tens GeV. 

In the framework of polarized DIS, the gluon polarization can be directly ac- 
cessed via the underlying photon-gluon fusion (PGF) mechanism, resulting 
in a quark- ant iquark pair. Experimental signatures to tag this subprocess are 
hadron pairs with high-p-r in the final states, and open-charm events where 
the qq pair is required to be a cc pair and an outgoing charmed meson is 
reconstructed. 

First measurements of longitudinal spin asymmetries in high-p^ hadron-pair 
production are given by HERMES [69J at DESY and more recently by COM- 
PASS [70] at CERN, both in the photoproduction regime, and by SMC [71] 
at CERN in the DIS region (Q 2 > 1 GeV 2 ). However, in such a reaction, the 
measured asymmetries receive contributions not only by pure PGF events, 
but also by a significant fraction of background events, mainly due to the 
two competing processes of gluon radiation by QCD Compton scattering 
(7*g — ► qg) and photon absorption at the lowest order of DIS (7*g — > q). 
High-pr hadron pair photoproduction is affected by large higher order cor- 
rections, and a full perturbative understanding is lacking. Estimates of the 
relative contributions from the background processes then rely on Monte- 
Carlo simulations, which are tuned to describe the data (see e.g. [70]). The 
average gluon polarization is probed within a restricted range of momentum 
fraction, i.e. roughly x ~ 0.1, and turned out, on the whole, to be consistent 
with a moderate gluon polarization (for a summary of such measurements 
see [70]). Phenomenological studies of hadron-pair polarized photoproduc- 
tion at HERMES and COMPASS can be found in [721 GS E3 • 



7.1.3 The open-charm method at COMPASS 

At variance with hadron-pair production, the open-charm approach is free of 
background, since the PGF subprocess is the main mechanism for producing 
charm quarks in polarized DIS. A clean perturbative description is thus possi- 
ble, and higher order corrections have been found to be reasonably small 
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This approach is currently used at COMPASS [75] to extract the ratio of the 
polarized and unpolarized gluon distribution Ag/g. Here, polarized muons 
with a beam energy of E^ = 160 GeV scatter off deuterons in a polarized 
6 LiD target, corresponding to a c.m.s. energy of roughly V5 ~ 18 GeV, in 
the photoproduction regime where the photon virtuality is Q 2 ~ (a descrip- 
tion of the experimental setup can be found in [75]). D mesons in the final 
states, originated from charm quark fragmentation, are reconstructed from 
their decay products in the channel D° — > Kir, with a branching ratio of only 
~ 4%, which results in low statistics. In addition, the channel D* — > D°7r s , 
where n s is a slow pion, is studied, increasing the signal-to-brackground ra- 
tio [77]. 

From the measured spin asymmetry All & preliminary determination of 
Ag/g has been obtained by a LO analysis in QCD, that relies on the ap- 
proximate relation [78] 

All = -Rpgf &ll V ^4bg , (7.1) 

9 

where clll is the partonic asymmetry, i.e. the ratio of spin-dependent and 
spin-averaged hard cross sections, -Rpgf is the fraction of PGF events in the 
selected sample and A^q is the background asymmetry. An estimation of 
(111 on a event-by-event basis is then provided as a function of the measured 
kinematical variables, and is based on the comparison of the data sample 
with Monte-Carlo studies (see [79] for details on the analysis method). 
The extracted value [79] of Ag/g for data collected from 2002 up to 2004 
turned out to be negative, i.e. < Ag/g >= — 0.57 ± 0.41 (stat)±0.17 (syst), 
but still consistent with zero, and represents an average over the probed ar- 
range, < x >~ 0.15. The hard scale is approximatively given by the charm 
quark mass /z 2 = 4(m 2 + p 2 D ) ~ 13 GeV 2 , where p t D is the D° meson trans- 
verse momentum with respect to the photon direction. 
However, the above Monte-Carlo approach entails a number of approxima- 
tions, thus missing a clean theoretical description of the process. On the 
other hand, as already noted, since the photon-gluon fusion is the main un- 
derlying mechanism for producing charm quarks by polarized lepton-nucleon 
scattering, the open-charm reaction can be reliably calculated in perturbation 
theory, with the charm quark mass as the hard scale required. 
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7.2 Theoretical analysis of charm asymme- 
tries at COMPASS 

A thorough analysis of open-charm production from deep inelastic muon scat- 
tering off nucleons, with longitudinally polarized beam and target, has been 
performed here in the photoproduction approximation, where only quasi- 
real photons (Q 2 ~ 0), radiated from the incoming leptons, are considered. 
The analysis is indeed tailored to the COMPASS experiment [751 E2| with 
relation to upcoming results on spin asymmetries in open-charm photopro- 
duction, and is aimed at assessing the sensitivity of the actually observed 
quantities to the polarized gluon distribution. 

As already noted, many earlier studies of heavy quark polarized photoproduc- 
tion, at both high and low energy up to NLO, are also available [67J I68|l80]. 
This issue is here re-analysed by constructing a close theoretical descrip- 
tion of the experimental spin asymmetry at COMPASS, as a function of the 
measured kinematic variables. 

7.2.1 General framework 

The spin-dependent and spin-averaged cross sections are defined by 

Aa = -[a^ -a n ] and a = - [a TT + a n ] (7.2) 

— _ 

respectively, where the arrows denote parallel and opposite helicities of the 
scattering particles. The spin-dependent photoproduction cross section for 
the process IN — > I'HX, where a longitudinally polarized lepton beam scat- 
ters off a longitudinally polarized nucleon target producing one observed 
hadron H in the final state, can be written, on the basis of factorization 
theorem of QCD, as 

dAa tN = y^ / dxj dxj dz k Af> 1 ' (x i: fi F ) Af^ N) (xj, fi F ) (7.3) 

ijk 

x dAaij(xi, Xj, S, P H /z k , (Ir, (J>f, a4) D k( z k, /4) • 

Here, A/- (x„ Hf) is the usual polarized parton density in the nucleon at the 
factorization scale ji F , and the non-pertubative function Dj?(zk, fi' F ) describes 
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the fragmentation process of the produced parton k into the observed hadron 
H at the scale n' F , with momentum fraction Zk- The polarized hard scattering 
cross section dAa^ can be calculated in QCD order-by-order in the strong 
coupling a s (fi R ), /ir being the renormalization scale, and the sum in Eq. 
(J7.3P is taken over all partonic channels. Finally, the function Af { (xi, fip) 
represents the parton density in the lepton. 

In the framework of the photoproduction approximation, which is the phys- 
ically interesting case in connection with the COMPASS experiment, the 
initial-state lepton is considered as a source of quasi-real photons (equivalent- 
photon approximation), with energy distribution given by the Weizsacker- 
Williams spectrum [81] that can be theoretically computed. The function 
Af!> is given by the polarized Weizsacker- Williams function [82] 



(i)( - a f 



A/W(y) 



2vr 



y)2 in g - ( 2 1 r y) +2m t y 



y m W \Qmax m 2 y 2 . 

(7.4) 
where mi is the lepton mass, y the lepton momentum fraction carried by 
the radiated photon and Q ma x is the allowed upper value to the photon's 
virtuality, to be fixed according to the experimental conditions [83]. The 
unpolarized counterpart of Eq. (j7.4p reads 



A l) (y) 



2ty 



i + (i-o)' hl ^(i-v) +2m? / i i-» 



y mfy 2 \Q 2 ma x m fa 2 , 

(7.5) 

Actually, quasi-real photons undergo interactions with partons in the nu- 
cleon either directly ("point-like" component) or via their partonic structure 
("resolved" or "hadronic" component). In the first case photons simply in- 
teract as elementary point-like particles, whereas in the second case photons 
fluctuate into a hadronic state, before the hard QCD interaction takes place, 
and the so-called "resolved" contributions to the cross section compete with 
the "direct" part. Therefore, the experimentally measurable cross section 
is the sum of the point-like and the resolved photon contributions (see for 
example [HU [851 EH] and references therein) 

dAa lN = dAa?£ + dAa r ™ • (7.6) 

Here, the direct component is given by Eqs. (17. 3p and (17.41) . and the resolved 
component can also be cast in the form of Eq. (17.31) . by properly defining the 
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parton density of the lepton as the convolution product 

-1 dy 



Jxi y \y / 



(7.7) 



with A/y (y) given by Eq. (17.41) . The partonic content of the polarized 
photon is unmeasured so far and models are usually invoked [57] in actual 
calculations. Furthermore, the direct photon contribution can be viewed as 
a particular case of Eq. (17.71) . by defining 

A/W =5 (l-^y (7.8) 

The corresponding spin-averaged cross sections are finally given by Eq. (17. 3p 
replacing polarized parton densities and hard cross sections by the related 
unpolarized quantities. The experimentally relevant observable is the spin 
asymmetry, i.e. the ratio of polarized and unpolarized cross sections 

A LL = d -p±. (7.9) 

Focussing on the determination of the gluon polarization from photopro- 
duction reactions, such resolved contributions act as a background. Ex- 
tensive analyses of their relevance with relation to the kinematic regime 
have been performed in studies on high-energy lepton-proton collisions (see 
e.g. [651 |85| [88]). in view of the polarized mode of HERA at DESY and the 
planned eRHIC project at BNL [89]. Lower energies (fixed-target) experi- 
ments are generally expected to be less sensitive to the unmeasured partonic 
content of the polarized photon (see [551158] and references therein), and LO 
estimates [68] have revealed that the resolved contributions are negligibly 
small in the case of polarized photoproduction of charm quarks at COMPASS 
energies. In the analysis of this experiment, one can thus safely neglect the 
resolved piece in Eq. (17. 6p . 

In the case of open-charm photoproduction, as already noted, the only par- 
tonic subprocess contributing in LO to the spin-dependent cross section Eq. 
(17.31) . is photon-gluon fusion 7(7 — » cc. However, in the NLO approxima- 
tion of QCD, along with 1-loop virtual corrections to PGF and real cor- 
rections with an additional gluon in the final state, new subprocesses, in- 
duced by a light quark (antiquark) replacing the gluon in the initial-state, 
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i.e. jq(q) — ► ccq(q), can also contribute, diluting the dependence of the 
asymmetry on Ag(x). 

All NLO contributions have been explicitly calculated in [80], for the polar- 
ized case lj, and found to be reasonably small (around 10%) in the energy 
range accessible at COMPASS (roughly a/S^v — 12GeV), and should be 
outside of the present level of experimental accuracy at COMPASS [75] . 
Furthermore, the background due to the light quark induced subprocesses 
turned out [80] to be fairly negligible. In particular, the PGF mechanism 
dominates in the COMPASS kinematical region, and a clear determination 
of Ag may be performed by detecting open-charm events. 
On this basis, the theoretical analysis of the spin asymmetries in open-charm 
photoproduction is performed below at the leading order of QCD, and by re- 
taining only the point-like photon channel. The calculation is then specialized 
to the COMPASS kinematics in order to reconstruct the actual experimen- 
tally observed quantities. 

7.2.2 Leading order cross section 

Open-charm events in the final states are a clear signature of the LO PGF 
channel, namely 

7 (g, A 7 ) + g{p, X g ) -► c(k) + c(k') , (7.10) 

where, A 7 , q and X g , p denote photon and gluon helicities and four-momenta 
respectively, such that q 2 = p 2 = 0. For the outgoing heavy quark pair k 2 = 
k' 2 = m 2 , m c being the charm quark mass. The corresponding Mandelstam 
variables are then defined by 

s = (q + p) 2 , t\ = (p — k) 2 — m 2 ., and u\ = (q — k) 2 — m 2 , (7.11) 

with s+ti+Ui = 0. To order 0(a em a s ) the differential partonic cross section 
can be written as [HZl E3] 

|^(s, t u u u A 7 , A,) = 2 ™"£^ [E + A 7 A 9 A] 5(s + t 1 + u x ) (7.12) 



lr The NLO QCD corrections in the unpolarized case have been calculated for the first 
time in Refs. [90l [91] . Photoproduction of heavy quarks in the unpolarized case at fixed 
target experiments has been extensively studied [92] . 
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The spin- dependent and spin-averaged hard cross sections are defined as 
usual by 



^9 = l [d&H ~ d&lH , d ^ 9 = I [d*H + d °\ 9 ] 



? L jg "7flj ' -""19 2 L"'"79 ' "' KJ ~19\ V / 

respectively, where the arrows denote parallel and antiparallel polarizations 
of incoming photons and gluon. The Born-level cross sections thus explicitly 
read 
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<J(s + ti + Mi). (7.17) 



Physical cross sections for the leptoproduction process fj,N — > fi'HX where if 
is a charmed meson in the final state, can be related to the photoproduction 
analogue by the Weizsacker- Williams function Eq. (17. 4p . To order 0(a cm a s ) , 
at the partonic level only the gluon distribution is relevant. Then Eq. (17.31) . 
in the unpolarized (polarized) case, can be now rewritten schematically as 



d(A)a 



/iN 



J/max 



dy(A)&\y)d(A)a^ N 



(7.18) 



? M 



where (A)f^'(y) is given by Eq. ([73]) (Eq. (I72j) ). q = yP^ with P M the 
incoming muon four-momentum, and the integration range over y is fixed by 
the experimental conditions. The photoproduction cross section in Eq. (I7.18P 
at LO is given by 



d(A)cr 7A r = / dx(A)g(x,/j 2 F ) / dz D^(z)d{A)a ig {x,z,y,S,T,U^ 2 R ) 



(7.19) 
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where the (photoproduction) hadronic invariants S, T, U are related to their 
partonic counterparts by 

xT U 

s = xS, t\ = — and U\ = — . (7.20) 

z z 

The only parton distribution involved at LO is the unpolarized (polarized) 
gluon density (A)g(x,/M F ), whereas D^{z) is the non-perturbative fragmen- 
tation function of a produced charm quark into an observed charmed meson. 
Using the expressions Eqs. ( 17. 16ft and ( 17.171) of the hard scattering cross 
sections, Eq. (IT. 19[) explicitly reads 



d 2 Aa lN _ 47ra em a s (/i|) f dx 2 f dz H 

■i\g{x,n F ) / — 7T JJ C [z) 



dTdU 9S 2 J x * v '™' J z 

( U xT\ / 2m 2 c z 2 S\ . / n xT U* 



d <t 7 7v 47ra em a s ( / u i j) f dx 



dTdU 9S 2 



x 



g(x,fi 2 F )J^D*(z) x 



U xT Am 2 z 2 S / m 2 z 2 S 
+ Jl^ 1 



5 ( X S + — + - ] , (7.22) 
xT U TU V TU \ « « 

in the spin- dependent and spin-averaged case respectively. 

Actually, if the momentum fraction y of the incoming muon carried by the 
photon is reconstructed event-by-event, the open-charm asymmetries do de- 
pend on y, as well as on the measured variables of the outgoing D meson. 
Within this framework, no average over the Weizsacker- Williams spectrum, 
Eq. (17.181) . is required, and the observed quantity is consistently described 
by the photoproduction asymmetry 

A° LL = ^ (7.23) 

with d(A)a lN given by Eqs. (I7T211 and (I7T22J) . 



7.2.3 COMPASS kinematics 

Specifying the above general discussion to the case of the (fixed target) open- 
charm experiment at COMPASS, the muon momentum fraction y carried by 
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the photon is indeed reconstructed event-by-event with no experimental cut. 
Thus, the polarized and unpolarized cross sections in Eq. ( 17.231) are double 
differential distributions, parametrized by the measured D meson variables 
in the nucleon rest frame, that is, the transverse momentum p t £> with respect 
to the photon direction, and the energy fraction Zd = E D /E 1 with respect 
to the photon energy E T They also depend on y through the incoming 
photon energy, given as a function of the muon beam energy E 1 = yE^ , 
with En = 160 GeV. The hadronic invariant system in the laboratory frame 



then reads 



T 



-Sz 



D- 



U 



s 2 

■ —(z d -xd) and S = 2ME 1 



with 



Xd 



z d~ 



4M 2 m 2 D 
S 2 



(7.24) 



(7.25) 



and rritD — \Iv\d + m Di i- e - ^ ne D meson transverse mass. Using the Ja- 
cobian dTdU = (2p t c,S/xD)dzDdptD to express Eqs. (17.211) and (I7.22p in 
the nucleon rest frame variables Eq. (17.241) . the relevant spin-dependent and 
spin-averaged physical cross sections become 



2^ i\/r2 



d A<t 7 at lGira em a s (fj, )M 



PtD 



dptD dzn 9S 3 

2z D w D 



X 



\Z D + W D 



zdXd(z d -xd) 
Am 2 c M 2 
zd{zd ~Xd)S 2 



/ — Ag(x,fi 2 ){z D + w D ) 

** -^niiii 



[z D + w D ) 



and 



2^ t\/t2 



a ro" T jv lQira cm a s (fx )M 



PtD 



dptD dz D 



9S 3 



i- 2 ^^ q + M ( 



zdXd(z d -xd) 



z D + w D . 



xS 



2m 2 c M 2 



z D (z D - Xd)S 2 
respectively, where the shorthand 



w D 



jzp - Xd)S 
2xM 2 



D (zd + w D ) 



(7.26) 



1 dx 



x 



-g{x,n ) {z D + w D ) x 



,z D + w D 



D (z D + w D ) 



(7.27) 



(7.28) 
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has been used. The actual observable at COMPASS is finally the ratio of 
Eqs. (EH and (172711 . according to Eq. (17231) . 

It is worth noting that in Eqs. (17.261) and (J7.27P integration over the frag- 
mentation variable z has been explicitly performed, being constrained by 
the 5-function. Moreover, the integration range over the gluon momentum 
fraction x depends on the measured variables through the relation 

E-y fz D - xd 



Xmin(PtD, Z D ) = — f - _ 1 ■ (7.29) 

On the other hand, if one explicitly integrates over x, the integration thresh- 
old over z is constrained by the relation 

Zmin(PtD, Z D ) = Z D + ^{ZD ~ Xd) ■ (7.30) 

This actually means that the gluon momentum fraction x and the fragmen- 
tation variable z are not independent variables. 

The strong coupling a s (fi^), as well as the polarized gluon density Ag(x, n 2 F ) 
are evaluated in NLO at the common fixed value f/,R — \ip — 2m c for the 
renormalization and factorization scales respectively. For the charm quark 
mass m c = 1.5 GeV is taken, M ~ 0.938 GeV for the nucleon mass, and the 
D mass is taken as an average of the experimentally detected D° and D* 
masses [Sj, roughly nip ~ 1.94 GeV. The strong coupling is normalized as 
a s (Mz) — 0.118 at the Z boson mass [8]. 

Note that the experimental cut zd ?t 0.2 is imposed on the fractional D 
energy, and the most populated bins are estimated [75] as 0.32 < zd < 0.66 
and 0.35 <p tD <lGeV. 

Finally, as far as Eqs. (17.291) and (17.301) are concerned some comments are in 
order. Indeed, for fixed values of the observed kinematic variables zd, ptD, 
the fragmentation variable z is completely fixed by the momentum fraction 
x, and thus only the integral over x remains as shown by the expression 
of the cross sections Eqs. (17.261) and (17. 27ft . Furthermore, the integration 
threshold over x is determined by the values of zq and ptD- This bound, 
given by Eq. (17.291) . is displayed in Figs. l7.1l and l7T2l as a function of each of 
the two variables. As expected, the low x region (x < 0.1) is unaccesible at 
COMPASS energies, and the integration region involved is stretched towards 
larger x with increasing p t £> (Fig. l7.1|) . The same situation is observed as the 
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D meson energy approaches the kinematical boundaries (Fig. 17.2ft . whereas 
in the physically meaningful region of zd ^min is nearly costant (x m i n ~ 0.1). 
As already noted, one could equivalently determine x for given z and then 
only the integration over z remains. In such a case, it is the fragmentation 
variable z which is bounded from below with z min determined by Zd and p t D 
according to Eq. (17.301) . The bound on z is displayed in Figs. 17.31 and 17.41 
versus z D and p t D respectively. Both figures show that a sizable amount of the 
charm quark momentum is carried by the D meson in the whole kinematical 
region. 



7.2.4 Pdfs and fragmentation function 

The sensitivity of the asymmetry on the spin-dependent gluon distribution 
Ag(x,fj/p) is investigated here by using three (NLO) parametrizations, ob- 
tained from fit to inclusive DIS data, namely AAC [3E], BB [58], and the 
most recent best-fit result BF-A (i.e. type-A) discussed in Sec. 6.5.1. 
The use of NLO parton distributions in a LO calculation has obviously no 
effect on the accuracy of the calculation which remains LO. Note also that 
the AAC and BB polarized gluons are given in the MS factorization scheme, 
whereas the BF-A gluon is given in the Adler-Bardeen scheme. However, 
as already noted, the polarized gluon densities are the same in the two 
schemes at NLO accuracy. The three parametrizations used for Ag(x,/j, F ) 
are displayed at the common scale \xf = 2m c in Fig. I7.5[ together with the 
MRST2004f4 [91] unpolarized gluon distribution, which is used as a rule in 
the spin-averaged cross section Eq. (j7.27p . 

In Fig. 17.61 the LO and NLO MRST2004f4 parametrizations for the unpo- 
larized gluon distribution g{x) are compared in the region of interest for x 
(x > 0.1), at the scale /i = 2m c . This will be used to estimate the dependence 
of the results on the unpolarized gluon distribution. 

Note also that some of the available LO parametrizations of Ag(x) (e.g. BB 
and LSS [28]) violate the positivity constraint imposed by the LO MRST 
determination of g(x) at intermediate and large x . This problem is also 
marginally present for the NLO BB polarized gluon as seen in Fig. I7.5[ but 
does not affect significantly the analysis. Indeed, the uncertainty on the 
unpolarized gluon at large-x makes violation of the positivity bound less 
serious. 
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Figure 7. 1 : Variation of the lower bound x m i n in the integration over the gluon momen- 
tum fraction x, as a function of ptD at three different values of zd- 




ptD=1.GeV 
ptD=0.5GeV 
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Figure 7.2: Variation of the lower bound x m - ln in the integration over the gluon momen- 
tum fraction x, as a function of zd at three different values of PtD- 
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Figure 7.3: Variation of the lower bound z m ; n in the integration over the fragmentation 
variable z, as a function of ptD at three different values of zd- 
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Figure 7.4: Variation of the lower bound z m ; n in the integration over the fragmentation 
variable z, as a function of zd at three different values oiptD- 
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Finally, for the fragmentation function of a charm quark into a D meson (D° 
and D*) the non-pertubative Peterson (P) function [95] 

Df (z) = . (7.31) 

z[l-l/z-e/(l-z)f 

is used as a baseline, with JV the normalization factor. The phenomenological 
parameter e, that is predicted to scale as e ~ AQ CD /mg in the fragmentation 
of a heavy quark Q, is fixed to the usual value e = 0.06 [96] (see also [97]). In 
Fig. 17.71 the Peterson fragmentation function Eq. (17.311) is compared to the 
slightly different parametrization (CN) given in Ref. [98], namely 

D*(z) = N(l-z) ai z a \ a t = 0.8, a 2 = 3.2. (7.32) 

The latter is also used in the present analysis in order to estimate the impact 
of the precise shape of the fragmentation function on the results. 



7.3 Phenomenological results 

The theoretically predicted asymmetries A£ L , given by Eqs. ( 17.231) . ( 1T.26I) 
and (I7.27p . are analysed here as a function of the relevant kinematic variables 
y, Z£>,PtD, so as to determine their sensitivity to the polarized gluon density. 
Next, the effect of including COMPASS results in a global fit of pdfs, along 
with inclusive DIS data, is discussed by using pseudo-data generated on the 
basis of the theoretical prediction and the expected statistical error. 



7.3.1 Asymmetries as a function of the measured vari- 
ables 

The asymmetries A^ L as a function of the incoming photon energy E 1 = yE^ 
are compared in Fig. 17.81 for the aforementioned parametrizations of Ag(x). 
The D meson trasverse momentum and fractional energy are fixed to the 
values ptD = 1 GeV and zd = Ed/E 1 = 0.4 , corresponding roughly to a 
maximum in both the polarized and unpolarized cross sections at the COM- 
PASS kinematics. With increasing photon energy, the asymmetry slowly 
decreases due to the faster growth of the unpolarized cross section, similarly 
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Figure 7.5: NLO paramctrizations of gluon distribution functions used in the analysis, 
namely MRST2004f4 [94] for g(x), AAC [38], BB [58] and the most recent best-fit result 
BF-A for Ag(x), at the scale n 2 F = 4m;?. 
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Figure 7.6: LO and NLO MRST2004f4 parametrizations [94] of the unpolarized gluon 
distribution in the region of interest, at the scale fip = ^m? c . 
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Figure 7.7: Fragmentation functions Df? (z) used in the analysis, namely Eq. (|7.31|) 
(P) [951 [96] and Eq. (TQ2l) (CN) [98] • 



to the case of total charm photoproduction [68]. The gluon polarization has 
a significant impact on the size of the asymmetry, in that the normalization 
is mostly determined by the value of Ag(x) for x > x m [ n ~ 0.1. Because for 
all the given parametrizations Ag(x) is rapidly decreasing in this region, the 
asymmetry is in fact essentially determined by the value of Ag(2 min ). This is 
indeed in substantial agreement with the estimated [79] average x =< 0.15 > 
over the probed region at COMPASS. The shape of the asymmetry is instead 
dominated by that of the underlying hard process. 

The y-shape of A® L is also displayed in Fig. 17.91 at different fixed values of 
Ptn, and for the same polarized gluon density BF-A, used as a benchmark 
in what follows. As the measured transverse momentum p t £> is increased the 
large x region of g(x) and Ag(x) is probed. Because the ratio Ag(x)/g(x) 
increases with x, the asymmetry correspondingly increases with growing p t £>. 

The asymmetry as a function of the fractional energy Zd is given in Fig. 17.101 
at fixed photon energy y = 0.5 and p t D — 1 GeV, for the three different 
parametrizations of Ag(x). Once again the normalization is determined 
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mostly by the value of Ag(x m i n ) and the shape by that of the underlying hard 
process. The whole kinematically allowed region is displayed, even though 
the COMPASS data are cut at z D > 0.2. Note that in the small z D < 0.2 
region x min grows as Zp decreases. Therefore, this kinematical region probes 
Ag at increasingly large values of x, and it is thus sensitive to the poorly 
known polarized gluon at large x. The asymmetry (and even the individual 
cross sections) in this region remains sizable despite the fact that the pdf is 
very small. 

Fig. ETT] displays again the effects of increasing p t D , that is, of sampling only 
events with larger x, with Ag{x) again given by BF-A. 

The ptD distribution is shown in Fig. 17.121 as before for the three parametriza- 
tions of Ag(x), at fixed values of yr> and Zd- The observed p t D, being the 
minimum transverse momentum of the parent charm quark, is directly related 
to the minimum gluon momentum required to initiate the process. Generally 
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Figure 7.8: The spin asymmetry as a function of the photon energy y = E 1 /E ll , at 
PtD = 1 GeV and zjj = Eq/E 1 = 0.4. Three (NLO) parametrizations for the A<? 
are used, namely AAC, BB and BF-A, whereas g is given by MRST2004f4 (NLO). The 
Peterson fragmentation function is used with e = 0.06 . 
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Figure 7.9: The spin asymmetry as a function of y at different p t D and Ag given by 
BF-A (z D , g(x) and Df (z) as in Fig. EH) 



speaking, at low p t D (less than 1 GeV) a larger sample of allowed partonic 
events contributes, but the measurable asymmetry is small. On the other 
hand, the high p t D region (greater than 3 GeV) corresponds to a smaller 
sample of subprocesses restricted to larger x, where the uncertainty in the 
unpolarized gluon distribution dominates. 

Figs. 17.131 and 17.141 show the dependence of the results on the unpolarized 
gluon distribution by comparing the paj-shape of A® L obtained using the LO 
and NLO MRST parametrizations of g(x) and two different polarized gluon 
distributions, AAC and BF-A respectively. In both cases, at LO the positiv- 
ity constraint is violated by the gluon distributions at large x. Specifically, 
the asymmetry exceeds 1 at p t o > 3 GeV if Ag is given by AAC (Fig. 17.131) . 
and at p tD > 4.7 GeV with BF-A (Fig. I7TT4D . This is due to the fact that 
both polarized and unpolarized cross sections become very small in this re- 
gion. Actually, the large difference in the asymmetries in the two cases, that 
is, with the LO and the NLO parametrizations of g(x), in both Figs. 17.131 
and 17.141 indicates that for p t £> > 2 GeV the uncertainty in the unpolarized 
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Figure 7.10: The spin asymmetry as a function of the fractional energy zd = Ejj/E^, 
at fixed photon energy y = 0.5 (ptD, g{%), A.g(x) and D^ (z) as in Fig. \T.. 
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Figure 7.11: The spin asymmetry as a function of zd at different ptD and A<? given by 
BF-A (y, g(x) and Df (z) as in Fig. lTTO)) 
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Figure 7.12: The spin asymmetry as a function of the D meson transverse momentum 
PtD, at fixed photon energy y = 0.5 (zd, <?(#), Ag(x) and Df? (z) as in Fig. [77 



gluon distribution makes the result completely unreliable. 

As an example, the dependence of the asymmetry on the precise shape of 
the fragmentation function is shown as a function of ptD in Fig. 17.151 at fixed 
y and zq. The BF-A and NLO MRST parametrizations are used for the 
polarized and unpolarized gluon respectively. The choice of the fragmenta- 
tion function has no significant impact on the ^^-distribution over the entire 
kinematical region. Similar results are found for the asymmetries as a func- 
tion of the other relevant variables. 

Finally, as already noted in [68], also a moderate scale-dependence of the 
asymmetry as a function of p t o is observed in the region of interest, i.e. 
for p tD < 2 GeV. This is displayed in Fig. 17.161 with Ag(x) and g(x) as in 
Fig. 17.151 for a simultaneous variation of the factorization and renormaliza- 
tion scales by a factor a = 0.25, 1, 4 around the value \j? f = [i 2 R = 4m;?. 
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Figure 7.13: Effect of the uncertainty in the unpolarizcd gluon distribution at large x 
on the asymmetry as a function of ptD- LO and NLO MRST parametrizations of g are 
compared, with Ag given by AAC and D^(z), y, zd as in Fig. 17.121 




Figure 7.14: Effect of the uncertainty in the unpolarized gluon distribution g as in 
Fig. I7TT51 with Ag given by BF-A. 
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Figure 7.15: The asymmetry as a function of p t D with two different fragmentation 
functions, P and CN. Ag, y and zd are as in Fig. 17.141 and g is given by MRST at NLO. 
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Figure 7.16: The asymmetry as a function of ptD at different scales fip = ji 2 R = a 4m 2 , 
a = 0.25, 1,4. Ag, g, y and zd arc as in Fig. 17.141 and D^ is given by the Peterson 
function. 
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7.3.2 Results from pseudo-data fit 

As a last step, pseudo-data for open-charm asymmetries A® L have been gen- 
erated by a Gaussian-distributed random shift around the theoretical values. 
The latter are based upon two different assumptions on the polarized gluon 
density, namely the AAC parametrization [38] and the best-fit result of type- 
A (BF-A) discussed in Sec. 6. Both are displayed in Fig. 17.51 at the scale of 
the process at hand. 

The primary goal is to determine to what extent such data can pin down 
the x-shape of Ag, so far unconstrained by inclusive DIS measurements, and 
potentially reduce the systematic error induced by the choice of the initial 
parametrization. 

To this end, twelve bins have been selected for the measured kinematic 
variables, with central values y = 0.35, 0.65, z D = 0.45, 0.65 and p t D = 
0.5, 1, 2 GeV. These are chosen such that they uniformly cover the experi- 
mentally relevant kinematical region. 

The related asymmetries A^ L are calculated in each bin, for the chosen po- 
larized gluon. On the basis of the total number of the reconstructed D 
mesons and their distribution with respect to each relevant variable [79] , 
a statistical error is assigned to the predicted asymmetries, i.e. SA^ L = 
■\/(l — (A® L ) 2 ) / N , where N is the estimated total number of events in each 
bin. The theoretical values are then randomly shifted with Gaussian distri- 
bution to produce the final set of 12 pseudo-data for the assumed polarized 
gluon parametrization. 

As expected, due to the sensitivity of A® L to Ag, very different values of 
pseudo-asymmetries are found by using the two distinct polarized gluons, 
i.e. AAC and BF-A; the two sets of related pseudo-data, when included in 
the global fit of pdfs along with inclusive DIS data (see Sec. 6), also produce 
very different results for the fitted gluon polarizations. 

Indeed, a first set of 12 pseudo-data has been generated by using the AAC 
parametrization for Ag(x) (Aaac i n what follows), and included in the global 
fit of type-A (Sec. 6). A glance to Fig. 17.171 reveals the significant impact 
of the Aaac pseudo-data on the x-shape of Ag. Here, the original best fit 
result BF-A (only inclusive data) for Ag at the initial scale Q 2 = 1 GeV 2 is 
compared with that obtained by including Aaac pseudo-data. 
Also shown in Fig. 17. 171 are the statistical errors propagated to characteristic 
points of xAg(x), namely x — 0.1 (roughly corresponding to the region 
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Figure 7.17: Best fit result (type-A) for Ag at Q 2 = 1 GeV 2 including ^aac pseudo- 
data, compared to BF-A (only inclusive data). 



probed by COMPASS), the position and the value of the peak in xAg(x). As 
can be seen, in particular the peak values in the two cases are not consistent 
within several standard deviations. 

Then, in order to perform a quantitative study of the impact of open-charm 
asymmetries on the x-shape of Ag, a second set of pseudo-data has been gen- 
erated by using the best-fit result of type-A for the polarized gluon (Abf-a)- 
^4bf-a pseudo-data are then included in two types of fit to world data, type- 
A and B respectively (see Sec. 6.5.2), that differ by the assumed functional 
form of pdfs at the initial scale. The results for the gluon polarization are 
summarized in Fig. 17.181 

As already pointed out, the two original best-fit gluons, type-A and B re- 
spectively, both provide equally good fit to inclusive DIS data, even though 
they significantly differ in shape. These are represented in Fig. I7.18l by black- 
solid (BF-A) and black-dashed (BF-B) curves respectively. The red-solid and 
red-dashed curves in Fig. 17.181 denote the type-A and B results respectively 
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(CH-A and CH-B), obtained by adding the Abf-a pseudo-data to the fit. 
As a measure of the systematic error induced by the choice of the input gluon 
density, one may compare the values of characteristic points (and the related 
errors) of the x-shape for the two types of best-fit gluons (type- A and B), 
with and without charm data. Specifically, for only inclusive data fits, the 
position of the peak of xAg(x) and the corresponding function value in the 
case of type-A parametrization turn out to be 



a^max = 0.122 ± 0.071, x max Ag(x r _ 
whereas for the type-B fit they read 
x m ax = 0.552 ± 0.169 x raa , x Ag(x m 



0.156 ±0.062 (BF-A), (7.33) 



0.421 ± 0.871 (BF-B) . (7.34) 



The quoted uncertainties are obtained by propagating the related statistical 
errors of the parameters from the fit. 




Figure 7.18: Best fit results for Ag at Q 2 = 1 GeV 2 for inclusive DIS data alone (black 
curves) and including ^4bf-a pseudo-data (red curves). 
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On the other hand, by including charm pseudo-data Abf-a, one has 

aw = 0.104 ±0.041, x max Ag(x max ) = 0.151 ± 0.047 (CH-A) (7.35) 

for the type-A parametrization, and 

x max = 0.057 ±0.132, x max As(x max ) = 0.100 ±0.118 (CH-B), (7.36) 

in the case of type-B, the main effect being obviously on the latter. 
Specifically, one may take e.g. the half- difference of the central values in the 
two fits, A and B respectively, as an estimate of the uncertainty on the posi- 
tion of the peak in the polarized gluon density; in the case of inclusive DIS 
data alone, from Eqs. (I7.33H7.34)) this amounts to roughly ~ 0.2, whereas 
by including charm pseudo-data in the two fits it turns out to be strongly 
reduced, i.e. ~ 0.02 by using Eqs. (I7.35H7.361) . Indeed, in the latter case, the 
positions of the peak of xAg(x) for type-A and B are now largely consistent 
within the errors. Similar consideration hold also for the corresponding func- 
tion values. 

This fact clearly demonstrates the capability of open-charm asymmetries in 
constraining the x-shape of the gluon polarization, significantly reducing the 
bias induced by the choice of input densities. 

Finally, the two best-fit results for all the pdfs parameters, including Abf-a 
pseudo-data are shown in Tab. 9 below. Comparing these results with the 
type-A and B fits to inclusive DIS data alone, respectively Tabs. 5 and 7, 
it can be seen that, as expected, charm pseudo-data have no sizable effect 
on quark distributions. Moreover, the statistical errors are generally slightly 
reduced, and the quality of the fit improved. 

In particular, as far as the gluon first moment is concerned, very close values 
arise from the two types of fits if charm pseudo-data are included. Averaging 
over the two best-fit values one would obtain Ag(l, 1 GeV 2 ) = 0.54 ± 0.08, 
to be compared with the second of Eqs. (16.91) . 

One may thus conclude that the much constrained x-shape of the spin- 
dependent gluon distribution, beside the more precise determination of the 
first moment, is the most relevant feature of a global analysis that includes 
open-charm asymmetries along with inclusive DIS data. Indeed, charm 
pseudo-data strongly reduce the impact of the initial parametrization on the 
results, that significantly biases inclusive data analyses. Thus, on the basis of 
the above results, real charm asymmetries from the COMPASS experiment 
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should be capable to discriminate among the existing different best-fit forms 
of the gluon polarization which equally well describe DIS data. 



Parameters (Q% = 1 GeV 2 ) 


type- A (CH-A) 


type-B (CH-B) 


Vx 


0.388 ± 0.018 


0.391 ± 0.016 


«£ 


1.184 ± 0.317 


2.670 ± 0.412 


0E 


1.958 ± 0.974 


3.614 ± 0.456 


7e 


-1.229 ± 0.198 


- 


Vg 


0.525 ± 0.123 


0.553 ± 0.095 


a g 


-0.481 ± 0.321 


3.529 ± 0.257 


& 


8. (fixed) 


2.281±1.360 


lg 


6.620 ± 14.070 


11.656 ± 12.533 


V3 


1.127 ± 0.034 


1.207 ± 0.034 


OiNS 


-0.315 ± 0.270 


1.624 ± 0.110 


Pns 


2.965 ± 0.328 


5.498 ± 0.170 


Ins 


8.966 ± 12.571 


-0.180 ± 0.050 


X 2 /d.o.f 


0.888 


0.902 



Tab. 9. Fits of type-A and B, including 12 charm pscudo-asymmctries Abf-a to the set 
of 238 inclusive DIS data points (cut at W > 3 GeV, see Sec. 6). 



Chapter 8 



Conclusions 



To summarize, the possibility of inferring the spin-dependent gluon distribu- 
tion from scaling violations has been thoroughly studied, by means of a full 
NLO phenomenological analysis of the available experimental inclusive DIS 
data, performed in the Adler-Bardeen factorization scheme. 

Special attention, in particular, has been paid on the selection of the lower- 
energy data, in order to avoid systematic errors induced by higher-twist cor- 
rections. Indeed, a bulk of data in polarized DIS are restricted to the region 
of the final-state invariant mass which involves large Bjorken-x at moderate 
values of the momentum transfer Q 2 , roughly 1 < W < 3 GeV. 
It has been shown that such data do not sizably affect the shape of quark and 
gluon densities, whereas the main effect of the lower energy points (mainly 
from the CLAS and the HERMES experiments) is essentially the improve- 
ment in the determination of the higher-twist contributions to the structure 
function g\. Then, for a consistent perturbative treatment of world data, 
lower energy points can be safely discarded; to this end, a quantitative cri- 
terion has been established by a phenomenological study of higher-twist cor- 
rections to g\. As a result, a lower bound has been imposed on DIS data at 
W — 3 GeV, since above this threshold higher-twist corrections are consistent 
with zero within 1.5 standard deviations, thus justifying a pure perturbative 
analysis. 

Furthermore, the bias produced by the choice of input densities has been 
investigated by performing global fits to the same subset of DIS data using 
different boundary conditions for pdfs (namely type-A and B fits). 
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It turns out that the spin-dependent gluon distribution exhibits a sizable 
dependence on the assumed functional form, whereas minor effects are ob- 
served in the quark sector. Indeed, best-fit gluons obtained with different 
input densities (type-A and B) both provide equally good fits to inclusive 
DIS data, even though they significantly differ in shape. On the other hand, 
the first moment of the polarized gluon turns out to be more stable with 
respect to the initial parametrization, and, in particular, smaller if compared 
to previous estimates. If all the theoretical errors are properly taken into 
account, the final estimate reads 

Ag(l, 1 GeV 2 ) = 0.68 ± 0.12 (exp) ± 0.29 (th), 

in reasonable agreement with the results of other recent phenomenological 

studies. 

Thus, even though the better coverage of the (x, Q 2 )-plane, mainly due to 

the recent precise COMPASS data, allows a more reliable determination of 

the gluon first moment, an analyisis based on inclusive DIS data alone is 

not capable to pin down the x-shape of the gluon distribution, which is still 

largely unknown. 

At present, further constraints are expected to come from direct measure- 
ments of Ag(x), which are based upon the observation of gluon- initiated pro- 
cesses at the partonic level, such as for example production of heavy-flavored 
hadrons and jets with large transverse momentum, from both hadro- and 
lepto-production. 

Actually, one of the most promising approach is the measurement of open- 
charm events from (fixed target) polarized lepton-nucleon scattering, since 
charmed mesons in the final states are a clean signature of the underlying 
photon-gluon fusion mechanism, with essentially no competing processes and 
moderate NLO corrections. Furthermore, asymmetry data of open-charm 
photoproduction from the COMPASS experiment are upcoming. Therefore, 
a LO analysis of this process has been performed by constructing a close 
theoretical description of the actually observed quantities at COMPASS. 

As a first result, this study has revealed in particular that, in the range 
of gluon momentum fractions accessible at COMPASS energies, that is re- 
stricted to x > 0.1, the absolute value of the gluon polarization has a signif- 
icant impact on the size of the asymmetry viewed as a function of each of 
the kinematic variables of the outgoing D meson. This is interesting in view 
of the fact that available polarized gluon distributions have similar values 
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of the first moment, but very different shapes and thus very different values 
at fixed x > 0.1. For instance, the best-fit result of type-A of the present 
analysis, BF-A, and the AAC parametrizations both have a first moment of 
Ag around 0.5 at Q 2 = 1 GeV 2 , but at x — 0.2 they differ by more than a 
factor two, as seen in Fig. 17.51 Open-charm data may thus resolve the gluon 
x-shape in this region. 

In order to assess to what extent such data can actually pin down the re- 
shape of Ag(x), so far unconstrained by scaling violations, and potentially 
reduce the systematic error induced by the choice of input densities, a phe- 
nomenological study has been performed by using open-charm pseudo-data, 
generated by a Gaussian-distributed random shift around the theoretical val- 
ues. Adding pseudo-asymmetries to global fits of world data with different 
input densities, type-A and B respectively, a much constrained x-shape for 
the fitted gluon turns out, as well as a better determined first moment. 
Specifically, the two best-fit forms for the gluon polarization become largely 
consistent in shape if charm pseudo-data are included in the analysis. 
This fact clearly emphasizes that, at variance with analyses based on scal- 
ing violations, the open-charm approach should be capable to significantly 
reduce the bias induced by the choice of the initial parametrization, which 
is one of the most relevant theoretical uncertainties in the determination of 
the spin-dependent gluon distribution, and thus to discriminate among the 
available different best-fit forms which equally well describe DIS data. 
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